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_ 1_ 

Executive  Summary 


This  report  presents  a  model  of  manganese  and  iron  flux  from  sediments.  The  reaction 
kinetics  for  iron  and  manganese  in  the  overlying  water  column  are  also  developed.  As 
shown  below,  the  chemistry  of  iron  and  manganese  are  sufficiently  similar  so  that  a 
common  model  formulation  can  be  used.  The  model  is  an  extension  of  the  nutrient 
and  oxygen  sediment  flux  model  developed  previously  for  application  to  Chesapeake 
Bay.  It  has  subsequently  been  appUed  to  freshwater  lake  sediments  as  well.  Thus  the 
model  developed  in  this  report  is  a  comprehensive  model  for  nutrients,  oxygen  and 
sulfide,  and  now  iron  and  manganese. 

The  model  is  calibrated  initially  using  two  data  sets  which  contained  measured 
sediment  manganese  fluxes  together  with  nutrient  and  oxygen  fluxes.  The  detailed 
calibration  used  the  results  from  an  extensive  three  year  nutrient  addition  experiment 
conducted  using  the  MERL  mesocosm  experimental  facility.  This  investigation  in¬ 
cluded  complete  water  column  and  sediment  flux  sampling.  No  comparable  iron  flux 
data  sets  could  be  located,  presumably  because  aerobic  sediment  iron  fluxes  are  small. 

The  model  is  calibrated  initially  in  stand  alone  mode,  in  which  the  overlying  water 
concentrations  and  the  depositional  fluxes  of  organic  matter  and  particulate  manganese 
are  specified.  The  final  stand  alone  calibration  is  made  with  the  model  coupled  to  an 
interactive  water  column  volume,  in  which  the  kinetics  of  manganese  oxidation  are 
explicitly  included.  The  interaction  of  the  water  column  and  sediment  turns  out  to  be 
a  critical  component  in  determining  the  magnitude  and  seasonal  distribution  of  the 
sediment  flux. 

The  iron  and  manganese  sediment  flux  model  is  combined  with  a  water  column 
eutrophication  model  and  apphed  to  two  lakes.  Iron  and  manganese  state  variables 
and  the  appropriate  reaction  kinetics  are  added  to  the  eutrophication  model.  The 
simulations  reproduce  the  annual  cycle  of  algal  growth  and  stratification,  summer 
anoxia  and  the  increase  of  iron  and  manganese  in  the  hypolimnion,  and  fall  overturn 
and  the  return  to  low  metals  concentrations. 

The  sediment  model  required  an  adjustment  to  the  anaerobic  layer  partition  coef¬ 
ficient  for  iron,  in  order  to  match  the  data  from  the  two  lakes.  No  other  iron  model 
parameter  needed  changing.  In  contrast,  the  manganese  model  parameters  were  un¬ 
changed  from  the  stand  alone  cahbration  values.  This  suggests  that  a  possible  exten¬ 
sion  to  the  model  framework  is  to  replace  the  empirical  partition  coefficient  for  iron 
with  a  formulation  that  exphcitly  calculates  the  pore  water  -  soHd  phase  partitioning 
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2  EXECUTIVE  SUMMARY 


using  a  chemical  model.  The  model  already  computes  the  iron  and  sulfide  concen¬ 
trations.  The  partitioning  formulation  would  be  replaced  by  an  explicit  computation 
of  the  formation  of  iron  sulfide.  An  initial  attempt  at  such  a  formulation  is  reported 
below  for  the  formation  of  manganese  carbonate. 


_ 2_ 

Chemistry  of  Manganese  and  Iron 


2.1  MANGANESE 

The  chemistry  of  manganese  in  surface  waters  and  sediments  has  been  studied  for  quite 
some  time  [Stumm  and  Morgan,  1970].  Manganese  exists  in  two  valances  states:  the 
+4  state  Mn{IV)  in  oxic  waters  and  the  +2  state  Mn{II)  in  anoxic  waters.  Mn{IV) 
is  very  insoluble  and  precipitates  to  form  manganese  oxide,  Mn02{s): 

Mn{IV)^+  +  2H2O  Mn02{s)  +  4fr+  (2.1) 

which  is  the  predominant  form  of  manganese  in  oxic  surface  waters.  It  usually  exists  as 
a  coating  on  particles  [Jeime,  1968].  As  the  particles  settle  to  the  sediment,  manganese 
is  transported  as  well  providing  a  source  of  manganese  to  the  sediments. 

In  the  oxic  layer  of  the  sediment,  Mn02{s)  is  stable.  However,  particle  Tniving 
causes  manganese  containing  particles  to  be  transported  to  the  anaerobic  layer  of 
the  sediment  where  manganese  oxide  is  thermodynamically  unstable  and  a  reduction 
reaction  occurs.  Mn{IV)  is  reduced  to  Mn{JJ).  For  this  to  occur,  two  electrons  are 
required  as  shown  by  the  reduction  half  reaction: 

Mn02{s)  +  2e-  ^  Mn{Ilf+  -  4H'+  +  2H2O  (2.2) 

The  primary  source  of  electrons  in  sediments  is  decaying  organic  matter  CH2O  and 
the  oxidation  half  reaction  is: 

CH2O  +  H2O  CO2  +  41?+  +  4e“  (2.3) 

The  overall  reduction  reaction  can  be  constructed  by  supplying  the  electrons  required 
in  eq.(2.2)  from  eq.(2.3)  to  form  the  redox  reaction: 

Mn02{s)  +  ^CH20  +  2H+  ^  Mn{Ilf+  +  +  |i?20  (2.4) 

In  contrast  to  Mn{IV),  Mn{II)  is  more  soluble  and  exists  in  the  mg/L  range  in 
sediment  pore  waters.  As  a  consequence  it  can  diffuse  to  the  oxic  layer  of  the  sediment 
where  it  is  subject  to  oxidation.  The  oxidation  of  Mn{II)  to  Mn{IV)  occurs  via  the 
loss  of  two  electrons: 

(2.5) 


Mn(JJ)2+  ^  Mn{IV)'^+  +  2e- 
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For  oxygen  as  the  electron  acceptor,  the  overall  reaction  can  be  found  using  the  re¬ 
duction  half  reaction  for  oxygen: 

O2  +  4H+  -t-  4e-  2H2O  (2.6) 

SO  that: 

Mn{Ilf+  +  ^02  +  ^  Mn(IV)^+  -f  H2O  (2.7) 

The  Mn{IVY'^  that  is  formed  precipitates  as  manganese  oxide,  eq.(2.1),  and  the 
overall  redox  reaction  is: 

Mn{II)^+  +  ^02  +  H2O  Mn02{s)  -f  2H'^  (2.8) 

This  is  the  reaction  that  occurs  in  the  aerobic  layer.  The  kinetics  of  this  reaction  have 
been  examined  [Morgan,  1967]  and  found  to  be  slow  in  the  normal  pH  ranges  of  surface 
waters  (Fig.2.1).  However  the  reaction  can  be  bacterially  mediated  and  proceed  more 
rapidly  [Jaquet  et  al.,  1982;  Dortch  and  Hamlin  Tillman,  1995]. 


Fig.  2.1  Oxidation  of  Mn(II).  Semi-logarithmic  plots  of  normalized  manganese  concentration 
versus  time 


2.1.1  Solubility 

Fig.2.2  presents  a  solubility  diagram  for  manganese  as  a  function  of  pH  -f-  pe  where 
pe  is  the  activity  of  the  electron  which  is  a  measure  of  the  redox  potential  [Lindsay, 
1979].  The  pH  of  sediment  pore  water  is  typically  pH  =  6-8.  The  oxic  layer  has  a 
pe  ~  13  in  the  presence  of  dissolved  oxygen,  whereas  for  the  anoxic  layer  pe  ~  —3,  the 
redox  potential  at  which  sulfate  reduction  takes  place. 

For  the  oxic  conditions  pH  -f  pe  ~  20,  and  manganese  oxide,  Mn02(s),  is  the 
predominant  species.  The  concentration  of  manganese  is  predicted  to  be  very  low  ( 
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pe  ♦  pH 

Fig.  2.2  Solubility  of  manganese  minerals  at  pH  =  7.  The  insert  shows  the  effect  of  varying 
pH  [Lindsay,  1979] 


<  see  Fig.2.2).  Thus  there  should  be  no  dissolved  manganese  in  the  oxic  layer 

at  equihbrium.  We  shall  see  below  that  the  kinetics  of  the  formation  of  manganese 
oxide,  Mn02{s),  controls  the  dissolved  manganese  concentration. 

For  anoxic  conditions,  pH  +  pe  <  10  and  the  concentration  is  probably  controlled 
by  the  solubility  of  manganese  carbonate  MnCOz{s)  (rhodochrosite).  It  is  unlikely 
that  manganese  sulfide  MnS{s)  is  present  since  iron  sulfide  FeS{s)  is  present  in  most 
sediments  and  it  is  more  insoluble  than  MnS{s)  [Emerson  et  al.,  1983].  Typical  pore 
water  manganese  concentrations  are  in  the  mg/L  range  or  10“"^  M  which  is  within  the 
stability  field  of  MnCOz  (s)  in  Fig.2.2. 

In  addition,  Mn  (//)  will  partition  to  sorption  sites  on  the  sediment  particles.  Thus 
not  all  of  the  Mn  (II)  that  is  formed  by  the  reduction  of  Mn02(s)  remains  in  dissolved 
form  even  if  manganese  carbonate  does  not  form.  Some  will  sorb  to  sediment  particles. 
Therefore  the  transfer  of  Mn  (II)  from  the  anoxic  to  the  oxic  layer  occurs  via  particle 
mixing  which  transports  particulate  Mn,  and  via  the  diffusion  of  soluble  Mn(II). 
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Fig.  2.3  Solubility  of  iron  as  a  function  of  pH  +p€  for  pH  =  7  [Lindsay,  1979]. 


2.2  IRON 

Like  manganese,  the  chemistry  of  iron  in  natural  waters  and  sediments  has  also  been 
extensively  studied  [Stumm  and  Morgan,  1970).  Iron  exists  in  two  valances  states:  the 
+3  state  Fe{III)  in  oxic  waters  and  the  +2  state  Fe{II)  in  anoxic  waters.  Fe{III)  is 
very  insoluble  and  forms  iron  oxyhydroxide,  FeOOH{s),  and  eventually  more  insolu¬ 
ble  iron  oxides.  Fig.2.3  presents  a  solubility  diagram  for  iron  as  a  function  of  pH  +  pe 
[Lindsay,  1979]  At  equilibrium  iron  is  more  insoluble  in  oxic  waters  than  manganese 
(Fig.  2.2).  However  iron  carbonate  (Siderite)  and  manganese  carbonate  have  approxi¬ 
mately  the  same  solubility.  Thus  pore  water  iron  concentrations  should  be  comparable 
to  manganese. 

Like  manganese,  iron  usually  exists  as  a  coating  on  particles.  However,  unlike 
manganese,  there  are  other  forms  of  particulate  iron  that  exist  in  natural  waters. 
Since  the  earth’s  crust  is  approximately  2%  iron  by  weight  [Lindsay,  1979],  particles 
that  run  off  into  natural  waters  contain  a  large  amount  of  iron.  As  the  particles  settle 
to  the  sediment,  iron  is  transported  as  well.  This  is  the  source  of  iron  to  the  sediments. 
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Not  all  iron  in  sediments  is  reactive.  It  is  convenient  to  denote  the  reactive  portion 
of  oxic  iron  as  F eOOH{s)  (Goethite),  and  to  assume  that  it  includes  iron  hydroxide, 
Fe{OH)z{s)  as  well  since  Fe{OH)^  differs  from  FeOOH  (s)  in  water  content  only: 

FeOOH  (s)  +  H2O  Fe  (OH)^  (s)  (2.9) 

The  term  iron  oxyhydroxide  is  meant  to  denote  the  sum  of  FeOOH{s),  Fe{OH)z{s), 
and  any  other  reactive  iron  solid  phases  that  are  present  [Canfield,  1989]. 

In  the  oxic  layer  of  the  sediment,  FeOOH{s)  is  stable.  However,  as  particle  mixing 
transports  iron  containing  particles  to  the  anaerobic  layer  of  the  sediment,  iron  oxy¬ 
hydroxide  is  thermodynamically  unstable  and  a  reduction  reaction  occurs.  Fe{III) 
is  reduced  to  F e(/J).  For  this  to  occur,  one  electron  is  required  as  shown  by  the 
following  reduction  half  reaction: 

FeOOH{s)  -I-  e"  Fe{Ilf+  -  +  2H2O  (2.10) 

Again  the  source  of  electrons  in  sediments  is  organic  matter,  CH2O  (eq.2.3)  and  the 
overall  reduction  reaction  can  be  written: 

FeOOH(s)  +  ^CH20  -I-  2H+  Fe{Ilf+  +  ^C02  +  ^i?20  (2.11) 

By  contrast  to  F e(III),  Fe(II)  is  more  soluble  and  exists  in  the  low  mg/L  range  in 
sediment  pore  waters.  As  a  consequence  it  can  diffuse  to  the  oxic  layer  of  the  sediment 
where  it  is  subject  to  oxidation.  The  oxidation  of  Fe(II)  to  Fe(III)  occurs  via  the 
loss  of  one  electron: 

Fe(//)2+  -  e-  ^  Fe{IIlf+  (2.12) 

With  oxygen  as  the  electron  acceptor,  the  overall  reaction  can  be  found  using  the  half 
reaction  for  oxygen,  eq.(2.6),  followed  by  the  precipitation  of  iron  oxyhydroxide: 

Fe{IIlf+  +  2H2O  FeOOH{s)  +  ZH+  (2.13) 

to  yield  the  overall  redox  reaction: 

Fe(//)2+  -f-  io2  +  ^H20  FeOOH(s)  -f  2H+  (2.14) 

This  is  the  reaction  that  occurs  in  the  aerobic  layer.  The  kinetics  of  this  reaction 
have  been  examined  and  found  to  be  rapid  in  the  normal  pH  ranges  of  surface  waters 
(Fig.2.4).Thus  the  chemistry  of  manganese  and  iron  are  quite  similar.  The  oxidized 
forms  are  both  insoluble  and  form  oxides.  The  reduced  forms  are  soluble  in  the  mg/L 
range.  Their  concentrations  in  pore  water  are  regulated  by  solid  phases  and  sorption 
to  particles.  Their  flux  to  the  overlying  water  is  controlled  by  the  rate  at  which  the 
reduced  forms  are  oxidized  in  the  aerobic  layer.  A  rapid  oxidation  rate  prevents  escape 
since  insoluble  particles  form.  A  reduced  oxidation  rate,  or  anoxic  conditions  in  the 
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Fig.  2.4  Oxidation  of  Fe(II).  Semi-logarithmic  plot  of  the  normalized  iron  concentration  versus 
time. 

overlying  water,  allows  dissolved  metal  to  escape  as  fluxes  to  the  overlying  water. 
Responses  to  lowered  dissolved  oxygen  appear  to  be  similar  [Sundby  et  al.,  1986]. 
Hence  it  is  expected  that  the  same  framework  for  sediment  flux  models  can  be  used 
for  both  manganese  and  iron.  This  is  pursued  in  the  next  chapter. 
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Steady  State  Partitioning  Model 


3.1  INTRODUCTION 

The  chemistry  of  manganese  discussed  in  section  2.1.  The  reactions  that  occur  in  the 
aerobic  and  anaerobic  layer  suggest  a  conceptual  model  for  manganese  fluxes  from 
sediments  that  is  illustrated  in  Fig.  3.1.  This  type  of  model  has  been  suggested  by  a 
number  of  investigators  [Sundby  et  al.,  1986].It  is  instructive  to  think  of  the  processes 
in  temporal  sequence,  corresponding  to  the  numbers  in  Fig.  3.1. 

1.  Particulate  manganese  dioxide  Mn02(s)settles  to  the  aerobic  layer  of  the  sedi¬ 
ment. 

2.  Particle  mixing  moves  the  particle  downward  into  the  anaerobic  layer  of  the 
sediment. 

3.  Manganese  dioxide  is  imstable  in  a  reducing  environment  so  it  is  reduced  to 
soluble  manganese  (II)  Mn{JI). 

4.  Mn  {II)  diffuses  to  the  aerobic  layer. 

5.  Oxidation  of  Mn  {II)  to  manganese  dioxide  occurs  in  the  aerobic  layer  which 
regenerates  manganese  dioxide  Mn02{s). 

6.  However,  the  oxidation  reaction  converting  Mn{II)  to  Mn02{s)  is  competing 
with  the  diffusion  of  Mn  {II)  to  the  overlying  water. 

7.  If  the  oxidation  of  Mn{II)  is  fast  relative  to  its  diffusion  to  the  overl3dng  water, 
then  the  only  other  sink  is  via  burial  by  sedimentation.  If,  on  the  other  hand, 
the  oxidation  of  Mn  {II)  is  slow,  then  Mn  {II)  escapes  to  the  overlying  water 
as  a  flux  of  soluble  manganese,  completing  the  cycle. 

This  conceptual  model  description  is  somewhat  idealized  since  the  processes  all 
occur  simultaneously.  However,  the  main  point  is  that  the  magnitude  of  the  flux  to 
the  overl5dng  water  is  due  to  a  competition  between  the  rate  of  oxidation  of  Mn  {II) 
to  Mn02{s)  in  the  aerobic  layer,  and  the  diffusion  of  Mn  {II)  from  the  aerobic  layer 
to  the  overlying  water. 
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Fig.  3.1  Conceptual  model  for  the  flux  of  manganese  from  sediment. 


3.1.1  Partitioning 

In  addition  to  the  redox  reactions  discussed  above,  the  chemistry  of  manganese  (II) 
itself  is  also  important.  Since  it  is  a  divalent  cation  it  can  adsorb  to  sorption  sites  that 
are  present  in  the  sediment.  In  particular  it  can  sorb  to  the  freshly  precipitated  hydrous 
iron  oxide  [Dzombak  and  Morel,  1990]  in  the  aerobic  layer  and  to  other  sorption  sites 
in  the  anaerobic  layer. 

Manganese  (II)  can  also  form  solid  species  (see  Fig.  2.2).  A  model  formulation  which 
explicitly  considers  the  solubility  of  manganese  carbonate  is  presented  in  Chapter  4. 
As  we  shall  see,  the  result  is  a  rather  formidable  set  of  equations. 

For  the  model  presented  below,  the  formation  of  particulate  manganese  is  modeled 
as  a  reversible  partitioning  reaction 

Mn(//)^^Mn(//)p 


(3.1) 
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where  Mn  (II)^  and  Mn  (II)p  are  dissolved  and  particulate  manganese  (II)  respec¬ 
tively.  The  rational  for  this  choice  is  the  same  as  that  employed  for  the  nutrients 
and  oxygen  flux  model  [Di  Toro  and  Fitzpatrick,  1993].  First,  the  resulting  equations 
can  be  solved  analytically  which  is  an  important  aid  to  understanding  the  model’s 
behavior.  Second,  linear  partitioning  can  sometimes  be  a  realistic  description  of  the 
relationship  between  dissolved  and  particulate  chemical.  It  is  the  limiting  case  of  the 
Langmuir  model  which  is  implicit  in  most  surface  complexation  model  formulations 
[Westall  and  Hohl,  1980;  Dzombak  and  Morel,  1990].  Finally,  the  general  problem  of 
computing  the  chemical  composition  of  pore  water  would  involve  using  a  numerical 
chemical  equilibrium  model  [Jahnke  et  al.,  1994].  Mass  balance  equations  are  required 
for  the  various  chenucals  that  affect  the  pore  water  chemistry — ^for  example  hydro¬ 
gen  ion,  carbon  dioxide,  and  so  on  (e.g.  [Di  Toro,  1976]).  Thermodynamic  data  are 
required  for  the  relevant  aqueous  complexes,  and  the  stable  and  metastable  mineral 
phases,  some  of  which  are  uncertain.  Finally,  sorption  as  well  as  precipitation  reac¬ 
tions  need  to  be  considered.  All  this  is  necessary  to  compute  the  fraction  of  a  chemical 
that  is  either  dissolved  or  particulate. 

The  equivalent  partitioning  model  employs  only  a  partition  coefficient  ttmu,  the 
ratio  of  particulate  to  dissolved  chemical  concentration. 


'^Mn  —  TT  /  tt\ - 


(3.2) 


where  Mn  {II)^  is  the  concentration  of  particulate  manganese  per  unit  sediment  vol¬ 
ume  and  m  is  the  concentration  of  solids  per  unit  sediment  volume.  Thus  Mn  {II)  /m 
is  the  concentration  of  manganese  per  unit  sediment  solids  (e.g.  Mn/g  solids).  Since 
Mn  {II)  is  the  concentration  of  dissolved  manganese  per  unit  volume  of  pore  water 
(e.g.  ixg  Mn/L),  ■kmu  has  the  units  of  L/kg. 

It  is  possible  to  build  some  level  of  chemical  realism  into  the  partitioning  formu¬ 
lation.  If  necessary,  can  be  varied  as  a  function  of  other  physical  and  chemical 
parameters  such  as  pH  in  order  to  produce  more  realistic  behavior.  The  practical  ques¬ 
tion  is:  does  the  added  difficulty  of  including  equilibrium  chemistry  into  the  model 
structure  result  in  added  realism?  Whatever  the  answer,  it  is  prudent  to  begin  the 
modehng  using  hnear  partitioning  and  examine  the  utility  of  the  results. 


3.1.2  Model  Formulation 

The  model  structure  is  shown  in  Fig.3.2.  It  is  an  elaboration  of  the  conceptual  model 
presented  in  Fig.3.1.  It  is  formulated  in  terms  of  manganese  but  it  applies  equally 
to  iron  since  the  redox  chemistry  of  these  two  metals  are  similar.  There  are  four 
dependent  variables:  Mn{II)  and  Mn02{s)  in  layers  1  and  2.  These  concentrations 
are  denoted  by  Mn(l)  and  Mn{2),  and  Mn02{l)  and  Mn02{2)  respectively.  They 
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Fig.  3.2  Schematic  of  the  manganese  flux  model 


correspond  to  the  total  Mn(^II)  and  Mn^IV^  in  each  layer.  The  source  of  manganese 
to  the  sediment  is  the  settling  of  particulate  manganese  oxide  from  the  overlying  water 
to  the  sediment.  The  flux  is  denoted  by  JMn02- 

Two  reactions  occur  in  the  aerobic  layer.  Mn(II)  partitions  to  form  particulate 
manganese.  This  reaction  is  parameterized  with  a  linear  partition  coefficient,  ttmti,!- 
In  addition,  dissolved  Mn(II)  is  oxidized  to  Mn02{s)  following  eq.(2.8)  with  first 
order  rate  constant  Thus  the  rate  of  Mn02  (s)  production  in  layer  1  per  unit 

area  is: 

SMn02  =  kMn,\HifdlMn{l)  (3.3) 

where  Hi  is  the  depth  of  the  aerobic  layer  and  fdi  is  the  fraction  of  Mn  (1)  that  is 
dissolved.  The  depth  of  the  aerobic  layer  Hi  is  estimated  from  the  ratio  of  the  diffusion 
coefficient  Di  in  layer  1  and  the  surface  mass  transfer  coefficient  s 
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where  the  surface  mass  transfer  coefficient  is  obtained  from  the  sediment  oxygen  de¬ 
mand  SOD  computed  by  the  oxygen  portion  of  the  sediment  flux  model,  and  the 
overlying  water  oxygen  concentration  O2  (0) . 


_  SOD 
"  ~  02(0) 

Thus  the  oxidation  rate  becomes 


(3.5) 


2 

=  kMn,l—  =  ^  (3.6) 

s  s 

where 

KMn,i  =  y/kMn,iDi  (3.7) 

which  is  the  definition  of  a  reaction  velocity  in  layer  1.  The  rational  for  this  definition 
and  eq.(3.4)  has  been  described  [Di  Toro  and  Fitzpatrick,  1993]. 

Two  reactions  also  occur  in  the  anaerobic  layer.  Mn(^II)  partitions  to  form  particu¬ 
late  manganese  which  is  parameterized  with  a  linear  partition  coefficient,  TrMn,2-  This 
may  be  different  from  the  aerobic  layer  partition  coefiicient  7iMn,i  due  to  differences  in 
pH,  the  nature  of  the  particles,  the  absence  of  freshly  precipitated  iron  oxy hydroxides, 
and  whatever  else  is  influencing  the  extent  of  the  formation  of  Mn  {II)  particulates. 
In  addition,  Mn02{s)  is  reduced  to  Mn{II)  following  eq.(2.4)  with  first  order  rate, 
kMn,2-  Hence  the  source  of  Mn{II)  in  the  anaerobic  layer  is: 


^Mn{II)  =  kMn,2H2Mn02{2)  =  KMn,2Mn02{2)  (3.8) 

where  H2  is  the  depth  of  the  anaerobic  layer  and 


I^Mn,2  —  kMn,2H2  (3-9) 

which  is  the  definition  of  a  reaction  velocity  in  layer  2  [Di  Toro  and  Fitzpatrick,  1993). 

The  mass  transport  between  the  overlying  water  and  layer  1  is  quantified  via  the 
surface  mass  transfer  coefiicient,  Klqi,  which  is  set  equal  to  s  =  50^/02 (0)  (eq.3.5) 

Klqi  =  s  (3.10) 

as  in  the  previous  models  for  oxygen  and  nutrient  fiuxes[Di  Toro  and  Fitzpatrick, 
1993].  Particle  mixing  with  mixing  velocity  w\2  and  diffusive  exchange  with  mass 
transfer  coefficient  Kli2  between  layers  1  and  2  are  included  as  before,  as  the  loss  of 
manganese  by  burial  with  sedimentation  velocity  W2. 


3.1.3  Equations  and  Solutions 

The  mass  balance  equations  for  the  model  follow  from  the  reactions  and  transport 
processes  discussed  above.  They  are: 
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Layer  1  Mn(II): 


0  =  Hi 


dMnjl) 


dt 


—  —  Mn  (0))  +  KLi2{fd2Mn{2)  —  fd\Mn{l)) 


+'Wi2{fp2Mn{2)  -  /piMn(l))  -  wsMn  (1) - 

Layer  2  Mn{II): 

0  =  =  -KLi2{fd2Mn{2)  -  /diMn(l)) 

-  wi2{fp2Mn{2)  -  fpiMn{l))  -  W2  (Mn(2)  -  Mn(l))  +  KMn,2Mn02{2) 
Layer  1  Mn02  (s) : 


(3.11a) 


(3.11b) 


0  =  ^  1j^f^^Mn{l)+wi2{Mn02i2)  -  Mn02{l))  +  JunO,  (3.11c) 

Layer  2  Mn02  (s): 


^dMn02{2) 


-KMn,2Mn02{2) 


—  wi2{Mn02{2)  —  Mn02{l))  —  W2  {Mn02{2)  —  Mn02{l)) 


(3.11d) 


The  particulate  fp  and  dissolved  fractions  are  computed  from  the  partition  coeffi¬ 
cients  and  the  concentration  of  sediment  solids,  m,  in  each  layer 


fdl  = 

1 

=  l-/pl 

(3.12a) 

1  -1- 

fd2  = 

1 

II 

1 

to 

(3.12b) 

1  -1-  TmtMn,2 

The  solutions  to  these  steady  state  equations  (3.11a-3.11d)  are  found  using  the 
symbolic  computation  program  MACSYMA  [Macsyma,  1993].  The  solutions  are  ex¬ 
tremely  complicated  if  the  dissolved  concentration  in  the  overlying  water  Mn  (0)  is 
included.  However  for  Mn  (0)  =  0  the  solutions  are  reasonably  concise.  The  reason 
is  that  for  this  case  only  one  equation  (eq.3.11c)  is  inhomogeneous^  with  an  external 
source  of  manganese.  The  rest  are  homogeneous  equations  with  terms  that  are  re¬ 
stricted  to  linear  functions  of  the  four  unknowns.  It  turns  out  that  the  equations  in 
such  a  set  have  relatively  concise  solutions. 


'An  equation  is  homogeneous  if  it  has  no  forcing  function.  An  inhomogeneous  equation  has  a  forcing 
function,  in  this  case  a  constant. 
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The  solutions  of  eqs. (3.1  la-3. lid)  for  the  case  Mn(0)  =  0  are; 

_  JMn02^Mn,2 


Mn(2)  = 


il^Mn,2  +  W2)  isfdiri2  +  W2)  +  fdl-^‘fl2‘W2 

S 


(3.13) 


where  ri2  is  the  ratio  of  the  layer  2  to  layer  1  concentrations.  It  is  found  to  be: 

_ fd2KLl2  +  fp2W\2 _ 


^  Mn(l) 
ri2  = 


f  (k  4. 

Jd\  (  Kli2  +  5-1 - 


(3.14) 


+  fpiwn  +  W2 


The  layer  1  solution  is  obtained  using  this  relationship: 

Mn(l)  =  ri2Mn{2) 


(3.15) 


The  manganese  flux  follows  from  the  relationship  between  surface  mass  transfer  and 
dissolved  layer  1  concentrations: 


J[Mn]  =  sfd\Mn{l) 


(3.16) 


Substituting  eqs.  (3. 13  and  3.15)  into  this  equation  yields: 

AMn\  =  - - -  p 

{K’Mn,2+W2)  (5/^1^2+102)  +  fdl—ri2W2 

s 

which  is  the  flux  of  manganese  to  the  overlying  water.  Various  special  cases  of  this 
equation  are  examined  below. 

Once  the  Mn  {II)  solutions  are  known  the  Mn02  {s)  concentrations  are  easily 
found.  They  are 


Mn02  (1)  = 


K*  \ 

fdi—^Mn  (1)  +  JmvO^  ]  {l^Mn,2  +  Wi2  +  W2) 


J 


(1012  +  t02)  {KMn,2  +  W2) 


(3.18) 


and 


Mn02  (2)  =  {fd2KL12  +  fp2Wi2  +  t02)  —  Mn  (1)  {fdlKii2  +  /plWl2  +  W2) 

^  Mn  j2 

(3.19) 

These  solutions,  eqs. (3. 13-3. 19)  can  be  checked  by  evaluating  the  global  mass  bal¬ 
ance  equation 

JMn02  =  J  [Mn]  +  W2  (Mn  (2)  +  Mn02  (2))  (3.20) 


which  states  that  the  input  of  particulate  manganese  to  the  sediment  JMn02  either 
returned  to  the  overlying  water  as  a  flux  of  dissolved  manganese,  or  is  buried  as 
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/*  List  the  solutions 


jmn=jmneqz; 


mn [1] =mnleq2; 


fdj  jmno2  kmno2  rl2  s 

knino2  rl2  +  fdj  rl2^  s  +  fdj  kmn rl2  +  knmo2  + 


jmno2  kmno2rl2 


^fdj  kmno2 rl2 +fdj  W2rl2^s  +  fdj  Icmn rl 2  +  kmno2  + 


mn [2] =mn2eq2 :mnleqz/rl2; 


iiDno2  kmno2 


^fdj  kmno2  rl 2  +  fdj  rl 2\  s  +  fd ^  kmn  rl 2  +  kmno2  +  w 


mno2 [2] =mno22z; 


01002^  = 


rl2=rl22; 


nin^(fd2kn3  +  fp2Wi^-.-w^)-mn|  (»,  +  W2) 


fdj  s  +fdj  ^12  *^^2  ''^12  '**''^2 


/*  Evaluate  the  mass  balance 

mb: jmno2-jmneqz-w[2] * (mn2eq2+mno22z) ; 

^  imno2  kmiio2 _ 

_  w  ^'2  s  +  fdj  kmn  rl2  +  knmo2  + 

(f<‘2  Ki2  ^  »2  ^12  ^  W2)  -  nm,  ^fd,  k„2  +  (p,  w^2  ■<■  W2) 
\  kmno2 

fdj  jmno2  kinno2  rl  2  s 

""  - - — r + jmno2 

^fdj  kmno2rl2+fdj  W2rl2^s+fdj  W2kiimrl2  +  W2kmno2 +W2 

/*  Substitute  the  solutions  for  Mn(l)  and  Mn(2) 

ratsimp(ev(mb,mn[l]=mnleqz,mn[2]=inn2eq2)  )  ; 

^  fdj  W2jmno2  rl2s  > 

+  [fdj  W2imno2  k^^^  +  fdj  jmno2  kmn  +  jmno2 j 

^  ♦  rl 2  -  fd^  jmno2  kj j2  “  ^P2  ^^2  ''’l2 ^ 

^fdj  kmno2  +  fd^  rl 2  s  +  fd ^  kmn  rl  2  +  kmno2  + 

/*  Substitute  for  rl2 


ratsimp (ev(%, rl2=rl2z) )  ; 


Fig.  3.3  Algebraic  evaluation  of  the  global  mass  balance  equation.  Listing  of  the  MACSYMA 
output  [Macsyma,  1993]. 
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either  Mn(II)  or  MnO^  (s) .  The  calculation,  which  is  performed  using  MACSYMA 
[Macsyma,  1993],  is  presented  in  Fig.3.3.The  solutions  are  Hsted,  the  mass  balance 
equation  is  formed,  and  the  results  are  simplified.  The  result,  in  the  last  line  of  the 
output,  is  zero  which  verifies  that  indeed  eq.(3.20)  is  satisfied. 


3.1.4  Simplified  solutions 


The  steady  state  solution  apportions  the  flux  of  manganese  dioxide  to  the  sediment 
JAfn02  fo  either  a  flux  of  Mn{II)  to  the  overlying  water  J  [Mn]  or  to  the  burial  of 
Mn02  (s)  and  Mn(II).  The  equations  can  be  simplified  considerably  by  examining  a 
special  case  which,  as  shown  subsequently,  appears  to  be  a  reasonable  representation 
of  what  actually  occurs. 

If  the  rate  of  reduction  of  Mn02  (s)  in  the  anaerobic  layer  is  rapid,  then  /CjUn02 
and  eq.(3.17)  becomes 

=  ^  (3.21) 

where 

Thus  p  controls  the  fraction  of  incoming  manganese  that  is  recycled  to  the  overlying 
water.  It  depends  on  the  ratio  of  the  burial  velocity  W2  and  the  surface  mass  transfer 
coefficient  s  modified  by  the  dissolved  fraction  in  layer  1  f^i,  and  ri2.  The  expression 
for  ri2  can  be  simplified  slightly  using  the  fact  that  the  most  of  the  manganese  (II)  in 
the  sediment  is  in  particulate  phases.  Thus 


fpi  ^  1 
fp2  ^  1 

and  ri2  becomes 


(3.23) 


ri2 


fdl 


_ fdzKLU  +  Wl2 


+  Wi2  +  W2 


(3.24) 


If  the  dissolved  fractions:  fdi  and  fd2  are  small  enough  and  s  is  not  small — note  the 
term — ^then 


ri2  ~ 


wn 


Wi2  +  W2 


(3.25) 


And  since  the  particle  mixing  velocity  1012  is  usually  large  relative  to  the  sedimentation 
velocity,  W2,  then  ri2  ^  1  as  indicated  in  eq.(3.25).  Thus  the  fraction  of  incoming 
particulate  Mn02  (s)  that  is  recycled  to  the  overlying  water  is  controlled  by 


W2 


18  STEADY  STATE  PARTITIONING  MODEL 


A  large  (3  corresponds  to  a  small  fraction  recycled  (eq.3.21).  A  large  /3  occurs  if 
the  partitioning  in  the  aerobic  layer  is  large  so  that  fd\  is  small.  To  understand  the 
magnitudes  involved  consider  a  typical  sedimentation  velocity  of  W2  =  10~®  m/d 
(0.36  cm/yr.).  A  typical  surface  mass  transfer  coefficient  is  s  =  SOD/O2  (0)  = 
(1  g/m^-d)/(l0  g/m^)  =0.1  m/d.  Thus  for  ~  1  which  corresponds  to  recycling 
one  half  of  the  incoming  flux,  fdi  ~  10“^  and  since  (eq.3.12) 


Mn 


(3.27) 


and  m  ~  1  kg/L  a  partition  coefficient  of  -kmu  ~  10'*  L/kg  is  required  if  fraction 
recycled  is  one  half.  If  the  partition  coefficient  is  much  larger,  then  none  of  the 
Mn02  (s)  flux  to  the  sediment  is  recycled  as  a  dissolved  flux  to  the  overlying  wa¬ 
ter  and  J  [Mn]  — >  0.  On  the  other  hand,  if  -nMn  is  much  smaller,  then  the  sediment 
traps  very  little  Mn,  all  the  incoming  flux  is  recycled,  and  J  [Mn]  ~  JMn02- 


3. 1.4.1  Anoxic  overlying  water  The  largest  manganese  fluxes  to  the  overlying  water 
are  observed  if  the  overlying  water  concentration  of  dissolved  oxygen  approaches  zero. 
Thus  it  is  important  to  examine  the  solutions  for  this  limiting  case.  For  O2  (0)  — v  0, 
s  =  SODf02  (0)  — >  00  and  equations  (3.13—3.16)  approach 


J[Mn]  =  JMnO.  f— V 

s-00  \KMn,2+W2j  \ 


fd2KL12  +  fp2Wi2 
fd2KLl2  +  fp2‘Wi2  -I-  W2 


) 


(3.28) 


In  this  case,  recycling  will  be  complete  if  both  of  the  bracketed  expressions  approach 
one.  The  first  necessary  condition  is  that 


«A/n,2  >  W2  (3.29) 

or,  using  the  definition  of  the  reaction  velocity  (eq.3.8),  if 

Wo 

kMn,2  »  TT  (3-30) 

rl2 

The  term  VJ2/H2  is  the  reciprocal  of  the  residence  time  of  Mn02  (s)  in  the  anaero¬ 
bic  layer.  Therefore  the  condition  in  eq.(3.30)  is  that  the  reduction  of  Mn02{s)  to 
Mn  (II)  occurs  rapidly  relative  to  the  residence  time  of  Mn02  (s)  in  H2. 

The  second  condition  is  that 


fd2l^L\2  +  fp2W\2  >  W2  (3.31) 

which  requires  that  the  rate  of  transport  of  dissolved  (fd2KL\2)  and  particulate  (/p2Wi2) 
manganese  (II)  to  layer  one  is  large  relative  to  the  rate  at  which  it  is  buried  W2.  This 
condition  is  necessary  so  that  the  manganese  produced  in  layer  2  by  the  reduction  of 
Mn02  (s)  to  Mn  (11)  is  transported  to  layer  1  rather  than  buried.  Once  it  is  in  layerl 
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it  escapes  to  the  overlying  water  without  further  modification  since  the  depth  of  the 
aerobic  layer  H\\s  approax;hing  zero 

lim  Hi=  lim  —  =  0  (3.32) 

s— ►OO  S—^OO  S  ^ 

and  no  re-oxidation  of  Mn  (II)  to  Mn02  (s)  can  occur. 

3.1.4.2  Aerobic  overlying  water  The  other  end  of  the  spectrum  is  the  behavior  of 
J[Mn]  as  s  0,  corresponding  to  low  SOD  and  high  overlying  water  dissolved 
oxygen.  To  examine  this  case  J  [Mn]  (eqs.3.13-3.16)  can  be  expanded  in  a  Taylor 
series  in  s: 

J  [Mn]  =  JmuO,  + . . .  (3  33) 

The  manganese  flux  decreases  quadratically  in  s  with  a  slope  determined  by  the  two 
reaction  velocities  and  the  sedimentation  velocity. 

The  results  of  this  section  have  shown  that  the  model  has  reasonable  limiting  be¬ 
havior  at  the  extremes  of  the  parameter  values.  It’s  ability  to  reproduce  observations 
will  be  examined  in  the  next  section. 

3.2  MANGANESE  FLUX  DATA 

Two  data  sets  will  be  examined  using  this  model.  This  first  is  a  relatively  small  num¬ 
ber  of  manganese  and  nutrient  flux  measurements  made  at  three  Long  Island  Sound 
stations  [Aller,  1980a].  The  second  is  a  large  number  of  manganese  [Hunt  and  Kelly, 
1988]  and  nutrient  flux  measurements  made  at  the  MERL  mesocosms  [Nixon  et  al., 
1986].  The  analysis  technique  is  to  attempt  to  reproduce  the  observed  relationship 
between  the  ammonia  J[NH^  and  manganese  J[Mn]  fluxes.  As  shown  in  Fig.3.4 
there  is  a  proportional  relationship  between  these  two  fluxes.  The  question  is:  what 
is  the  causal  hnkage. 

3.2.1  Relationship  between  Manganese  and  Ammonia  Fluxes 

The  parameter  in  the  manganese  model  equations  that  can  ultimately  be  related  to 
the  ammonia  flux  is  the  surface  mass  transfer  coefficient:  s  =  SOD /O2(0)-  This 
parameter  controls  both  the  rate  of  mass  transfer  from  the  pore  water  to  the  overlying 
water  (eq.3.11a)  and  the  depth  of  the  aerobic  layer,  iyi(eq.3.6).  As  SOD  increases  or 
the  overlying  water  dissolved  oxygen  decreases,  one  would  expect  that  J[Mn]  would 
increase  since  both  the  rate  of  surface  mass  transfer  increases  and  the  depth  of  the 
aerobic  zone  decrease.  The  latter  effect  decreases  the  residence  time  in  the  aerobic 
layer  and  makes  the  oxidation  of  Mn(II)  to  Mn02(s)  less  rapid  (eq.3.6),  enabling 
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Fig.  3.4  Manganese  versus  ammonia  flux.  Data  from  Long  Island  Sound  (A)  and  the  MERL 
mesocosm  experiment  (B). 

more  Mn{II)  to  escape  to  the  overlying  water.  As  shown  in  section  3. 1.4.1  page  18, 
this  is  indeed  how  the  steady  state  solution  behaves. 

The  dependency  of  s  on  ammonia  flux  occurs  because  the  ammonia  flux  and  SOD 
are  related  [Di  Toro  and  Fitzpatrick,  1993].  For  the  simple  case  with  zero  overlying 
water  ammonia  concentration,  the  relationship  between  ammonia  diagenesis  Jjv — ^the 
rate  at  which  organic  matter  is  mineralized  and  ammonia  is  released  to  the  pore 
water — and  ammonia  flux  to  the  overlying  water  J  [^^4]  is 

J[NHi]  =  -  (3-34) 

where  s  is  the  surface  mass  transfer  coefficient  and  knh4,i  is  the  reaction  velocity  for 
nitrification  in  the  aerobic  layer.  The  SOD  is  related  to  ammonia  diagenesis  since  as 
particulate  organic  matter  decays  to  produce  ammonia  at  a  rate  Jy,  it  also  produces 
reduced  organic  carbon  at  a  rate  Jc-  This,  in  turn,  is  eventually  oxidized  and  results 
in  SOD. 

The  temporal  variation  of  SOD,  and  the  fraction  of  reduced  end-products  that 
are  buried,  are  the  determined  by  the  full  sediment  model  equations  (Di  Toro  and 
Fitzpatrick,  1993].  However  an  approximate  analysis  can  be  made  as  follows.  The 
ratio  between  carbon  and  nitrogen  produced  by  diagenesis  is  closely  approximated  by 
the  Redfield  ratio  ac,N,  the  stoichiometric  composition  of  reactive  organic  matter 

Jc  =  o.c,nJn  (3.35) 

If  all  the  organic  carbon  were  oxidized,  the  SOD  would  be 


SOD  =  do^cJc  =  ^o-zCO-c,nJn 


(3.36) 
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where  aoa.c  is  the  Redfield  ratio  of  O2  to  carbon.  It  is  the  oxygen  required  to  oxidize 
the  organic  carbon  to  C02-  Hence  the  SOD  can  be  computed  from  ammonia  diagenesis 
Jjv  from  eqs.  (3.35-3.36).  The  equation  which  gives  ammonia  diagenesis  as  a  function  of 
ammonia  flux  and  overlying  water  oxygen  concentration  is  found  by  solving  eqs.  (3.34- 
3.36)  simultaneously.  The  result  is 


where: 


di 


+  J[NH4]  + 


(3.37) 


^  +  27J[NH4\  [i<'NHA,i02  (0)]^ 

+32  J[iViy4]K7v-/r4,lO2(0).^4a^2,iV‘^[-^^4]^  +  27[/C7Viy4,lO2(0)]2 


(3.38) 


This  approximation  assumes  that  all  carbon  diagenesis  eventually  becomes  SOD,  i.e. 
that  losses  due  to  sulfide  fluxes  and  burial  are  negligible.  Prom  an  analysis  of  Chesa¬ 
peake  Bay  model  fluxes,  these  losses  amount  to  no  more  than  25%  [Di  Toro  and 
Fitzpatrick,  1993]. 

The  second  source  of  error  is  due  to  the  time  lags  that  occur  between  the  production 
of  oxygen  equivalents  by  carbon  diagenesis  and  their  eventual  oxidation.  These  are 
caused  by  the  formation  and  the  subsequent  oxidation  of  FeS(s).  Prom  an  analysis 
of  the  steady  state  version  of  the  sediment  model  it  is  known  that  the  time  lag  effect 
causes  an  error  of  approximately  a  factor  of  two  between  the  carbon  diagenesis  esti¬ 
mated  from  SOD  assuming  steady  state,  and  the  actual  flux  [Di  Toro  and  Fitzpatrick, 
1993].  Thus  although  these  errors  are  not  negligible,  the  approximations  can  be  used 
so  long  as  the  magnitude  of  the  error  involved  is  recognized. 

Hence,  for  a  specified  ammonia  flux  J  [NH^l ,  ammonia  diagenesis  Jn  is  found 
(eq.3.37).  The  SOD  is  computed  using  eq.(3.36)  and  eq.(3.5)  yields  s  which  is  used  in 
eqs.  (3. 14-3. 17)  to  compute  J  [Mn] .  This  establishes  the  relationship  between  ammonia 
and  manganese  fluxes. 


3.2.2  Comparison  to  Data 

The  model  computations  are  compared  to  the  Long  Island  Sound  data  in  Fig.3.5.The 
model  parameters  are  listed  in  Table  3.1.  Two  cases  are  considered.  The  first  results 
from  assigning  a  slow  oxidation  rate  (small  KMn,i)-  The  input  flux  of  particulate 
manganese  JMn02  specified  to  be  150  mg  Mn/m^— d  for  the  Long  Island  Sound  case 
since  the  flux  of  manganese  from  the  sediment  J  [Afn]  is  approximately  this  magnitude 
at  the  high  ammonia  fluxes  (Fig.3.5A).  The  pore  water  and  sediment  data  are  plotted 
arbitrarily  at  an  ammonia  flux  of  10  mg  N/m^-d  for  convenience.  In  order  to  reproduce 
the  pore  water  and  sediment  manganese  concentrations  (Fig.3.5C,D)  it  is  necessary  to 
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Table  3.1  Parameter  Values  for  linear  partitioning  model’s  application  to  Long  Island  Sound 
(LIS)  and  MERL 


Parameters 

Values 

Units 

Small 

Large  K^n,! 

mg  Mn/m^-d 

LIS 

MERL 

LIS  MERL 

150 

20 

50  10 

0.25 

2.0 

m/d 

«Afn02.2 

0.1 

0.1 

m/d 

'RMn,l 

200 

200 

L/kg 

'RMn,‘2 

100 

100 

L/kg 

mi 

1.0 

1.0 

kg/L 

m2 

1.0 

1.0 

kg/L 

Kin 

0.05 

0.05 

m/d 

W12 

0.0012 

0.0012 

m/d 

W2 

0.5 

0.5 

cm/yr 

0.15 

0.15 

m/d 

[02(0)1 

5.0 

5.0 

mg/L 

ao2,c  X  o-c,N 

2.54  X  5.68 

2.54  X  5.68 

mg  02/mg  N 

recycle  most  of  the  incoming  manganese  to  the  overlying  water  (Fig.3.5B).  Hence  the 
oxidation  reaction  velocity  KMn.i  is  set  so  that  only  a  small  fraction  of  the  manganese 
(II)  flux  to  the  aerobic  layer  is  oxidized  to  Mn02  (s) .  As  a  result  very  little  is  trapped 
and  the  burial  fraction  is  small  (Fig.3.5B).  Since  most  of  the  flux  is  being  recycled, 
the  variation  in  s  which  occurs  as  J  [NH4]  varies  has  little  effect  on  the  manganese 
flux  as  shown  in  Fig.3.5A. 

An  alternate  calibration  is  shown  in  Fig.3.6.  The  oxidation  reaction  velocity  /CAfn.i 
is  increased  from  0.25  to  2.0  m/d  so  that  a  larger  fraction  of  the  manganese  (II)  flux  to 
the  aerobic  layer  is  oxidized  to  MnO^  (s) .  As  a  consequence,  the  sediment  manganese 
concentration  is  larger  Fig.3.6C  and  so  also  is  the  pore  water  concentration  Fig.3.6D 
because  of  the  partitioning  relationship  (eq.  3.2).  The  oxidation  rate  is  now  in  the 
range  where  the  variation  in  s  is  affecting  the  fraction  that  is  buried  versus  that 
recycled  Fig.3.6  B.  The  burial  fraction  is  predominant  at  low  ammonia  flux  and  the 
flux  to  the  overlying  water  predominates  at  high  ammonia  flux. 

A  similar  analysis  for  the  MERL  flux  data,  shown  previously  in  Fig.  3.4,  is  presented 
in  Fig.3.7  with  similar  results.  The  parameters  are  listed  in  Table  3.1.  The  flux 
of  particulate  manganese  to  the  sediments  appears  to  be  smaller  JMn02  ~  10“20 
mg  Mn/m^-d  than  the  case  of  the  Long  Island  Sound  data  JMn02  —  50-150  mg 
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Fig.  3.5  Manganese  linear  pj^titioning  model:  slow  cmdation  rate.  Comparison  of  model 
predictio^  (line)  to  data  (points)  from  Long  Island  Sound  as  a  function  of  the  ammonia 
flux.  (A)  Manganese  flux;  (B)  Mn  flux  normalized  to  depositional  flux,  (C)  pore  water  Mn 
concentration,  (D)  solid  phase  Mn  concentration. 


Fig.  3.6  Manganese  linear  partitioning  model:  rapid  oxidation  rate.  Comparison  of  model 
predictions  (line)  to  data  (points)  from  Long  Island  Sound  as  a  function  of  the  ammonia 
flux.  (A)  Manganese  flux;  (B)  Mn  flux  normalized  to  depositional  flux,  (C)  pore  water  Mn 
concentration,  (D)  solid  phase  Mn  concentration. 
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Fig.  3. 7  Manganese  linear  partitioning  model.  Ck^mparison  of  model  predictions  (line)  to 
data  (points)  from  the  MERL  experiment  as  a  function  of  the  ammonia  flux.  Top  group: 
slow  oxidation  rate.  Bottom  group:  rapid  oxidation  rate.  (A)  Manganese  flux;  (B)  Mn 
flux  normalized  to  depositional  flux,  (C)  pore  water  Mn  concentration,  (D)  sohd  phase  Mn 
concentration. 
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Mn/m^-d.  As  a  consequence,  the  pore  water  and  sediment  manganese  concentrations 
are  also  smaller. 

3.3  CONCLUSIONS 

The  steady  state  linear  partitioning  model  is  reasonably  successful  in  reproducing  a 
major  feature  of  the  Long  Island  Soimd  and  MERL  data  sets,  namely  the  correlation 
between  increasing  ammonia  and  manganese  fluxes.  However,  it  cannot  do  so  and  be 
consistent  with  the  pore  water  and  sediment  manganese  concentrations  at  the  same 
time.  Perhaps  the  problem  is  the  constant  linear  partition  coefficient  (eq.  3.2).  For 
a  steady  state  condition,  there  are  only  two  possible  exit  pathways  for  the  incoming 
manganese:  either  it  escapes  as  a  flux,  or  it  is  buried.  Therefore,  in  order  for  the 
model  to  reproduce  the  observations  it  is  necessary  that  it  predicts  a  higher  degree 
of  burial  at  low  ammonia  fluxes  and  a  larger  manganese  flux  to  the  overl3dng  water 
at  high  ammonia  fluxes.  If  the  actual  partitioning  were  such  that  more  manganese 
were  sequestered  at  high  ammonia  fluxes,  then  perhaps  the  solid  phstse  and  pore  water 
observations  could  be  satisfied  as  well  as  the  fluxes. 

Consider  the  following  possibility.  Suppose  the  fraction  of  manganese  that  is  par¬ 
ticulate  in  the  anaerobic  layer  is  determined  not  only  by  the  partitioning  of  Mn  {II) 
to  sediment  solids  but  also  by  the  solubility  of  MnC03{$),  This  is  not  a  linear  parti¬ 
tioning  process.  Rather  it  is  controlled  by  a  chemical  equilibrium  between  manganese 
and  carbonate  ion  COl”  concentrations.  Suppose  further  that  the  concentra¬ 
tion  of  manganese  carbonate  increases  with  increasing  ammonia  flux.  This  mechanism 
might  produce  the  necessary  variation  in  partitioning. 

The  reason  that  the  MnCOs{s)  concentration  might  increase  is  as  a  consequence 
of  increased  alkalinity  production  associated  with  increased  sulfate  reduction.  Sulfate 
reduction  increases  because  more  organic  matter  is  being  mineralized — ^this  must  be 
happening  for  the  ammonia  flux  to  increase — and  the  reduced  organic  carbon  is  being 
oxidized  using  sulfate  as  the  electron  acceptor. 

CH2O  +  ^SOf-  ^  ^C02  +  +  ^HS-  +  ^H20  (3.39) 

This  reaction  produces  alkalinity  HCO^  as  a  consequence  of  the  destruction  of  a 
strong  acid  cation  50|~ .  The  increased  alkalinity  causes  an  increased  precipitation  of 
MnCOz  (s) . 

In  the  next  chapter  we  shall  explore  the  behavior  of  such  a  model  using  a  much 
simplified  yet  reasonably  realistic  model  which  includes  chemical  equilibrium  control 
of  the  solid  phase. 
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Solubility  Control  of  Particulate  Manganese 

4.1  INTRODUCTION 

The  model  constructed  in  the  previous  chapter  used  linear  partitioning  to  determine 
the  fraction  of  manganese  in  the  anaerobic  layer  that  formed  a  solid  phase.  This  is 
a  critical  mechanism  since  it  influences  the  extent  of  burial.  And  the  more  burial 
that  occurs  the  less  is  the  flux  to  the  overlying  water.  In  this  chapter  we  examine 
the  behavior  of  a  similar  flux  model,  but  one  in  which  the  precipitation  of  manganese 
carbonate  is  explicitly  computed.  One  of  the  first  models  for  the  concentration  of 
manganese  in  the  pore  water  of  a  sediment  was  based  on  this  idea  [Holdren  et  al., 
1975;  Holdren,  1977). 

In  order  to  model  the  chemistry  of  manganese  carbonate,  we  need  to  be  able  to 
model  the  chemical  control  of  the  carbonate  concentration  in  pore  water.  Since  calcium 
carbonate  is  the  predominant  inorganic  csurbon  specie  in  sediments,  its  chemistry  is 
considered  flrst. 

4.2  CALCIUM  CARBONATE 

Most  freshwater  and  marine  sediments  contain  large  concentrations  of  calcium  car¬ 
bonate.  Typical  values  are  10-100  mg  CaCOzjg  or  1  to  10%  of  the  dry  weight.  This 
is  in  the  same  order  of  magnitude  as  the  organic  carbon  concentration  in  sediments. 
Calcium  cjirbonate  provides  a  buffer  system  for  the  pH  in  sediment  pore  water  and 
also  influences  the  concentration  of  carbonate  ion.  The  objective  of  the  model  formu¬ 
lated  in  this  section  is  to  reproduce  the  observations  of  pore  water  and  solid  phase 
concentrations  of  calcium,  alkalinity,  and  calcium  carbonate.  The  model  structure  is 
shown  in  Fig.4.1. 

The  only  chemical  mechanism  included  in  the  model  is  the  precipitation  of  calcium 
carbonate.  Chemical  speciation  in  the  pore  water  is  greatly  simplified.  Thus  the  model 
should  be  viewed  more  as  a  feasibility  investigation  into  the  methods  and  difficulties 
of  including  chemical  equilibrium  reactions  into  the  sediment  flux  framework  rather 
than  a  definitive  model  of  calcium  carbonate  formation  in  sediments. 
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Fig.  4.1  Schematic  diagram  of  the  calcium  and  alkahnity  flux  model. 


4.3  CHEMISTRY  AND  SIMPLIFICATIONS 

A  model  of  a  chemical  system  is  specified  by  the  chemical  components  and  the  species 
formed  by  the  components  [Morel,  1983].  This  type  of  system  is  termed  a  closed  system 
since  the  concentration  of  the  components  are  thought  of  as  being  constant  and  fixed. 
This  is  the  normal  situation  considered  in  typical  chemical  equilibrium  calculations. 
The  use  of  components  in  model  formulations  also  applies  to  systems  that  are  open 
[Di  Toro,  1976],  i.e.  systems  which  have  external  sources  and  sinks  of  the  components. 

The  equiUbrium  chemistry  of  carbon  dioxide,  alkalinity  and  calcium  in  natural  wa¬ 
ters  is  well  understood  [Stumm  and  Morgan,  1981;  Loewenthal  and  Marais,  1976; 
Morse  and  Mackenzie,  1990;  Butler,  1991].  It  is  specified  by  three  components:  cal¬ 
cium  Ca,  alkalinity  Alk,  and  total  inorganic  carbon  TIC.  However,  using  all  these 
three  components  results  in  an  equation  set  that  is  not  conveniently  solved.  Therefore 
a  simplified  set  of  components  and  equations  are  used  in  the  analysis  presented  below. 
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The  equation  that  determines  the  solubility  of  calcium  carbonate  is 

[Ca^^][COl-]  =  Ksi  (4.1) 

where  the  brackets  [  ]  denote  molar  concentrations  and  Kg\  is  the  solubility  product  of 
calcium  carbonate  CaCOz  (s).  Ionic  strength  corrections  are  assumed  to  be  included 
in  Ksi.  Since  the  objective  of  this  exercise  is  to  produce  a  tractable  set  of  model 
equations,  additional  approximations  are  introduced.  The  concentration  of  is 
approximated  with  the  concentration  of  total  dissolved  calcium  denoted  by  [Ca],  This 
ignores  the  contributions  of  complexes  of  calcimn  with  bicarbonate  CaHCO^  and 
other  ligands  to  the  total  dissolved  calcixim. 


[Ca]  =  [Ca^+]  +  [CaHCOt]  +  •  ■  •  [Ca^+j  (4.2) 

The  concentration  of  carbonate  can  be  obtained  using  the  definition  of  alkalinity 

[Aik]  =  [HCO^  ]  +  2[CO\~]  +  [OH~]  —  [iJ+j  (4.3) 


In  the  pH  range  that  is  typical  of  sediment  pore  waters  pH  =  7-8  and  for  large 
enough  alkalinity,  [All^  >  0.1  meq/L,  and  the  bicarbonate  ion  makes  up  essentially 
all  the  alkalinity 

[Aik]  ~  [HCO^]  (4.4) 

The  deprotonation  reaction  for  bicarbonate  HCO^  which  yields  carbonate  COf"  is 


Hco^  ^  fr+  +  co|~ 

and  the  mass  action  equation  is 


[H^][COl-] 

[HCO^] 


(4.5) 

(4.6) 


where  K2  is  the  equilibrium  constant  for  the  reaction  4.5.  Using  eq.(4.4)  for  HCO^ 
yields 


[H^WCOt] 

[Aik] 


=  K2 


(4.7) 


Thus  the  carbonate  concentration  becomes 


K2[Alk] 

[H+] 


(4.8) 


Substituting  this  equation  in  the  solubility  mass  action  equation  (4.1)  yields 


[Co\[Alk\  =  ^[H*\  A  KcMk 


(4.9) 


30  SOLUBILITY  CONTROL  OF  PARTICULATE  MANGANESE 


where  KcaAik  is  the  apparent  solubility  constant  of  CaCO^is)  in  terms  of  total  dis¬ 
solved  calcium  [Cd\  and  alkalinity  It  is  an  apparent  or  conditional  solubility 

constant  because  it  is  a  function  of  pH.  The  effect  of  decreasing  pH,  which  increases 
H'^,  increases  the  apparent  solubility  constant,  as  it  should  since  [CO^"]  is  decreasing 
as  pH  decreases  (eq.4.8). 

This  simplification  reduces  the  number  of  components  to  two:  Ca  and  Aik,  and 
requires  only  that  the  pH  be  specified.  This  is  a  significant  reduction  in  the  complexity 
of  the  equations. 

4.4  CLOSED  SYSTEM 

Consider,  initially,  a  closed  system  with  no  inflows  or  outflows.  Let  Car  and  Alkx  be 
the  total  concentrations  of  calcium  and  alkalinity  in  the  system.  The  mass  balance 
equations  are 

[Ca]  +  [CaC03(s)]  =  Cax  (4-10) 

and 

[Aik]  +  2[CaC03(s)]  =  Alkr  (4.11) 

where  [CaC03(s)]  is  the  concentration  of  calcium  carbonate  that  forms.  The  term 
2[CaC03(s)]  in  the  alkalinity  equation  arises  from  the  two  equivalents  of  alkalinity  in 
each  mole  of  CaCOz{s) — see  eq.(4.3).  Substituting  these  equations  into  the  mass  ac¬ 
tion  eq.(4.9),  dropping  the  brackets,  and  denoting  CaCOz  (s)  by  CaCOz  for  notational 
convenience  yields 


{Cax  —  CaCOz)  {Alkx  —  2CaCOz)  =  KcaAik 


(4.12) 


which  is  a  quadratic  equation  that  can  be  solved  for  CaCOz- 

It  is  more  convenient  and  instructive  to  solve,  instead,  for  dissolved  calcium  Ca^ 

Cad  =  Cax  —  CaCOz  (413) 

Finding  the  equation  for  Cad  yields 

Cal  -  (Car  -  \AlkT)Cad  -  =  0  (4.14) 

The  point  is  that  the  concentration  of  dissolved  calcium  is  a  function  only  of  the 
diflFerence  Car  —  ^ Alkr 

(4.15) 


Ca^  4  Car  - -Alkr 


which  we  define  as  the  excess  calcium  Ca^.  Hence,  the  equation  for  Cad  is 

Ca^-(Ca,)Carf-:^^  =  0 


(4.16) 
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The  reason  that  the  dissolved  calcium  concentration  Caj,  is  only  a  function  of  ex¬ 
cess  calcium  is  a  reflection  of  the  chemical  fact  that  changing  the  concentration  of 
solid  phase  calcium  carbonate  in  a  solution  that  is  already  saturated— i.e.  already  in 
equilibrium  with  CaCOz  (s) — does  not  affect  the  concentration  of  dissolved  calcium. 
The  system  is  already  saturated  with  respect  to  calcium  carbonate  so  that  solid  cal¬ 
cium  carbonate  is  present  and  adding  more  solid  calciiun  carbonate  does  not  cause 
any  chemical  change.  The  equation  for  calcium  carbonate  does  not  have  this  property 
since  the  concentration  of  calcium  carbonate  depends  on  both  the  concentration  of 
total  calcium  and  alkalinity 

To  see  that  excess  calcium  is  independent  of  CaCOz,  note  that  adding  a  quantity 
ACaCOz  of  calcium  carbonate  to  the  system  increases  Cax  by  ACaCOz  and  increases 
Alkr  by  2ACaCOz  since  each  mole  of  calcium  carbonate  is  two  moles  of  alkalinity 
However  the  excess  calcium,  Cox  (eq.4.15)  is  unchanged  and  therefore,  so  is  the  dis¬ 
solved  calcium  concentration  since  it  depends  only  on  the  excess  calcium  concentration 
(eq.  4.16). 

The  solution  for  dissolved  calcium  is  found  by  solving  eq.(4.16) 

Once  the  dissolved  calcium  is  known,  the  concentration  of  calcium  carbonate  follows 
from  eq.(4.10) 

CaCOz  =  Cax  —  Cad  (4- 18) 

We  shall  see  in  section  4.5  that  the  form  of  the  solutions  for  the  equations  that  result 
from  the  sediment  flux  model  are  similar. 

4.4.1  Solution  Behavior 

The  concentration  of  calcium  carbonate  is  determined  by  the  difference  between  the 
total  and  dissolved  calcium,  eq.(4.10)  which  in  turn  is  determined  by  the  excess  calcium 
(eq.4.17).  It  is  instructive  to  examine  the  solution  in  terms  of  the  two  independent 
variables;  total  calcium  and  alkalinity,  Co/r  and  Alkx- 

A  simple  way  of  understanding  the  solution  is  to  obtain  the  relationship  between  the 
concentrations  of  Car  and  Alkx  for  a  flxed  concentration  of  CaCOz-  The  solubility 
constraint  is 


[Ca][Alk\  =  {Car  -  CaCOz){AlkT  -  2CaCOz)  =  KcaAik  (4.19) 


For  a  fixed  concentration  of  CaCOz  this  equation  gives  a  relationship  between  Ca^ 
and  Alkx 


Co/j'  = 


KcaAlk 

Alkr  -  2CaCOz 


+  CaCOz 


(4.20) 
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Fig.  4.2  Concentration  of  CaCOz  (5)  as  a  function  of  total  alkalinity  and  total  calcium  con¬ 
centrations.  CaCOz  =  0  in  the  shaded  region.  The  dotted  lines  are  the  stoichiometric  limits 
to  the  quantity  of  CaCOz  (s)  that  can  form  as  Car  or  Alkr  — ►  00. 

For  various  concentrations  of  CaCOs^  the  relationship  between  Car  and  Alkx  is  shown 
in  Fig.4.2. 

The  lack  of  symmetry  in  Fig.4.2,  which  is  implied  by  eq.(4.9),  is  due  to  the  require¬ 
ment  that  two  moles  of  alkalinity  are  required  to  form  one  mole  of  calcium  carbonate 
whereas  only  one  mole  of  calcium  suffices.  If  the  graph  were  constructed  with  the 
axes:  Cax  and  2Alkx,  then  the  graph  would  be  symmetric  in  Cax  and  Alkx^ 

There  is  a  range  of  Cax  and  Alkx  over  which  no  precipitation  occurs.  This  is 
signaled  by  the  mathematics  if  either  of  the  resulting  dissolved  concentrations  are 
negative.  That  is,  if  either  dissolved  concentration,  computed  from  the  relationships 

Cad  =  Car  -  CaCO^  <  0  (4.21) 

or 

Alkd  =  Alkr  -  2CaCO^  <  0  (4.22) 

are  negative,  that  signals  that  the  solubility  constraint,  eq.(4.9)  cannot  be  satisfied 
with  chemically  reasonable,  i.e.  positive,  concentrations  and  the  only  realistic  solution 
is  CaCOz  —  0.  The  solubility  constraint  (eq.4.19)  becomes 

{Car  -  CaCOz)  {Alkr  -  2CaCOz)  =  {Car)  {Alkr)  <  KcaAik  (4.23) 

which  means  that  the  solution  is  undersaturated  with  respect  to  calcium  carbonate. 
The  shaded  region  in  Fig.4.2  is  where  no  CaCOz  can  form.  For  larger  concentrations, 
calcium  carbonate  forms  as  shown. 

This  example  illustrates  an  important  point.  When  computing  the  solutions  of  a 
chemical  equilibrium  problem  from  the  algebraic  equations  that  result  from  the  mass 
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balance  and  mass  action  equations,  the  solution  must  be  tested  to  insure  that  it  is 
physically  realistic,  namely  that  all  the  computed  concentrations — ^the  solid  phase 
and  dissolved  concentrations — are  positive.  The  mathematical  solutions  of  sets  of 
equations  where  the  inequalities  that  express  the  range  of  possibilities — saturated  and 
unsaturated  conditions,  e.g. 


(Car)  {Alkr)  <  KcaAik  (4.24) 

are  replaced  by  equalities,  e.g. 

{Car)  {Alkx)  =  KcaAik  (4.25) 

which  are  only  true  if  the  solid  phase  precipitates,  do  not  necessarily  result  in  positive 

concentrations.  We  shall  see  that  the  need  to  test  the  realism  of  the  solutions  extends 

to  the  more  complicated  settings  examined  below. 

4.5  SEDIMENT  MODEL  EQUATIONS  AND  SOLUTIONS 

The  model  formulation  is  shown  in  Fig.4.1.  The  dependent  variables  are  Alk{l)  and 
Alk(2),  the  total  alkalinity  in  layer  1  and  layer  2  respectively;  Ca{l)  and  Ca{2),  the 
total  calcium  in  the  same  sequence;  and  CaC03(2),  the  calcium  carbonate  in  layer  2. 
Calcium  carbonate  precipitation  in  layer  1  is  not  considered  since  the  layer  is  usually 
much  thinner  than  layer  2.  The  mass  balance  and  mass  action  equations  are  as  follows. 


Layer  1  Alkalinity 

0  =  s  [Alk(0)  -  Alk(l)]  +  Kli2  [Alk(2)d  -  Alk{l)]  (4.26) 

Layer  2  Alkalinity 

0  =  K[,i2  [Alk(2)d  —  Alk(l)]  —  W2  [2C7aC03(2)]  +  JAik  (4-27) 

Layer  1  Calcium 

0  =  s  [C7a(0)  -  <170(1)]  +  Kli2  [Ca{2)d  -  <7a(l)]  (4.28) 

Layer  2  Calcimn 

0  =  Kli2  [Ca{2)d  -  C7o(l)]  -  W2  [CaC03{2)]  (4.29) 

CaCOs  (s)  Solubility 

Ca(2)dAlk{2)d  =  KcaAik  (4.30) 

where 

Ca{2)^  =  Ca{2)-  CaCOz{2)  (4.31) 
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Aik  (2)^  =  Alk(2)  -  2  CaCOz{2)  (4.32) 

axe  the  dissolved  calcium  and  alkalinity  in  layer  2.  Since  no  calcium  carbonate  is 
allowed  to  precipitate  in  layer  1,  the  total  concentrations:  Ca{l)  and  Alk{l)  are  equal 
to  the  dissolved  concentrations. 

The  solution  to  these  equations  can  be  found  as  follows.  The  variables:  Alk(i'j, 
Alk{2),  C'a(l),  and  Ca(2)  are  eliminated  from  the  five  simultaneous  equations  (4.26- 
4.30).  The  resulting  equation  for  CaCOz{2)  is  a  quadratic  equation  of  the  form 

aCaCOs  "b  bCaCOs  -h  c  =  0  (4.33) 

where 

a  =  2(s  -I-  Kii2)‘^W2  (4-34) 

b  =  -W2{s  +  KLi2)[{‘2.Ca{Q)  +  Alk{Qi))sKLn  +  JAik{s  +  KLn)\  (4.35) 

^  =  -sKlx2  {sKl\2  [KcaAik  “  Alk(0)Ca{0)]  -  JAikCa(0)  {s  +  Km))  (4.36) 

It  is  remarkable  that  the  solution  to  these  five  simultaneous  equations  is  reduced  to 
the  solution  of  a  quadratic  equation,  as  it  is  for  the  closed  system  analyzed  in  the 
last  section  (eq.4.14).  This  suggests  that  so  long  as  the  chemistry  can  be  simplified  to 
the  point  that  the  equations  for  a  closed  system  reduce  to  a  manageable  result,  then 
the  open  system  equations  can  also  be  solved  at  the  same  level  of  algebraic  difficulty. 
We  shall  see  subsequently  that  this  observation  continues  to  be  true  for  more  complex 
systems. 

4.5.1  Simplified  Solution 

Before  examining  the  solution  of  eq.(4.33)  it  is  instructive  to  consider  the  limiting  case 
where  KcaAlk  0  which  corresponds  to  a  completely  insoluble  calcium  carbonate. 
In  fact,  this  is  a  somewhat  unrealistic  assumption  as  can  be  seen  in  Fig.4.2  where 
the  dotted  lines  represent  the  contours  for  the  KcaAik  0  case.  Nevertheless  the 
resulting  solution  is  instructive.  For  KcaAik  =  0,  eq.(4.30)  requires  that  either  Ca  (2)^ 
or  Aik  (2)^  be  zero  or,  using  the  definitions  (eqs.4.31-4.32) 

CaCOsi2)  =  \Alk{2)  (4.37) 

or 

CaCOz{2)  =  Ca{2)  (4.38) 

The  conditions  are  that  either  all  the  alkalinity  (eq.4.37)  or  all  the  calcium  (eq.4.38) 
becomes  calcium  carbonate,  whichever  is  the  least  abundant.  Assuming  that  alkalinity 
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is  in  the  shortest  supply  so  that  eq.(4.37)  is  true,  eqs.  (4.26-4.27)  become 

0  =  s  {Alk(0)  -  Alk{l)]  +  Kli2  [-Alk{l)]  (4.39) 

and 

0  =  Kli2  [-Alk{l)]  -  W2[2CaCOz{2)]  +  JAik  (4.40) 

Solving  eq.(4.39)  for  Alk{l)  3delds 

Aim  =  j:^^^Alk{Q)  (4.41) 

and  solving  eq.(4.40)  for  CaCOz{2)  yields 

CaCO,i2)  =  ^  +  (i  +  ^)  ■■  .4/fc(0)j  (4.42) 

To  understand  this  solution  consider  the  burial  flux  of  calcium  carbonate 

wzCaCOz  (2) 

JcaC03{,s)  =  W2CaCOz{2)  =  -  JAik  Alk(0)  (4-43) 

The  concentration  of  calcium  carbonate  is  determined  by  a  balance  between  the  loss 
by  burial  W2  CaCOz{2)  and  the  sources  of  alkalinity  to  the  sediment  on  the  right 
hand  side  of  eq.(4.43).  As  the  right  hand  side  of  eq.(4.43)  increases,  the  concentration 
CaCOz  (2)  increases  so  that  eq.(4.43)  remains  true. 

A  portion  of  the  alkalinity  transferred  to  the  sediment  is  via  diffusion  from  the 
overlying  water.  The  reciprocal  of  the  sum  of  the  reciprocals  of  the  mass  transfer 

coefficients  ^7  +  7c^)  total  mass  transfer  resistance  for  transferring  overl5dng 

water  alkalinity  Alk{0)  to  the  anaerobic  layer.  This  source  is  added  to  the  source  of 
alkalinity  due  to  sulfate  reduction  JAik-  Since  CaCOz  is  assumed  to  be  completely 
insoluble,  both  of  these  sources  of  alkalinity  are  transformed  entirely  into  calcium 
carbonate.  The  ^  reflects  the  two  moles  of  alkalinity  that  are  required  to  make  one 
mole  of  calcium  carbonate,  eq.(4.3). 

If  the  calcium  is  in  shortest  supply,  then  eq.(4.38)  is  true.  Using  this  result  and 
solving  eqs. (4.28-4.29)  for  JcaCOsis)  yields 

JcaCOs(s)=W2CaC03{2)=(^^  +  ^^  'c'a(O)  (4.44) 


where  the  calcium  carbonate  concentration  is  now  determined  by  a  balance  between 
the  loss  of  calcium  via  sedimentation  of  calcium  carbonate  W2  CaC03(s)  and  the 

source  of  calcium  from  the  overlying  water  ^5  +  Ca(0). 
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4.6  APPLICATION  TO  LONG  ISLAND  SOUND 

The  data  for  this  application  comes  from  observations  of  fluxes  from  three  stations  in 
Long  Island  Sound  ([Aller,  1980b]-[Aller  and  Yingst,  1980]).  The  parameters  required 
for  the  calcium-alkalinity  flux  model  are  the  usual  transport  parameters  as  well  as 
those  specific  to  these  components.  These  are  listed  in  Table  4.1. 

Table  4.1  Transport  parameter  values  for  the  calcium  carbonate  model 


Variable 

Value 

Units 

Ca(0) 

9 

mM 

Alk{0) 

2 

mM 

KcaAlk 

10-30 

mM^ 

In  addition  to  the  sotuce  of  alkalinity  from  the  overlying  water  there  is  also  a 
sediment  source  of  alkalinity,  JAik-  This  is  computed  from  the  ammonia  flux  using 
eqs. (3.5-3.38)  given  above  to  obtain  s  and  Redfield  stoichiometry  (eq.3.35)  is 
used  to  compute  the  carbon  diagenesis,  J^.  If  it  is  assumed  that  all  carbon  diagenesis 
reduces  sulfate  to  sulfide  then  the  following  formula  applies 

CH2O  +  ^SOf-  ^  HCO3  +  ^H2S  +  H2O  (4.45) 

Therefore  each  mole  of  organic  carbon  reacted  produces  one  equivalent  of  bicarbonate 
alkalinity.  However,  this  conclusion  depends  on  which  end-product  is  assumed  for  the 
sulfide  formed.  For  example  if  HS~  forms  then 

CH2O  +  ^SOl-  ^  ^C02  +  \hCO^  +  ^HS-  +  ^H20  (4.46) 

and  one  concludes  that  only  1/2  mole  of  alkalinity  forms.  If  iron  sulfide  is  the  final 
repository  of  the  sulfide,  which  is  assumed  in  the  sulfide  and  oxygen  flux  model  [Di  Toro 
and  Fitzpatrick,  1993],  then  the  reaction  is 

CH2O  +  ^FeOOH  +  ^50|-  IC02  +  IhCO^  +  ^FeS  +  \h20  (4.47) 

y  y  y  y  y  y 

which  is  almost  a  1  to  1  molar  ratio  between  CH2O  mineralized  and  alkalinity  HCO^ 
produced.  The  difference  between  eq.(4.45)  and  eq.(4.47)  is  the  1/9**'  mole  of  organic 
matter  required  to  reduce  the  Fe(III)  to  Fe(II)  which  does  not  produce  alkalinity. 
An  analysis  of  the  pore  water  relationships  in  Chesapeake  Bay  sediments  indicates 
that  a  1  to  1  stoichiometry  is  actually  observed  [Di  Toro  and  Fitzpatrick,  1993].  We 
adopt  this  stoichiometry  for  the  calculations  presented  below. 
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Fig.  4.3  Pore  water  concentration  of  (A)  dissolved  calcium  Ca,  (B)  alkalinity  Aik,  and  (C) 
the  apparent  solubility  product  Ca  x  Aik  from  three  Long  Island  Sound  stations. 

4.6.1  Chemical  Parameters 

The  only  chemical  parameter  required  KcaAik  is  the  apparent  solubility  of  CaCOz. 
Fig.4.3  presents  the  pore  water  data  for  dissolved  calcium  Ca  and  alkalinity  Aik  from 
three  Long  Island  Sound  stations.  The  product:  Ca  x  Aik  —  Kcaj^k  is  also  presented. 
It  ranges  from  approximately  20  (mM)^  to  over  40  (mM)^.  This  sets  the  range  for 
^CaAik-  With  the  parameters  established,  the  model  is  evaluated  by  specifying  an 
ammonia  flux  and  computing  the  resulting  concentrations.  The  additional  parameter 
values  are  presented  in  Table  4.2. 

Table  4.2  Parameter  values  for  the  calciiun  carbonate  model 


Variable 

Value 

Ca(0) 

9  mM 

Alk{0) 

2  mM 

KcaAlk 

10-30  mM2 

4.6.2  Results 

Fig.4.4  presents  the  results  for  KcaAik  =  10-30  (mM)^.  Alkalinity  flux  increases  as 
ammonia  flux  increases  (Fig.4.4A).  This  is  due  to  the  increased  production  of  alkalinity 
as  a  consequence  of  the  increase  in  sulfate  reduction  in  response  to  the  increase  in 
organic  matter  deposition  that  must  have  occurred  to  fuel  the  increase  in  ammonia  flux 
(eq.4.45).  For  KcaAik  =  30  (mM)^  the  solubility  of  CaCOz  (s)  is  exceeded  for  J[NH4] 
>  10  mg  N/m^-d  and  CaCOz{s)  starts  to  form  (Fig  4.4E).  At  that  point  the  calcium 
flux  changes  from  zero  to  negative  (Fig.4.4B)  which  is  the  flux  into  the  sediment  that 
is  required  to  support  the  burial  of  CaCOz  {s).  Also  the  flux  of  alkalinity  from  the 
sediment  to  the  overlying  water  increases  less  rapidly  (Fig.4.4A)  since  a  portion  is 
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Fig.  4.4  Comparison  of  model  results  to  Long  Island  Soimd  pore  water  data  for  various 
apparent  solubility  products  KcaAik- 


being  buried.  For  smaller  solubility  products,  KcaAik  =  10-20  (mM)^,  CaCOz{s) 
forms  over  the  entire  range  of  J[NH/^  investigated  (Fig.4.4E)  and  therefore  there  is  a 
corresponding  calciiun  flux  to  the  sediment  (Fig.4.4B). 

Pore  water  concentrations  are  compared  in  Fig.4.4C-D.  The  data  au’e  plotted  at 
J[NH^  =  10  mg  N/m^-d  for  convenience  only.  The  annual  average  ammonia  flux 
is  actually  somewhat  larger.  Alkalinity  and  calcium  are  reproduced  for  KcaAik  =  30 
(mM)^.  Calcium  concentrations  as  essentially  equal  to  the  overlying  water  concentra¬ 
tions  whereas  the  alkalinity  is  larger,  due  to  the  additional  source  that  results  from 
sulfate  reduction. 

The  concentration  of  CaCOz{s)  is  computed  to  increase  from  less  than  10  mg 
CaCOz  (s)/g  (1%  of  dry  weight)  to  over  100  mg  CaCOz  (s)/g  (10%).  Long  Island 
Sound  sediments  contain  between  25  and  200  mg/g.  The  model  reproduces  the  ob¬ 
servations  for  the  lower  values  of  the  solubility  constant.  This  is  in  contrast  to  the 
pore  water  results.  Perhaps  there  is  an  additional  source  of  calcium  carbonate  to 
the  sediments,  e.g.,  C'aC'03(s)  from  bivalve  shells,  which  accounts  for  the  additional 
CaCOz{s). 

Nevertheless,  the  model  is  reasonably  successful  in  reproducing  the  general  features 
of  the  data.  Pore  water  calcium  concentrations  are  predicted  to  be  close  to  the  over- 
lying  water  value,  9  mM  =  326  mg/L,  whereas  alkalinity  is  predicted  to  be  larger 
than  the  overlying  value  of  2  meq  =  100  mg/L.  Calcium  carbonate  concentrations  are 
predicted  to  be  in  the  range  of  10  to  100  mg/g  which  is  the  range  of  the  observations. 
These  results  are  obtained  using  transport  parameters  that  are  calibrated  from  Chesa¬ 
peake  Bay  sediment  flux  data,  suggesting  that  these  parameters  are  also  representative 
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for  these  Long  Island  Sound  sediments.  In  particular,  the  sedimentation  velocity  for 
these  sediments  is  quite  close  to  0.5  cm/yr.  used  for  these  calculations. 

4.7  MANGANESE  -  CALCIUM  -  ALKALINITY  FLUX  MODEL 

A  model  is  formulated  in  which  the  pore  water  and  solid  phase  concentration  of  man¬ 
ganese  are  controlled  by  the  solubility  of  manganese  carbonate.  To  this  end,  the 
calcium  carbonate  model  is  included  as  part  of  the  formulation.  The  rest  of  the  model 
parallels  the  linear  partitioning  model  presented  in  Chapter  3. 

4.8  CHEMISTRY  AND  SIMPLIFICATIONS 

The  simplification  of  the  manganese  carbonate  chemistry  parallels  that  for  calcium. 
The  concentration  of  Mn(//)2+  is  approximated  by  the  concentration  of  total  dissolved 
manganese,  Mn.  This  ignores  the  complexes  of  manganese  with  bicarbonate  and 
other  hgands.  The  concentration  of  carbonate  can  be  approximated  n.sing  alkalinity 
as  formulated  in  section  4.3,  eq.(4.8).  The  mass  action  equation  for  the  solubility  of 
MnCOz{s)  is: 


[Mn\[Mk\  =  A  (4.48) 

Pore  water  data  firom  Chesapeake  Bay  sediments  [Bricker  et  al.,  1977]  and  Long  Island 
Soimd  [Aller,  1980b]  are  shown  in  Fig.4.5.  MnAlk  apparent  solubility  products  are 
similar  for  both  sets,  ranging  from  0.1  to  1.0  mM^. 

4.8.1  Equations  and  Solutions  for  Manganese  and  Calcium  Carbonate 

The  mass  balance  and  mass  action  equations  for  a  closed  system  which  include  man¬ 
ganese  carbonate  as  well  as  calcium  carbonate  are: 

Alkalinity  mass  balance: 


Alkr  -  2CaAlk  -  2MnAlk  -  Alk^  =  0 

(4.49) 

Calcium  mass  balance: 

Cax  —  Cad  —  CaAlk  =  0 

(4.50) 

Manganese  mass  balance: 

Mut  ““  Mud  —  MnAlk  =  0 

(4.51) 

CaCOs{s)  Solubility: 

CadAlkd  =  KcaAlk 

(4.52) 
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Fig.  4.5  Top:  Pore  water  data  from  Chesapeake  Bay  sediments  [Bricker  et  al,  1977].  Bottom: 
Pore  water  data  from  Long  Island  Sound  [AUer,  1980a] 
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MnCOz{s)  Solubility: 

MudAlkd  =  KMnAik  (4.53) 

where  Alkx  and  Alkd  are  the  total  and  dissolved  alkalinity;  Mut  and  Mud  are  the 
total  and  dissolved  manganese,  and  Cclt  and  Ca,d  are  the  total  and  dissolved  calcium. 
Eliminating  all  the  independent  variables  except  MnAlk  yields  a  quadratic  equation 
of  the  form: 

aMnAlk^  +  bMnAlk  4-  c  =  0  (4.54) 

with  coefficients: 

a  =  —A{KMnAlk  +  KcaAlk)  (4.55) 

b  =  2{2KMnAlkMnT  +  AKcaAlkMriT  -  2CaTKMnAlk  +  AlkTKMnAlk)  (4.56) 

^  ~  ~‘^(.^^CaAlkAfn^  —  ^Ca^K MnAlk AIut  +  AlkxKMnAlkMnx  —  KudnAlk)  (4.57) 
The  equation  for  CaAlk  is  also  quadratic: 

aCaAlk"^  +  bCaAlk  +  c  =  0  (4.58) 

with: 

a  =  2  (4.59) 

b  =  2MnAlk  -  Cax  -  Alkx  (4.60) 

c  =  —2C(iTMnAlk  —  KcaAik  +  AlkxCax  (4.61) 

It  is  remarkable  that  adding  manganese  as  a  component  and  MnCOz{s)  as  a  solid 
phase  does  not  materially  complicate  the  solution.  It  requires  only  the  solutions  of 
quadratic  equations. 

The  solution  procedure  is  to  solve  the  quadratic  eq.(4.54)  for  MnAlk  which  yields 
two  solutions.  For  each  of  these  solutions,  two  solutions  for  CaAlk  are  found  from 
eq.(4.58).  For  each  of  the  four  possible  solutions  the  dissolved  concentrations  are 
computed.  Then  the  solutions  are  checked  for  the  following  conditions. 

1.  Are  all  concentrations  positive? 

2.  Are  there  no  oversaturated  solids? 

3.  If  both  of  these  conditions  are  true,  then  the  feasible  solution  has  been  found. 

4.  If  no  feasible  solution  is  foimd,  then  either  CaAlk  and/or  MnAlk  are  zero  and 
the  solution  is  undersaturated  with  either  or  both  of  these  solids. 
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4.8.2  Results 

An  example  computation  is  presented  below  in  Tables  4.3  and  4.4.  The  calcium  con¬ 
centration  is  kept  constant,  Ca  =  300  mg/L  =  7.5  mM  which  is  approximately  the 
concentration  in  overlying  water  for  the  Long  Island  Sound  sediment  data.  The  solu¬ 
bility  products  are:  KcaAlk  =  10  (mM)^  and  KmuAIL  =  0.4  (mM)^.  Total  alkalinity  is 
varied  from  1  to  20  mM  and  total  manganese  from  0.05  to  0.5  mM.  These  ranges  are 
representative  for  the  data  considered  above.  Note  that  the  alkalinity  is  much  larger 
than  the  manganese  concentrations.  The  tables  4.3  and  4.4  have  alkalinity  varying  in 
the  vertical  direction  and  manganese  in  the  horizontal  direction. 

The  computed  solid  phase  concentrations,  MnAlk  and  CaAlk,  and  the  computed 
solubility  products  are  given  in  Table  4.3.  CaAlk  forms  as  Alkr  is  increased.  Its 
concentration  is  not  affected  by  the  concentration  of  Mrix-  MnAlk  also  forms  as  Alk^ 
and  Mut  are  increased.  Its  concentration  is  affected  by  both  independent  variables: 
Alkx  and  Mnx-  Note  that  if  the  aqueous  phase  is  imdersaturated  with  either  solid — 
the  solubility  product  is  less  than  KcaAik  or  KMnAlk—tiien  no  solid  phase  forms.  The 
dissolved  concentrations  are  listed  in  Table  4.4.  They  are  all  positive  as  required  by 
the  feasibility  conditions. 

4.9  SEDIMENT  MODEL  EQUATIONS  AND  SOLUTIONS 

The  structure  of  the  manganese  flux  model  is  illustrated  in  Fig.4.6.  The  oxidation  of 
Mn{II)  to  Mn02{s)  occurs  in  the  aerobic  layer,  and  the  reduction  or  Mn02{s)  to 
Mn(II)  occurs  in  the  anaerobic  layer.  The  partitioning  of  Mn(II)  in  the  anaerobic 
layer  is  controlled  by  the  solubility  of  MnCOz{s).  In  order  to  calculate  the  carbon¬ 
ate  concentration,  the  influence  of  CaCOz{s)  must  also  be  considered.  Both  these 
reactions  are  shown  as  the  solubility  of  MnAlk  and  CaAlk,  as  discussed  above. 

Hence,  the  equations  for  the  sediment  model  are  a  combination  of  the  equations 
for  the  linear  partitioning  model  for  manganese,  with  the  addition  of  the  equations 
for  calcium,  from  the  calcium-alkalinity  model,  and  a  solubility  equation  for  MnAlk. 
The  equations  and  the  solution  procedures  are  listed  in  the  Appendix.  Once  again, 
the  solutions  are  found  from  the  roots  of  quadratic  equations.  They  are  very  similar 
to  the  closed  system  equations,  although  the  coefficients  of  the  equations  are  much 
more  involved.  Nevertheless,  no  essential  complexity  is  introduced  and  the  numerical 
solutions  are  straightforward. 

4.9.1  Results 

A  calibration  to  the  Long  Island  Sound  data  is  shown  in  Fig.4.7.  The  additional  param¬ 
eter  values  are  presented  in  Table  4.5.  The  observed  pore  water  alkalinity  (Fig.4.7C), 
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Table  4.3  Concentration  and  solubility  products  of  MnAlk  and  CaAlk  as  a  function  of  Alkr 
and  MtiT 


Concentrations 


Alhr  (mM) 

MnAlk  (mM) 

1.0000 

0 

0 

0 

0 

0 

0.699 

1.8206 

0 

0 

0 

0 

0.0309 

0.2094 

3.3145 

0 

0 

0 

0 

0.0777 

0.2566 

6.0342 

0 

0 

0 

0.0389 

0.1525 

0.3325 

10.9856 

0 

0.0045 

0.0504 

0.1230 

0.2382 

0.4207 

20.0000 

0.0198 

0.0490 

0.0952 

0.1685 

0.2847 

0.4688 

Alkx  (mM) 

CaAlk  (mM) 

1.0000 

0 

0 

0 

0 

0 

0 

1.8206 

0.4101 

0.4101 

0.4101 

0.4101 

0.3844 

0.2348 

3.3145 

1.6152 

1.6152 

1.6152 

1.6152 

1.5548 

1.4146 

6.0342 

3.5210 

3.5210 

3.5210 

3.4971 

3.4266 

3.3132 

10.9856 

5.6321 

5.6309 

5.6190 

5.5998 

5.5689 

5.5184 

20.0000 

6.7445 

6.7429 

6.7404 

6.7364 

6.7299 

6.7195 

Mut  (mM) 

0.0500 

0.0792 

0.1256 

0.1991 

0.3155 

0.5000 

Solubility  products  of  MnAlk  and  CaAlk 


Alkx  (mM) 

[Mn][Alk]  (mM)2 

1.0000 

0.0500 

0.0792 

0.1256 

0.1991 

0.3155 

0.4000 

1.8206 

0.0705 

0.1118 

0.1771 

0.2808 

0.4000 

0.4000 

3.3145 

0.0850 

0.1347 

0.2134 

0.3382 

0.4000 

0.4000 

6.0342 

0.1257 

0.1992 

0.3156 

0.4000 

0.4000 

0.4000 

10.9856 

0.2677 

0.4000 

0.4000 

0.4000 

0.4000 

0.4000 

20.0000 

0.4000 

0.4000 

0.4000 

0.4000 

0.4000 

0.4000 

Alkx  (mM) 

[Ca\[Alk]  (mM)^ 

1.0000 

7.5000 

7.5000 

7.5000 

7.5000 

7.5000 

6.9756 

1.8206 

10.000 

10.000 

10.000 

10.000 

10.000 

10.000 

3.3145 

10.000 

10.000 

10.000 

10.000 

10.000 

10.000 

6.0342 

10.000 

10.000 

10.000 

10.000 

10.000 

10.000 

10.9856 

10.000 

10.000 

10.000 

10.000 

10.000 

10.000 

20.0000 

10.000 

10.000 

10.000 

10.000 

10.000 

10.000 

Mnx  (mM) 

0.0500 

0.0792 

0.1256 

0.1991 

0.3155 

0.5000 
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Table  4.4  Dissolved  concentration  of  Mn,  Aik  and  Ca  as  a  function  of  Alkr  and  Mut 


Alkr  (mM) 

Mn  (mM) 

1.0000 

0.0500 

0.0792 

0.1256 

0.1991 

0.3155 

0.4301 

1.8206 

0.0500 

0.0792 

0.1256 

0.1991 

0.2846 

0.2906 

3.3145 

0.0500 

0.0792 

0.1256 

0.1991 

0.2378 

0.2434 

6.0342 

0.0500 

0.0792 

0.1256 

0.1601 

0.1629 

0.1675 

10.9856 

0.0500 

0.0748 

0.0752 

0.0760 

0.772 

0.0793 

20.0000 

0.0302 

0.303 

0.304 

0.0305 

0.0308 

0.0312 

Alkx  (mM) 

Aik  (mM) 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.9301 

1.8206 

1.4105 

1.4105 

1.4105 

1.4105 

1.4105 

1.3764 

3.3145 

1.6993 

1.6993 

1.6993 

1.6993 

1.6820 

1.6433 

6.0342 

2.5132 

2.5132 

2.5132 

2.4982 

2.4550 

2.3885 

10.9856 

5.3535 

5.3502 

5.3163 

5.2627 

5.1784 

5.046 

20.0000 

13.236 

13.208 

13.164 

13.095 

12.985 

12.812 

Mnx  (mM) 

0.0500 

0.0792 

0.1256 

0.1991 

0.3155 

0.5000 

Alkr  (mM) 

Ca  (mM) 

1.0000 

7.5000 

7.5000 

7.5000 

7.5000 

7.5000 

7.5000 

1.8206 

7.0899 

7.0899 

7.0899 

7.0899 

7.1156 

7.2652 

3.3145 

5.8848 

5.8848 

5.8848 

5.8848 

5.9452 

6.0854 

6.0342 

3.9790 

3.9790 

3.9790 

4,0029 

4.0734 

4.1868 

10.9856 

1.8679 

1.8691 

1.8810 

1.9002 

1.9311 

1.9816 

20.0000 

0.7555 

0.7571 

0.7596 

0.7636 

0.7701 

0.7805 

Mnx  (mM) 

0.0500 

0.0792 

0.1256 

0.1991 

0.3155 

0.5000 
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Fig.  4.6  Schematic  diagram  of  the  Mn  —  Ca—  Aik  Model 


calcium  (Fig.4.7D),  and  CaCOz{s)  (Fig.4.7E)  are  reproduced  by  the  model  calcula¬ 
tions.  The  data  are  plotted  at  JlNH^]  =  10  mg/m^-d  for  convenience.  The  fit  is 
actually  slightly  better  than  the  model  for  CaCOz{s)  alone,  Fig.4.4. 

The  manganese  results  are  shown  in  Fig.4.8.  Pore  water  manganese  is  reproduced 
reasonably  closely  (C)  while  the  computed  concentration  of  MnCOz{s)  is  higher  than 
the  observation  (D).  The  manganese  flux  varies  with  respect  to  ammonia  flux  (A,B). 
However  it  is  still  not  as  pronounced  as  the  observations  (E). 

4.9.2  Conclusions 

The  inclusion  of  manganese  carbonate  precipitation  as  the  speciation  reaction  that 
controls  manganese  partitioning  in  the  anaerobic  layer  has  slightly  improved  the  model 
results.  However  if  the  results  are  compared  to  the  linear  partitioning  case — compare 
Fig.  3.6  to  Fig.4.8 — ^it  is  clear  that  the  model  is  producing  substantially  the  same 
results.  The  addition  of  the  chemical  control  on  partitioning  in  the  anaerobic  layer  has 
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Calcium  Flux  Pore  Water  Calcium 


J(NH4  (mg  N/m^-d)  JINH4  (mg  N/m^-d) 


Fig.  4.7  Comparison  of  mociel  and  data  from  Long  Island  Sound.  Pore  water  calcium,  alka¬ 
linity,  and  sediment  calcium  carbonate  concentrations. 


Fig.  4.8  Comparison  of  model  and  Long  Island  Sound  data.  Pore  water  and  sediment  man¬ 
ganese  concentrations.  Manganese  versus  ammonia  flux. 
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Table  4.5  Parameter  values  for  the  MnCOz  and  CaCOz  model 


Variable 

Value 

Units 

0.03 

m/d 

KAfn,2 

10 

m/d 

pil 

300 

L/kg 

KcaCOs 

20 

mM2 

KcaCOa 

0.5 

mM2 

not  fundamentally  changed  the  model’s  behavior,  at  the  cost  of  an  order  of  magnitude 
increase  in  complexity  of  the  resulting  equations.  Thus,  this  approach  is  abandoned 
and  the  linear  partitioning  formulation  is  retained. 

The  next  step  is  to  examine  the  time  variable  version  of  the  model,  which  is  pre- 
sented  in  the  next  chapter. 
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Appendix:  Equations  and  Solutions  for  MnCOa,  CaCOa,  Alkalinity  and  TIC 

The  mass  balance  equations  are: 

Layer  1  Alkalinity 

Kli2  i-2MnAlk  -  2CaAlk  +  Aik  (2)  -  Aik  (1))  -  (Aik  (1)  -  Aik  (0))  s  =  0  (A.l) 
Layer  2  Alkalinity 

0  =  -W2  {2MnAlk  +  2CaAlk) 

Kij\2  (  2MnAlk  —  2  CaAlk  -I-  Aik  (2)  —  Aik  (1))  ■!“  Jail 

Layer  1  Calcium 

Cd  (1)  s  +  CcL  (0)  s  +  ( — CcbAlk  -|-  Ca  (2)  —  Ca  (1))  Kli2  ~  0  (•^•3) 

Layer  2  Calcium 

~  {—CaAlk  +  Ca  (2)  —  Ca  (1))  Kl\2  —  W2  CaAlk  =  0  (•^•4) 

Layer  1  Mn{II) 

0  =  -  (Mn  (1)  fd  -  Mn  (0))  s  +  wn  (MnAlk  -  Mn  (1)  /p) 

+Kli2  (-MnAlk  -  Mn  (1)  fa  +  Mn  (2))  -  ki  Mn  (1)  fa  (A.5) 

Layer  2  Mn(II) 

0  =  —wi2  (MnAlk  —  Mn  (1)  fp)  —  W2MnAlk 

—Kli2  (—MnAlk  —  Mn  (1)  fd  +  Mn  (2))  +  K2  MnOa  (2)  (A.6) 

Layer  1  Mn02(s) 

JMn02  +  ^<^1  Mn  (1)  fd  +  (Mn02  (2)  —  MnOa  (1))  wi2  =  0  (A.7) 

Layer  2  Mn02(s) 

—  (Mn02  (2)  —  Mn02  (1))  rui2  —  Mn02  (2)  W2  —  K2  Mn02  (2)  =  0  (A.8) 

The  mass  action  equations  are 
Layer  2  CaCOz(s)  Solubility 

(Ca  (2)  -  CaAlk)  (-2MnAlk  -  2  CaAlk  +  Aik  (2))  -  K^p^CaCOz  =  0  (A.9) 
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Layer  2  MnCOz{s)  Solubility 

(-2MnAlk  -  2  CaAlk  +  Aik  (2))  (Mn  (2)  -  MnAlk)  -  Ksp^MnCOa  =  0 


(A.10) 


The  first  step  in  the  solution  is  to  solve  for  all  the  dependent  variables  except  CaAlk  (2) 
and  MnAlk{2)  using  the  mass  balance  equations.  The  solutions  are 


Alk{l) 


AlkjO)  s-2KLi2MnAlk-2CaAlkKLi2+Alk(2)  Kli2 

S  +  Kli2 


Ca{l) 


Ca  (0)  s  -  CaAlk  Klu  +  Ca  (2)  Kli2 
S  +  Kli2 


_  Afn(O)  s  —  Ki,i2MnAlk+wi2  MnAlk+Mn(2)  Kli2 
fds  +  wi2  fp  +  fd Kiri  +  fd 


(A.11) 

(A.12) 

(A.13) 


Aik  {2)  = 


2  Ki,i2MnAlk  s  —  2  W2MnAlk  s  +  2  CaAlk  Klu  s 
+AZA:(0)  Kli2 s  —  2w2 CaAlk s  +  Jahc s 
-2  W2  Ki,i2MnAlk  —2  W2  CaAlk  Km  +  JMk  Km 
Kms 


Ca{2)== 


K2 


Mn  (2)  = 


Mn02  (1)  = 


CaAlk  Km  s  +  Ca  (0)  Km  s  —  W2  CaAlk  s-W2  CaAlk  Km 

Kms 

Mn  (0)  Wi2  fpS  +  fd  JMn02  s  -  Wi2  fd  MnAlk  s 
W2  fdMnAlk  s  +  wi2  JmtiO^  fp  -  '^2  wi2MnAlk  fp 
f  Mn  (0)  wi2  fpS  —  Wi2  fd  MnAlk  s  —  W2  fdMnAlk  s 
\  —W2  wi2MnAlk  fp 

+/C2  fd  Km  {MnAlk  s  +  Mn  (0)  s  +  JMn02  ~  'W2MnAlk) 
+W2  fd  Km  {MnAlk  s  +  Mn  (0)  s  —  W2MnAlk) 

+Ki  «2  fd  {Mn  (0)  s  +  JMn02  -  W2MnAlk) 

+K1W2  fdKmM''T'Alk—KiW2  fd  {wi2  +W2)MnAlk 


(A.14) 

(A.15) 


) 


+W2 


) 


W2  fd  Km  s  +  K2fd  Km  s  +  kiW2  fd  Km 

{■W12  +W2+  K2) 

{JMn02  S  +  Mn  (0)  Kl  S  +  Ki  JMn02  ~  W2MnAlk) 
W\2  {W2  s  +  K2  S  +  /Cl  W2) 


M  O  (2)  =  S+Mn{0)  Kis  +  ki  JMn02  -  W2MnAlk 

W2  S  +  K2  s  +  K1W2 


(A.16) 

(A.17) 

(A.18) 


The  mass  action  equations  are  used  to  solve  for  CaAlk  and  MnAlk.  Only  the  layer 
two  solutions  that  are  needed  in  these  equations.  They  are  written  in  the  following 
form  in  order  to  substitute  for  various  parameter  groups  and  to  isolate  the  dependency 
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on  CaAlk  and  MnAlk: 


Alk(2')  —  JaO  "I"  ,^acCciAlk(2^  +  JamAfnAlk(2,') 
Ca(2)  =  JcO  +  JccCaAlk{2)  +  JcmMnAlk{2) 

Mn{2)  =  Jmo  +  JmcCaAlk{2)  +  JmmMnAlk{2) 


(A.19) 

(A.20) 

(A.21) 


where  the  notation  denotes  the  Jacobian  {J)  of  the  equation  (a,  c,  m  for  Aik,  Ca,  Mn) 
with  respect  to  the  variables  c,m  for  CaAlk,  MnAlk,  or  0  for  the  constant  term. 


JaO  = 


A/fc(0)  Kli2  s  +  Jail  s  +  Jaiic  Km 

Kli2  s 


Jcc  — 


^  _  2  {Kl\2  S-W2S-W2  Kli2) 

c/OC  -  7-^ 

Ktns 

^  _  2  (Kli2  S—W2S—W2  Kli2) 

dam  — 

Kl12S 

JcQ  =  Ca{0) 

_  K112  S  -W2S  -W2  Kli2 
dec  ”” 

Kl12S 
Uem  —  0 

Mn(0)  W2  Wi2  fp  S+Mn(0)  K2  W12  fpS  +  K2  fd  JMn02  S 
+Mn(0)  W2  fd  Kli2  S+Mn(0)  K2  fd  Km  s+Mn{0)  ki  K2  fd  s 
_  +/C2  W12  JmuO^  fp  +  ^2  fd  Km  JMnOj  +  Ki  K2  fd  JmtiO^ 

~  fd  Km  {W2S  +  K2S  +  Ki  W2) 

ffmc  ~  0 

W2  fd  Km  S  +  K2fd  Km  S-W2  Wi2  fdS-K2  Wi2  fd  S 
-W2  fdS  -  K2W2fdS  -  wl  'W12  fp  -  K2  W2  W12  fp 

-wl  fd  Km  -  K2  W2  fd  Km  +  ki  W2  fd  Km  -  «i  W2  wi2  fd 

_ -KlW^  fd  -  Ki  K2W2fd _ 

”  fd  Km  {W2S  +  K2S  +  Ki  W2) 

An  important  relationship  among  these  coefficients  is; 

Sac  —  ^Scc 


SwTi  — 


(A.22) 


(A.23) 

(A.24) 

(A.25) 

(A.26) 

(A.27) 


(A.28) 

(A.29) 


(A.30) 


(A.31) 


O^am  — 


(A.32) 


Using  these  equations  and  substituting  into  the  mass  action  equations  yields  quadratic 
equations  for  both  CaAlk  and  MnAlk.  The  form  is: 

aCaAlk^  +  bCaAlk  +  c  =  0  (A.33) 
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and  similarly  for  MnAlk.  The  coefficients  for  MnAlk(2)  are  denoted  by  the  subscript 
Mn: 


^Mn  “^KcaAlk  +  ‘^^mm^CaAlk  +  Jmm{~^KcaAlk  ~ 


'\~^KMnAlk  >^cci  "^^MnAlk  "b  "^J^mmKji^fnAll^  (A.34) 

bjifn  =  <Jmo{—^KcaAlk  +  ^SmmKcaAlk  ~  “^K^nAlk 
MnAlk  MnAlk) 

"b <Jao(,  K MnAlk  "b  3mmK MnAlk  (A.35) 

CMn  =  —(—‘^JmO^CaAlk  +  KMnAlk 

+JTnQ{—{,JaQK MnAlk)  +  “^JdoKMnAlk))  (A.36) 

and  for  CaAlk(2)  which  are  denoted  by  the  subscript  Ca: 

aca  =  2-  4Jcc  +  2J^  (A.37) 

bca  =  Jaoi-l  +  Jcc)  +  (-2  +  2Jcc)J^  +  2MnAlk(2) 

-4JccMnAlk{2)  +  2J^MnAlk(2)  (A.38) 

cca  =  -{-(JooJdo)  +  KcaAik  +  Jdo{2MnAlk{2)  -  2JccMnAlk(2))  (A.39) 


The  solution  procedure  is  to  solve  the  quadratic  equation  for  MnAlk(2)  which  yields 
two  solutions.  For  each  of  these  solutions,  two  solutions  for  CaAlk{2)  are  found.  For 
each  of  the  four  possible  solutions  the  rest  of  the  variables  are  computed.  Then  the 
solutions  are  checked  for  the  following  conditions.  Are  all  concentrations  positive?  Are 
there  no  oversaturated  sohds?  If  both  of  these  conditions  are  true,  then  the  feasible 
solution  has  been  found.  If  no  feasible  solution  is  foimd,  then  either  CaAlk  and/or 
MnAlk  are  zero  and  the  solution  is  undersaturated  with  either  or  both  of  these  solids. 


5 

Time  Variable  Linear  Partitioning  Model 


5.1  INTRODUCTION 

This  chapter  presents  the  formulation  and  application  of  the  time  variable  version  of 
the  linear  partitioning  manganese  flux  model  presented  in  Chapter  3.  The  MERL  data 
set  is  used  for  calibration  and  comparison. 


5.2  MODEL  FORMULATION 

The  time  variable  manganese  flux  model  is  a  straightforward  extension  of  the  steady 
state  model  discussed  in  Chapter  3.  The  time  variation  introduces  an  additional 
feature  besides  the  need  to  consider  the  derivatives  in  eqs.(3.11a-3.11d).  It  is  related 
to  the  variation  in  the  thickness  of  the  aerobic  and  anaerobic  layers  and  the  mass 
balance  consequences. 


5.2.1  Entrainment  Flux 


The  time  variable  form  of  the  mass  balance  equations  can  be  derived  by  considering 
the  rate  of  change  of  the  mass  of  the  constituent  being  considered  in  each  layer.  For  a 
concentration  ci  and  a  depth  Hi  the  product  HiCi  is  the  mass/unit  area,  e.g.  mg/m^, 
in  layer  1.  The  mass  balance  equations  equate  the  rate  of  change  of  the  Tna.<;s  to  the 
sources  and  sinks 


djHiCi) 

dt 


Si 


(5.1a) 


d{H2C2) 

dt 


S2 


(5.1b) 


where  Si  and  S2  represent  all  the  source-sink  terms.  Consider  what  happens  if  the 
layer  thickness  changes  and  there  are  no  sources  or  sinks.  The  situation  is  illustrated 
in  Fig.  5.1.  Eqs.(5.1)  become 


dt 

(5.2a) 

d{H2C2)  _ 

dt 

(5.2b) 

53 

54  TIME  VARIABLE  LINEAR  PARTITIONING  MODEL 


t  t  +  At 

Fig.  5.1  Volume  changes  from  i  to  t  +  At. 


t  t  +  At 


Fig.  5.2  Entrainment  of  layer  2  mass  (shaded  region)  into  layer  1  during  the  time  step  At  as 
Hi  increases. 

and  the  solutions  are 

=  ffi(0)ci(0)  (5.3a) 

H2{t)c2(t)  =  H2{0)c2{0)  (5.3b) 

That  is,  the  mass  in  each  layer  remains  constant  over  the  time  interval  At.  Therefore 
the  volumes  are  not  coimected  in  any  way.  As  the  volume  expands  or  contracts,  the 
concentrations  decrease  or  increase,  respectively,  to  keep  the  mass  constant. 

However,  this  is  not  the  correct  result  if  the  depth  of  one  layer  is  changing  at  the 
expense  of  the  other.  In  the  sediment,  the  motion  of  the  H1—H2  boundary  causes  an 
entrainment  of  mass  into  and  out  of  the  layers.  This  is  illustrated  in  Fig.  5.2.Consider 
the  sediment  column  at  time  t.  Suppose  that  during  the  time  At  the  H1-H2  boundary 
moves  downward  as  illustrated.  For  this  case  Hi  is  increasing  and  dHi/dt  >  0.  Then 
the  shaded  layer  of  sediment,  which  was  in  H2  at  time  t,  moves  into  Hi  at  time  t  +  At. 
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This  amounts  to  a  loss  of  mass  from  H2  and  a  gain  for  Hi-  Thus  the  mass  balance 
equations  become 


d(HiCi) 

dt 


=  +C2- 


(5.4a) 


(5.4b) 


d{H2C2)  ^  dBj 

For  the  case  that  the  boimdary  moves  upward  and  dHi/dt  <  0  the  source  and  sink 
terms  are  reversed 


djHiCi) 

dt 

d{H2C2) 

dt 


=  -Cl 


=  +Ci 


=  -Ci- 


(5.5a) 


(5.5b) 


Mass  is  lost  from  top  layer  and  gained  by  the  bottom  layer.  The  absolute  value  signs 
are  employed  to  emphasize  that  the  amoimt  of  mass  entrained  is  determined  by  the 
magnitude  of  dHi/dt  and  not  its  sign. 

It  is  possible  to  write  eqs.(5.5)  as  one  set  of  equations  by  defining  the  functions 


Hr  = 


if  ^>0 

dt  dt 

o  T  dHi  ^ 


if  ^<0 

dt  dt 

0  if  ^>0 

at 


(5.6a) 


(5.6b) 


where  the  dot  denotes  the  time  derivative.  A  straightforward  way  of  computing  this 
function  is  to  use  the  magnitude  of  the  derivative  so  that  positive  and  negative  terms 
subtract  as  appropriate  to  give  zero^ . 


1  _  dHi 

2  \  dt  dt 


H-  Hi 


(5.7b) 


straightforward  implimentation  of  eq.(5.6)  involves  the  use  of  if()  statements.  For  compilers 
that  parallelize  computations  it  is  imperative  that  if  ()  statements  be  avoided.  The  abs()  function 
is  faster  and  accomplishes  the  same  thing. 
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Using  these  functions  the  mass  balance  equations  which  include  entrainment  can  be 
written  as 

d(HiCi) 

dt  -  -<=.«!- 

(5.8a) 

d(H2C2) 

(5.8b) 

and  employing  the  chain  rule  for  the  derivatives  yields 

„  dci  dHi  .  , 

dt  = 

(5.9a) 

„  dc2  dH2  •  , 

Cliff 

(5.9b) 

and  finally 

^  =  +C2H^  -  Cl  (Hi  +  H~ ^ 

(5.10a) 

H2^  =  -C2{H2  +  H+^+CiH{ 

(5.10b) 

where  Hi  =  dHidt  and  H2  =  dH2dt. 


5.2.2  Layer  Thickness 


The  time  derivative  of  the  aerobic  layer  thickness  Hi  is  computed  from  the  equation 
specifying  f/'i(eq.3.4)  and  from  substituting  s  =  S0Df02  (0) 


Hi 


Di  [O2  (0)1 
SOD 


Taking  the  derivative  of  this  equation  yields 


(5.11) 


dHi 

dt 


(A  [02(0)]) 

_ dt 


-(A  [02(0)]) 


50Z)2 


dSOD 

dt 


(5.12) 


1  /d(Oi[O2(0)])  dSOD^\ 

SOD  V  dt  dt  7 


(5.13) 


The  variation  in  H2  comes  from  the  assumption  of  a  constant  active  layer  depth  Ht 


Ht  =  Hi  +  H2 

Differentiating  with  respect  to  t  yields 

dt  dt 


(5.14) 


(5.15) 
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dH2  dHi 
dt  ~  dt 

which  simply  states  that  as  Hi  increases  H2  decreases  and  vise  versa. 


5.2.3  Model  Equations 

The  hnear  partitioning  model  as  described  in  Chapter  3  is  the  basis  for  the  time 
variable  version.  The  mass  balance  equations  for  the  model  are  based  on  eqs.  (3.11a- 
3.11d)  with  the  addition  of  the  entrainment  terms. 

Layer  1  Mn{II): 

„  dMn(l) 

Jt  =  (/diMn(l)  -  Mn  (0))  +  KLi2ifd2Mn{2)  -  fdiMn{l))  {5.17a) 


■^wi2{fp2Mn{2)  -  /piMn(l))  -  W2Mn  (1) - 

+Mn  (2)  H+  -  Mn  (1)  j 

Layer  2  Mn{II): 

dMn(2) 

It  ~  -KLi2{fd2Mn{2)  -  /diMn(l))  (5.17b) 

-wi2(fp2Mn(2)  -  fpiMn(l))  -  W2  (Mn(2)  -  Mn(l)) 
+i^Mn,2Mn02(2)  -  Mn{2)  +  Mn{l)H- 

Layer  1  Mn02  (s): 

^  -  — ^/diMn(l)  -  w2Mn02(l)  (5.17c) 

+Wi2(Mn02(2)  -  Mn02(l))  +  J^nOa 


Layer  2  Mn02  (s): 

dMn02(2) 

- 


+Mn02(2)ij+  -  Mn02(l)  (Hi  +  i/f  j 
-KMn,2Mn02i2) 


(5.17d) 


-u;i2(Mn02(2)  -  Mn02(l))  -  W2  {Mn02{2)  -  Mn02{l)) 
-Mn02{2)  (h2  +  +  Mn02{l)Hi 
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The  numerical  scheme  used  to  solve  these  equations  is  the  same  implicit  integration 
technique  that  is  used  for  the  nutrient  model  equations  [Di  Toro  and  Fitzpatrick, 
1993].  Since  there  are  four  equations  matrix  notation  will  be  used.  Equations  (5.17) 
are  written  as 

H(t)^  =  A(t)c(t)  +  b(t)  (5.18) 

where  the  bold  face  symbols  denote  vectors  (lower  case)  and  matrices  (upper  case). 
The  manganese  concentration  state  variable  vector  is 


Mn(l) 

Mn(2) 

Mn02(l) 

Mn02(2) 


(5.19) 


The  matrix  A{t)  =  [aij]  has  elements 

ail  =  -sfdi  -  /i2  -W2-  -^fdi  -  [Hi  +  iff)  (5.20a) 

ai2  =  hi  +  Ht  (5-20b) 

^21  =  /l2  +  '^2  (5.20c) 

022  =  -/2i-«;2- (^2  +  ^+)  (5.20d) 

024  =  X-Mnfi  (5.20e) 

031  =  fdl  (5.20f) 

s 

a^z  =  — ioi2  —W2—  [Hi  +  Hi  ^  (5.20g) 

034  =  1012  +  ^1+  (5-20h) 

043  =  W12  +  W2  +  Hi  (5.20i) 

O44  =  —KMn,2  —  '*012  —  102  —  (i?2  +  H^  ^  +  ijj  (5.20j) 


where  /12  and  /21  are  the  effective  mass  transfer  coefficients  from  layer  1  to  layer  2 
and  from  layer  2  to  layer  1  respectively 


fl2  =  KL12fdl  +  1012/pl 
/2I  =  KL12fd2  +  10l2/p2 


(5.21a) 

(5.21b) 
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and  all  other  Oy’s  are  zero.  The  forcing  function  is 

sMn  (0)  (t) 

0 

JMn02  (^) 

0 

and  H(t)  is  a  diagonal  matrix  of  layer  depths. 

H(t)  =  diag[  ffi(t)  H2(t)  Hiit)  H2{t)]  (5.23) 

A  finite  difference  equation  is  used  to  solve  eq.(5.18).  An  implicit  forward  in  time 
scheme  [Hamming,  1962]  is  used  for  the  concentrations  since  some  of  the  coefficients 
in  A  can  become  quite  large,  e.g.  if  s  is  small,  and  an  explicit  scheme 

would  require  an  unacceptably  small  At  to  remain  stable.  The  A  and  H  matrices 
are  evaluated  at  time  level  t  for  the  terms  coupled  to  the  other  state  variables  such 
as  s.  Otherwise,  the  equation  set  would  be  nonlinear  and  more  difficult  to  solve.  The 
resulting  finite  difference  equation  is 

which  can  be  factored 

(H(t)  -  AtA(t))  c(t  +  At)  =  H(t)c(t)  +  Atb(t)  (5.25) 

and  solved 

c(t  +  At)  =  (H(t)  -  AtA(t))-^  (H(t)c(t)  +  Atb(t))  (5.26) 

by  inverting  the  matrix  as  shown,  or,  using  a  more  efficient  linear  equation  solver 
[Press  et  al.,  1989]  applied  to  eq.(5.25).  It  requires  the  solution  of  four  simultaneous 
linear  equations  instead  of  the  two  required  for  the  other  variables  such  as  ammonia 
or  oxygen  [Di  Toro  and  Fitzpatrick,  1993]  because  the  equations  for  Mn{II)  and 
Mn02{s)  in  the  two  layers  are  coupled.  This  presents  no  essential  difficulty,  however, 
and  the  implementation  is  straightforward. 

5.3  REVIEW  OF  MERL  DATA 

The  MERL  mesocosms  are  large  outdoor  tanks  approximately  1.8  meters  in  diameter 
and  have  a  water  depth  of  5.0  meters.  Narragansett  Bay  water  flows  through  at  a  rate 
to  establish  a  detention  time  of  approximately  30  days.  The  tanks  have  a  mixer  to 
reproduce  the  vertical  mixing  regime  in  Narragansett  Bay  and,  in  particular,  the  lack 
of  stratification.  Sediments  are  obtained  using  a  large  box  cores  which  maintains  the 


(5.22) 
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vertical  orientation  and  the  top  40  cm  are  placed  into  containers  in  the  bottom  of  the 
tanks  [Nixon  et  al.,  1986]. 

The  data  to  be  analyzed  below  comes  from  the  Nutrient  Addition  Experiment.  Its 
purpose  was  to  examine  the  consequences  of  nutrient  enrichment  to  coastal  estuaries. 
The  duration  was  approximately  2  1/4  years  over  three  calendar  years.  The  nutrient 
dosing  was  increased  in  a  geometric  series,  IX,  2X,  4X,  8X,  16X,  and  32X,  in  addition 
to  three  control  tanks,  denoted  as  OX.  The  nutrients  added  were  inorganic  nitrogen, 
phosphorus,  and  silica  in  a  molar  ratio  of  12.8  N  :  1.0  P  :  0.91  Si  to  match  the 
stoichiometry  of  sewage  entering  the  bay  [Nixon  et  al.,  1986).  Areal  loading  rates  of 
total  nitrogen  to  the  tanks  varied  from  23  mg  N/m^-d  for  the  controls,  63  mg  N/m^-d 
for  IX,  103  mg  N/m^-d  for  2X  and  so  on  geometrically  to  1308  mg  N/m^-d  for  the  32X 
[Kelly  et  al.,  1985).  As  a  result,  mean  annual  water  column  dissolved  inorganic  nitrogen 
(DIN)  concentration  increased  from  56  to  4200  pg  N/L;  mean  annual  chlorophylla 
ranged  from  4  to  70  fig/L,  and  total  system  carbon  production  ranged  from  0.55  to 
2.2  g  C/m^-d  [Nixon  et  al.,  1986]. 

Sediment  processes  were  also  examined  during  the  experiment.  Sediment  oxygen 
demand,  nutrient  fluxes,  pore  water  and  solid  phase  concentrations  were  also  measured. 
In  addition,  manganese  flux  and  sediment  compositional  data  were  collected  [Hunt  and 
Kelly,  1988].  These  are  the  data  that  were  compared  to  the  steady  state  model  results 
in  Figs.3.4  and  3.7. 

Since  the  manganese  flux  model  is  coupled  to  the  nutrient  and  oxygen  flux  model— 
the  manganese  flux  model  requires  s  and  the  depth  of  the  aerobic  layer  frTi— the 
nutrient  and  oxygen  data  and  model  performance  is  evaluated  below. 

5.4  APPLICATION  OF  THE  NUTRIENT  AND  OXYGEN  FLUX  MODEL 

The  sediment  model  is  applied  in  stand  alone  mode  when  it  is  being  initially  calibrated. 
That  is,  it  is  not  dynamically  coupled  to  the  overlying  water,  as  it  would  be  if  it  were 
being  used  as  part  of  a  larger  coupled  water  quality-sediment  model.  Therefore  the 
overlying  water  concentrations  are  specified  externally.  In  addition,  the  average  annual 
depositional  flux  of  particulate  organic  nitrogen  to  the  sediment  Jpon  is  specified  and 
Redfield  stoichiometry  is  assumed  in  specifying  the  carbon,  phosphorus  and  silica 
fluxes. 

For  application  to  the  MERE  dataset,  the  magnitudes  of  JpoN  s.re  chosen  to  fit 
the  measured  ammonia  fluxes.  The  kinetic  and  transport  parameter  values  are  set  to 
the  values  obtained  from  the  calibration  to  the  Chesapeake  Bay  data  set  [Di  Toro  and 
Fitzpatrick,  1993].  The  required  depositional  fluxes  are  shown  in  Fig.5.3.  Fluxes  vary 
from  less  than  50  to  130  mg  N/m^-d,  less  than  a  3  fold  variation. 
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0  Year  2 
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Fig.  5.3  Depositional  flux  of  particulate  organic  nitrogen. 


In  order  to  run  the  model  it  is  necessary  to  specify  the  initial  state  of  the  sediment, 
i.e.  the  concentrations  of  all  the  state  variables.  Since  these  concentrations  were  not 
measured  ^indeed  it  would  be  diflScult  if  not  impossible  to  measure  the  concentrations 
in  the  aerobic  layer— they  must  be  arrived  at  in  another  way.  The  objective  is  to 
produce  the  initial  conditions  in  the  sediments  at  the  start  of  the  experiment.  The 
method  chosen  is  to  equilibrate  the  model  by  cychng  the  model  until  a  periodic  steady 
state  is  achieved.  The  same  forcing  functions  are  used  for  each  annual  cycle.  Since 
the  objective  is  to  simulate  the  state  of  the  sediments  at  the  start  of  the  experiment, 
the  obvious  choice  is  the  inputs  for  the  first  year  from  the  control  tanks.  That  is,' 
the  overlying  water  concentrations  observed  in  a  control  tank  and  the  depositional 
flux  which  reproduces  the  ammonia  flux  for  that  tank  are  used  to  initialize  all  the 
tanks.  This  is  equivalent  to  assuming  that  the  state  of  the  sediments  in  the  year  when 
they  were  collected  is  similar  to  that  which  is  produced  by  a  cyclical  repetition  of  the 
conditions  which  were  observed  in  the  control  tank  for  the  first  year  of  the  experiment. 
Since  the  controls  are  a  reasonable  replication  of  Narragansett  Bay,  this  seems  not  to 
be  an  unreasonable  assumption. 

It  is  also  a  very  convenient  and  parsimonious  way  of  proceeding.  The  initial  condi¬ 
tions  are  not  externally  specified  by  some  arbitrary  means.  Therefore,  they  are  not  an 
additional  set  of  calibration  parameters  that  need  to  be  determined  by  fitting  to  the 
observations.  This  greatly  reduces  the  degrees  of  freedom  in  the  model  and  renders 
the  calibration  less  of  a  curve  fitting  exercise.  Also  the  initial  conditions  are  consis- 
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tent  with  all  the  parameters  used  in  the  model.  Thus,  initially  all  the  sediments  have 
the  same  concentrations  as  the  control.  The  differences  in  subsequent  behavior  occur 
solely  due  to  the  increased  nutrient  loading  applied  during  the  experiment. 

The  resulting  ammonia  fluxes  are  shown  in  Fig.5.4.  The  seasonal  variation  is  reason¬ 
ably  well  reproduced  except  for  32X.  However,  the  negative  fluxes  during  the  winter, 
which  increase  as  loading  increases,  are  not  captured,  and  we  have  no  good  explanation 
for  these  observations. 

The  oxygen  fluxes  are  shown  in  Fig.5.5.  The  results  are  remarkably  good.  The 
absolute  magnitudes  and  the  seasonal  variations  are  reproduced  up  to  32X.  Since  the 
manganese  model  requires  the  surface  mass  transfer  coefficient  s  =  SOD/[O2{0)],  it 
is  important  that  the  oxygen  flux  be  reproduced  by  the  model. 

5.5  MANGANESE  MODEL  APPLICATION 

The  manganese  flux  model  is  run  as  follows.  The  nutrient  portion  is  run  as  described 
above.  The  overlying  water  concentrations  of  dissolved  manganese  are  shown  in  Fig. 
5.8  for  all  the  tanks.  For  the  equilibration  phase  the  dissolved  manganese  concentra¬ 
tions  for  all  the  tanks  are  set  to  the  average  of  the  overlying  water  concentration  of 
the  control  tanks,  which  is  shown  in  the  top  two  panels  of  Fig.5.8.  These  conditions 
are  cycled  periodically  until  the  sediment  has  come  into  equilibrium.  For  the  par¬ 
ticulate  manganese,  the  sediments  are  equilibrated  using  a  constant  depositional  flux 
which  is  calculated  from  a  water  column  settling  velocity  and  the  average  particulate 
concentrations,  shown  in  Fig.  5.9  for  all  the  tanks.  The  overlying  water  particulate 
concentration  is  assumed  to  be  constant  as  shown. 

Once  equilibrium  has  been  reached,  the  nutrient  loadings  are  increased  in  accor¬ 
dance  with  the  various  loading  levels  and  the  model  is  run  for  the  three  years  of  the 
experiment.  The  overlying  water  dissolved  manganese  concentrations  are  set  at  the 
observed  values  for  the  three  years  (Fig.  5.8).  The  particulate  manganese  concentra¬ 
tions  are  continued  at  the  same  constant  value  as  used  during  the  equilibration  phase 
(Fig.  5.9). 

Thus,  it  is  assumed  that  water  column  processes  during  the  experiment  do  not 
change  the  depositional  flux  of  manganese  to  the  sediment.  This  is  only  an  approxi¬ 
mation  since  the  overlying  water  concentrations  vary  as  dosing  varies.  However,  overall 
the  variation  is  not  too  large  so  that  the  assumption  of  a  constant  flux  from  the  over- 
lying  water  appears  to  be  a  reasonable  first  approximation.  The  consequences  of  this 
assumption  is  evaluated  below  when  the  overlying  water  is  explicitly  included  in  the 
model. 
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A  comparison  of  the  model  and  observed  fluxes  is  shown  in  Fig.  5.10.  The  param¬ 
eters  are  listed  in  Tables  5.1  and  5.2.  The  general  shape  of  the  seasonal  variation  of 
the  manganese  fluxes  is  reproduced  quite  nicely.  This  variation  is  due  solely  to  the 
variation  in  s  since  the  flux  of  particulate  manganese  to  the  sediment  is  held  constant 
(Fig.  5.9).  This  is  one  reason  for  using  the  constant  manganese  depositional  flux. 

The  variation  in  annual  average  manganese  flux  due  to  differing  loading  rates  is 
examined  in  Fig.  5.6.  The  model  captures  the  trend  of  increasing  flux  as  nutrient 
loading  increases  for  the  first  (Fig.  5.6A)  and  third  (Fig.  5.6C)  years.  However 
the  sharp  increases  in  the  second  year  for  the  OX  to  4X  tanks  (Fig.  5.6B)  are  not 
reproduced. 


20 


Fig.  5.6  Average  annual  manganese  flux  for  the  three  years  of  the  experiment.  Comparison 
of  model  and  data. 

The  variation  of  manganese  fluxes  with  respect  to  temperature,  overl3dng  water 
dissolved  oxygen,  and  ammonia  flux  for  both  the  data  and  the  model  are  presented  in 
FiS-  5-11.  The  data  exhibit  a  more  pronounced  variation  with  respect  to  temperature 
and  ammonia  flux  than  does  the  model.  This  is  consistent  with  the  findings  of  the 
steady  state  model  (Fig.  3.7). 

The  pore  water  and  solid  phase  concentrations  are  compared  in  Fig.5.7.  The  model 
concentrations  are  shghtly  higher  than  the  observations,  as  was  the  case  with  the 
steady  state  model  (Fig.  3.7). 

5.6  CONCLUSIONS 

The  time  variable  version  of  the  model  captures  many  of  the  features  exhibited  by 
the  data:  a  strong  seasonal  variation  (Fig.5.10);  an  increasing  manganese  flux  with 
increasing  nutrient  inputs  (Fig.5.6);  and  a  slight  depletion  of  the  sediment  manganese 
concentration,  as  a  consequence  of  the  increased  manganese  flux  (Fig.5.7).  However, 
the  modeled  relationship  between  manganese  and  ammonia  fluxes  is  still  less  than 
observed  (Fig.5.11),  and  the  inter-annual  variations  are  not  quite  resolved. 

The  problem  may  be  that  the  pH  variation  exhibited  in  the  MERL  experiments, 
and  which  is  known  to  significantly  affect  the  oxidation  rate  of  manganese  (II)  to 
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Fig.  5. 7  Comparison  of  pore  water  and  solid  phase  concentrations. 


manganese  (IV),  has  not  been  taken  into  account.  This  is  addressed  in  the  next 
section. 
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Fig.  5.8  Dissolved  manganese  concentrations  in  the  overlying  water 


Fig.  5.9  Particulate  concentrations  in  the  overlying  water.  Line  is  the  concentration  used 
the  model  for  the  depositional  flux. 
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Fig.  5.10  Comparison  of  observed  and  modeled  manganese  fluxes. 
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Fig.  5.11  Comparison  of  modeled  and  observed  fluxes  versus  temperature,  overlAdng  water  O2 
concentration,  and  ammonia  flux  J  [7V/f4] . 
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Table  5.1  Sediment  Model  Parameters  for  MERL  Application 


mi 

1712 

W2 

H2 

Hq 

to 

(kg/L) 

(kg/L) 

(m/d) 

(m) 

(m) 

d 

0.5 

0.5 

6.85  X  10"® 

0.10 

5. 

27.4 

Dd 

Dp 

Ddo 

(mVd) 

- 

(mVd) 

(mVd) 

5.0E>03 

1.08 

1.2  X  10-4 

1.15 

1.0  X  10-4 

1.08 

^NHa. 

Km,NHa 

^7Cm,NH4 

Km,NHa,02 

(m/d) 

(L/kg) 

- 

(mg  N/L) 

- 

(mg  O2/L) 

0.131 

1.00 

1.123 

0.728 

1.125 

0.74 

K'N03,2 

^NOa 

(m/d) 

(m/d) 

- 

0.10 

0.25 

1.080 

T^HzS,! 

^H2S 

Km,H2S,02 

(na/d) 

(m/d) 

(L/kg) 

(L/kg) 

- 

(mg  02/Lj 

0.2 

0.4 

100.0 

100.0 

1.0800 

4.0 

fcsi 

[Si]  sat 

Attsi,! 

’’■31,2 

^Si 

Km,P  Si 

(/d) 

(mg  Si/L) 

- 

(L/kg) 

- 

(mg  Si/L) 

0.5 

40. 

10.0 

100.0 

1.100 

5.0  X  104 

A7rpc>4,i 

‘^POi,2 

^2  (0);r,PO4 

Km,Dp,02 

C>2  (0),,si 

— 

(L/kg) 

(mg  O2/L) 

(mg  O2/L) 

(mg  O2/L) 

20.0 

20.0 

2.0 

4.0 

1.0 

ATTA/n.l 

T^Mn,2 

^2  (^)n,Mn 

- 

(L/kg) 

(mg  O2/L) 

3. 

700. 

2.0 

^Afn02,l 

^Atn02,l 

I^Mn02,0H,l 

^Mn02,0H,\ 

1^Mn022 

^Mn022 

(m/d) 

- 

(m/d) 

- 

(m/d) 

— 

0.0 

1.08 

0.01 

1.08 

0.1 

1.08 

Alko 

pH2 

Alk2 

/colloidal 

Ks 

-Dp, min 

(meq/L) 

- 

(meq/L) 

- 

(/d) 

(mVd) 

2.0 

7.5 

3.0 

0.5 

0.03 

3  X  10-® 

'^Mn02iin 

^Jn 

(m/d) 

- 

2.0 

1.15 
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Table  5.2  Diagenesis  Parjimeters  for  MERL  Application 


fpONi 

kpONi 

GpONi 

fpONa 

kpONi 

^PONa 

— 

(/d) 

- 

- 

i/d) 

0.65 

0.035 

1.10 

0.25 

0.0018 

1.15 

fpoci 

kpoci 

GpOCi 

fpOCi 

kpoCi 

^POCi 

— 

(/d) 

- 

- 

i/d) 

0.65 

0.035 

1.10 

0.20 

0.0018 

1.15 

fPOPi 

kpOPi 

OpoPi 

fpOPa 

kpOPa 

^POPi 

(/d) 

- 

- 

i/d) 

_ 

0.65 

0.035 

1.10 

0.20 

0.0018 

1.15 

fpONs 

kpONs 

^PONa 

- 

(/d) 

- 

0.10 

0.0 

1.17 

fpOCs 

kpOCs 

^POCs 

- 

i/d) 

- 

0.15 

0.0 

1.17 

fpOPs 

kpOPz 

^POPa 

- 

(/d) 

- 

0.15 

0.0 

1.17 

0'C,N 

oc.p 

0■C,Si 

(mg  C/mg  N) 

(mg  C/mg  P) 

(mg  C/mg  Si) 

5.68 

41 

2.5 
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Effect  of  pH 


6.1  INTRODUCTION 

The  purpose  of  this  chapter  is  to  examine  the  effect  of  pH  variations  on  the  flux  of 
manganese  from  the  sediment.  The  rate  of  the  oxidation  reactions  for  both  manganese 
and  iron  are  strongly  affected  by  the  pH  (Fig.2.1-2.4)  so  it  is  reasonable  to  expect  that 
this  effect  is  important.  Since  the  oxidation  reaction  takes  place  in  the  aerobic  layer 
of  the  sediment,  this  is  where  the  pH  is  required.  Since  the  layer  1  pH  depends  on  the 
conditions  in  the  overlying  water  and  in  the  anaerobic  layer,  and  on  the  magnitudes 
of  the  mass  transfer  coefficients  that  affect  the  aerobic  layer,  it  must  be  computed  as 
part  of  the  model. 

6.2  pH  IN  THE  AEROBIC  LAYER 

The  pH  in  the  aerobic  layer  (layer  1)  is  determined  by  the  layer  1  alkalinity  Alkiand 
inorganic  carbon  TICi  concentrations.  These,  in  turn,  are  determined  by  the  mass 
fluxes  of  alkalinity  and  inorganic  carbon  from  the  adjacent  layers  and  internal  sources 
and  sinks.  For  the  sake  of  simplicity,  we  assiune  that  neither  internal  sources  nor  .sinks 
oi  Aik  and  TIC  are  significant.  Any  other  assumption  would  require  a  consideration  of 
the  acid-base  reactions  in  the  layer  as  well  as  any  TIC  and  alkalinity  sources.  Because 
the  layer  is  quite  thin,  it  is  hkely  that  their  effects  are  not  large.  Therefore  alkalinity 
and  TIC  are  assmned  to  be  conservative  variables  in  layer  1.  The  schematic  diagram 
is  presented  in  Fig.  6.1. 

The  mass  balance  equations  that  result  for  this  model  are 

s  (Alko  -  Alki)  -f  Kli2  {Alk2  -  Alk{)  =  0  (6.1) 

and 

s  {TICo  -  TICi)  +  Kli2  {TIC2  -  TICi)  =  0  (6.2) 

where  s  and  K112  are  the  mass  transfer  coefficients  for  layer  0-1  and  layer  1-2  respec¬ 
tively  and  subscripts  are  used  to  identify  the  layer.  Solving  these  eqs.(6.1-6.2)  for  the 
layer  1  concentrations  yields 


(6.3) 


Alki  =  aAlko  -k  (1  —  a)  Alk2 
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Fig.  6.1  Schematic  diagram  of  the  model  for  the  aerobic  layer  1  pH. 


and 


TICi  =  aTICo  +  ( 1  -  a)  TIC2  (6.4) 


where 


s 

a  = - 

s  +  K112 

Once  Alki  and  TICi  are  known,  the  pH  is  computed  using  the  equations 
the  carbon  dioxide  system  [Stumm  and  Morgan,  1996]. 

The  mass  balance  equations  for  the  carbon  dioxide  system  are 


(6.5) 
that  describe 


[TIC\  =  ICO2I  +  [HCO;]  +  [COf-] 

(6.6) 

\Mk]  =  [HCOl]  +  2  [col-] 

(6.7) 

and  the  mass  action  equations  are 

[HCOi]  p+l 

ICO2I  - 

(6.8) 

[col]  {H*\ 

[J/CO3-]  " 

(6.9) 

Eliminating  [iy+j,  and  [CO^~]  yields 

K1K2  {Aik  -  2TIC)  +  KiH  {Aik  -  TIC)  +  H^Alk  =  0 

(6.10) 

which  relates  the  pH  to  the  alkalinity  and  inorganic  carbon  concentrations.  The  square 
brackets  and  charges  are  dropped  from  the  notation  for  convenience. 

The  aerobic  layer  alkalinity  and  TIC  concentrations  are  given  by  eqs.(6.3-6.4).  For 
the  alkalinity,  both  AlkQ  and  Alk2  are  known  by  measurement.  However,  the  TIC’s 
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Fig.  6.2  Aerobic  layer  pHi  as  a  function  of  overlying  water  pHo.  Assiuned  conditions:  anaer¬ 
obic  layer  pH2  =  7.5;  alkalinity  Alki  =  3  meq/L;  and  layer  1-2  mass  transfer  coefiBicient 
Kli2  =  0.01  m/d.  Remaining  parameters:  s  and  Alko,  are  noted  in  the  figure. 


are  not  measured.  They  can  be  computed  from  the  respective  pHs  and  alkalinities. 
Using  eq.(6.10)  to  solve  for  TIC  3delds 


TIC  = 


Alk{Ki{K2  +  H)  +  H^) 
Ki{2K2  +  H) 


Using  the  relationship  between  pH  and  hydrogen  ion  concentration  yields 


(6.11) 


Hq  =  io-p^° 


(6.12) 


which,  together  with  AlkQ,  determine  TICq  using  eq.(6.11).  A  similar  procedure  is 
applied  to  layer  2  to  obtain  TIC2  using  pH2  and  Alk2.  Then  eq.(6.4)  yields  TICi. 
Since  the  layer  1  alkaUnity  and  TIC  are  now  known,  eq.(6.10)  can  be  solved  for  Hi 
to  obtain  the  pHi.  The  result  is 


Ki  {TICi  -  Alki)  ±  ^JiAlki  {2TICi  -  Alki)  K1K2  +  {TICi  -  Alkif  Kl 
^  2^1 

(6.13) 


and 


pHi  =  -\og^o[Hi\ 


(6.14) 


An  example  computation  is  given  in  Fig.6.2  which  presents  the  pH  in  layer  1  as  a 
function  of  the  pH  in  the  overlying  water,  for  various  overlying  water  alkalinities,  for 
two  values  of  the  surface  mass  transfer  coefficient  s. 

It  is  useful  to  understand  what  is  determining  pHi.  A  part  of  the  behavior  has  a 
straightforward  rationalization.  For  a  large  surface  mass  transfer  coefficient  s  =  0.1 
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m/d  the  aerobic  layer  pH i  tracks  the  overlying  water  pHo  since  the  mixing  transports 
both  the  overlying  water  alkalinity  and  TIC  to  the  aerobic  layer.  For  the  smaller 
surface  mass  transfer  coefficient  s  =  0.01  m/d  the  aerobic  layer  pHi  remains  close  to 
the  anaerobic  layer  pH2  =  7.5  since  it  is  more  isolated  from  the  aerobic  layer. 

What  is  puzzling,  however,  is  why  the  aerobic  layer  pHi  varies  as  it  does  with 
varying  overlying  water  alkalinity  Aik  (0).  Increasing  the  alkalinity  from  Aik  (0)  =  1 
meq/L  to  2  meq/L  brings  the  overlying  water  alkalinity  closer  to  the  alkalinity  of  the 
anaerobic  layer  Aik  (2)  =  3  meq/L.  One  would  expect  that  making  the  alkalinity  in 
the  overlying  water  closer  to  the  alkalinity  of  the  anaerobic  layer  would  bring  the  pH 
in  the  aerobic  layer  closer  to  the  anaerobic  layer  pH  as  well.  But  the  reverse  happens. 
The  aerobic  layer  pHi  becomes  more  like  the  overlying  water  pHo. 

The  key  to  understanding  this  result  in  particular,  and  the  CO2  system  in  general, 
is  to  form  the  ratio  of  eq.(6.7)  to  eq.(6.6) 


[Aik] 

[HCO^]+2[COl-] 

(6.15a) 

[TIC] 

[cd2]  +  [Hco^]  +  [cdi-] 

(6.15b) 

(■•S) 

(6.15c) 

so  that  the  pH  is  determined  by  the  ratio  of  alkalinity  to  inorganic  carbon  [Aik]  /  [TIC] 
and  not  the  individual  values. 

This  relationship  is  shown  in  Fig.  6.3. 

Note  that  small  changes  in  /  [TIC]  result  in  large  changes  in  pH  due  to  the 
lack  of  buffering  capacity  in  the  pH  =  7-9  range  [Stumm  and  Morgan,  1996].  As  pH 
approaches  pKi  or  pK2  biaffering  is  increased  and  the  relationship  flattens  out. 

The  role  of  [A/A;]  /  [T/C]  in  determining  the  pH  is  illustrated  in  Fig.  6.4.  The  two 
panels  represent  the  alkalinity,  TIC,  and  [Aik]  /  [TIC]  ratio  for  a  change  in  overlying 
water  alkalinity  from  Aik  (0)  =  1  meq/L  to  2  meq/L.  An  extreme  case,  pHo  =  10,  is 
chosen  corresponding  to  the  right  hand  side  of  Fig.  6.2  to  amplify  the  differences. 

The  flrst  two  groups  of  bars  present  the  concentration  of  Aik  and  TIC  in  the 
three  layers.  The  center  bars  present  the  alkalinity  and  TIC  concentrations,  A/A:i  and 
TICi,  computed  from  the  adjacent  A/A:o  and  Alk2,  and  TICq  and  TIC2  concentrations 
respectively  using  eqs.(6.3-6.4).  For  the  third  group  of  bars,  the  unfilled  bars  present 
Alko/TICo  and  Alk2/TIC2  respectively.  Note  that  they  are  the  same  for  the  two  cases 
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Fig.  6.3  pH  as  a  function  of  [j4ZA:]  /  [T/C] .  Parameters:  Ki  =  10”®  ®,  K2  =  10”^°  ®. 


since  pHo  and  pffj  are  held  constant  and  there  is  a  one-to-one  relationship  between 
pH  and  Alk/TIC.  The  filled  bars  present  Alki/TICi  for  the  two  cases.  Increasing 
AlkQ  results  in  an  increase  in  Alki/TICi  and,  therefore  an  increase  in  pHi  as  shown 
in  Fig.  6.2.  The  difference  between  the  two  cases  is  small  but  because  pH  is  such  a 
sensitive  function  of  Alk/TIC  the  pH  varies  appreciably. 

Since  the  graphical  presentation  does  not  clearly  indicate  why  the  behavior  is  oc¬ 
curring,  perhaps  the  algebraic  relationships  can  provide  the  rationale.  The  equations 
that  determine  the  result  are 

Alkx  __  aAlko  -I-  (1  —  a)  Alk2 
TI^i  ~  aTICQ  +  {l-oc)TIC2 


(6.16) 


aAlko  +  (1  —  Q;)  Alk2 

arQ^Alko  +  (1  -  a)  r2^Alk2 


(6.17) 


where  tq  =  Alko/TICo  and  r2  =  Alk2/TIC2.  Factoring  this  expression  and  grouping 
the  terms  yields 

Alki  _  1  +  13 


TICi 


where 


■^1.12  Alk2 


The  limiting  cases  are 


Alkx  _  _  Alko 

tICx  0^0  ~  tIc^ 


(6.18) 


(6.19) 


(6.20) 


pHx  r*pHo 


(6.21) 
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Aik  TIC  Alk/TIC  Aik  TIC  Alk/TIC 

O.U  0.1.2  0.1.2  0.1.2  0.1.2  0.1.2 


Fig.  6.4  The  effect  of  overlying  water  alkalinity.  Example  computation  for  pHo  =  10,  pH2  = 
7.5,  Alk2  =  3.0  meq/L,  s  =  Klii  =  0.01  m/d.  Beirs  represent  Aik,  TIC,  and  Alk/TIC  for  layers 
0,1,2.  The  filled  bar  represents  Alk/TIC  for  layer  1  which  determines  pHi. 


Alternately 


Alki  _  Alk2 

TICi  ^-00  “  TJC2 

(6.22) 

pHi  pH2 

/3— ►©© 

(6.23) 

Thus  the  behavior  of  the  aerobic  layer  pH  is  determined  by  the  magnitude  of  /?.  A 
large  /?  can  occur  if  K112  is  large  relative  to  s  so  that  layer  2  is  transported  into  layer 
1  and  pHi  pH2-  This  is  obvious  case. 

Alternately  a  large  Alk2  relative  to  Alko  aJso  increases  p.  What  is  happening  in 
this  case  is  more  subtle  and  explains  the  results  in  Fig.6.2.  Decreasing  Alko  causes  a 
decrease  in  TICq  since  pHo  is  assumed  to  be  fixed.  Hence  more  TIC2  is  transported 
into  layer  1,  making  layer  1  more  hke  layer  2  so  that  again  pHi  pH2- 

What  prompted  this  exercise  was  the  need  to  understand  the  results  in  Fig.  6.2: 
namely  that  increasing  Alko  —*  Alk2  caused  the  pHi  to  more  closely  resemble  pHo- 
The  reason  is  that  increasing  Alko  causes  an  increase  in  TICo  which  makes  TIC\  more 
hke  TICo  s^nd  therefore  Alki/TICi  more  like  Alko/TICo  and  hence  pHi  more  like 
pHo. 


6.3  EFFECT  OF  pH  ON  THE  MANGANESE  FLUX 

The  pH  in  the  aerobic  layer  can  now  be  calculated  using  the  equations  developed  in 
the  previous  section,  eqs.(6.3-6.5)  and  (6.13).  The  results  are  shown  in  Fig.  6.5.  The 
solid  line  is  an  interpolation  of  the  data  points  in  the  overlying  water,  which  are  also 
shown.  The  dashed  line  is  the  computed  pH  in  the  aerobic  layer.  The  dotted  line 
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is  at  the  pH  of  the  anaerobic  layer,  which  is  assumed  to  be  at  pH  —  7.5,  consistent 
with  pore  water  observations.  The  alkalinities  are  Alka  =  2.0  meq/L  and  Alk^  —  3.0 
meq/L,  which  are  also  estimated  from  observations.  The  remaining  parameters:  s  and 
Klvi  are  from  the  simulation  to  be  presented  below. 

For  the  controls,  IX,  and  2X  loadings,  there  is  little  computed  variation  in  pHi .  The 
reason  is  that  surface  mass  transfer  coefficient  s  =  SOD/O2  (0)  is  small,  corresponding 
to  a  comparatively  small  SOD  (Fig.5.5),  so  the  aerobic  layer  resembles  the  anaerobic 
layer.  However,  for  the  more  heavily  loaded  tanks,  the  SOD  increases,  s  increases 
and  the  aerobic  layer  starts  to  be  influenced  by  the  overlying  water  pH,  which  is 
varying.  Hence  the  variation  in  pHi  is  more  substantial  for  these  tanks,  sometimes 
approaching  one-half  pH  imit,  which  would  correspond  to  one  order  of  magnitude  in 
oxidation  reaction  rate  variation,  as  shown  subsequently.  It  is  for  this  reason — the 
expected  importance  of  the  pH  variation  in  the  aerobic  layer — ^that  this  refinement  is 
added  to  the  model. 

6.3.1  Reaction  Kinetics 

The  incorporation  of  pH  dependent  reaction  kinetics  is  straightforward 
layer  reaction  rate 

^2 

_J^O2(0)/diMn(l) 
s 

is  replaced  by  a  pH  dependent  expression 

_  [OH-f02  (0)  fdiMn{l)  (6.25) 

0 

where  the  dependency  on  [OH~]  is  from  an  analysis  of  experimental  results  [Morgan, 
1967].  This  is  expressed  as  follows 

(l02(pHi-7)^  (0)  /diMn(l)  (6.26) 

The  pH  is  normalized  to  pH  =  7  so  that  the  numerical  value  of  the  reaction  velocity 
corresponds  to  pHi  =  7.  This  procedure  is  useful  so  that  the  magnitude  of 
has  some  specific  meaning.  The  +2  in  the  exponent  results  for  the  dependency 
on  [OH-f 

=  (i^)'  =  (l^)'  =  (6.2T) 

where  Kw  =  [H"'']  [OH~]  ,the  ionization  constant  for  water.  With  the  kinetics  mod¬ 
ified,  it  remains  only  to  rerun  the  model  using  the  computed  pHi  and  examine  the 
results. 


.  The  aerobic 
(6.24) 
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pH  in  HO.  HI.  and  H2 
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Fig.  6.5  pH  in  the  overlying  water  (solid  line  is  interpolated  to  the  data  points),  aerobic  layer 
(computed  -  dashed  line),  and  anaerobic  layer  (assumed  constant  -  dotted  hne). 
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6.3.2  Application  to  MERL 

The  consequences  of  this  modification  are  shown  in  Figs.  (6.6-6.7).  There  are  almost 
no  discernible  differences.  Perhaps  this  should  not  have  been  too  surprising  if  one 
considers  the  steady  state  solutions  for  manganese  flux  equations:  eqs.(3.21-3.22  and 
3.26)  in  section  3.1.4  and  eq.(3.28)  in  section  3.I.4.I.  These  highlight  the  importance 
of  s  and  /^land  not  the  oxidation  rate  KMn,i  except  in  the  case  that  s  is  small  (eq.3.33). 

Actually  the  pH  is  important,  as  we  shall  see  in  the  next  chapter.  But  it  is  important 
in  the  overlying  water  where  the  variation  is  larger  and  so  also  is  its  effect  on  the 
oxidation  rate. 

6.4  EFFECT  OF  VARYING  DEPOSITIONAL  FLUX 

There  is  one  more  source  of  variation  that  has  not  been  taken  into  account  yet:  the 
variation  in  the  overlying  water  particulate  manganese  concentrations.  The  flux  of 
manganese  to  the  sediment  is 


JMn02  —  ^0  [Mtip  (0)]  (6.28) 

where  wq  is  the  settflng  velocity  from  the  overlying  water  to  the  sediment  and  [Mup  (0)] 
is  the  concentration  of  particulate  manganese  in  the  overlying  water.  Fig.  (6.8) 
presents  [Mnp(O)]  and  the  interpolated  function  used  to  compute  the  depositional 
flux  (eq.  6.28).  The  results  are  shown  in  Figs.  (6.9  and  6.10).  Again  the  differences 
axe  small,  suggesting  that  the  time  variation  of  [Mnp  (0)]  is  not  critical. 

But,  as  we  shall  see  in  the  next  chapter,  the  magnitude  of  the  particulate  flux  JMn02 
is  of  critical  importance.  The  question  to  be  addressed  next  is:  how  is  JMn02  be 
specified  in  model  applications  to  natural  bodies  of  water.  Presumably,  through  the 
use  of  eq.(6.28).  So  the  question  becomes:  how  is  [Mnp(O)]  to  be  computed.  We 
examine  this  in  the  next  chapter. 


JIMnI  (mg  Mn/m  -d)  ^  JlMnl  (mg  Mn/m  -d)  JlMnl  (mg  Mn/m^~d) 
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02(0)  (mg/L)  02(0)  (mg/L) 


Fig,  6.7  Compaxison  of  modeled  and  obsdrved  fluxes  versus  temperature,  overly  water  O2 
concentration  and  ammonia  flux  J[NH4]. 


Pirt.  Mn  (mg/L)  Part.  Mn  (mg/L)  Part,  Mn  (mg/L)  Part.  Mn  (mg/L) 


EFFECT  OF  pH 


Fig.  6.9  Manganese  flux  to  the  overlying  water.  Run  with  pH  dependent  oxidation  rate 
the  aerobic  layer  and  observed  overlying  water  particulate  naanganese  concentrations. 


JlMnl  (mg  Mn/m^-d)  JIMnI  (mg  Mn/m^-d)  J(Mn)  (mg  Mn/m^-d) 
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Manganese  Flux 
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Fig.  6.10  Manganese  flux  versus  T,  O2  (0) ,  and  J[NH4] .  Run  with  pH  dependent  oxidation 
rate  in  the  aerobic  layer  and  observed  overlying  water  particulate  manganese  concentrations. 
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Water  Column  and  Sediment  Interactions 

7.1  INTRODUCTION 

The  flux  of  dissolved  manganese  from  the  sediment  to  the  overlying  water  is  fueled 
by  the  depositional  flux  of  particulate  manganese  from  the  overlying  water  to  the 
sediment.  In  the  previous  chapters  this  flux  JMnOi  was  either  specifled  as  a  constant 
input  parameter  (Chapters  2-5)  or  computed  using  a  settling  velocity  wq  applied  to  the 
observed  overlying  water  particulate  concentration  Mup  (Chapter  6).  This  is  a  useful 
approach  for  model  development  and  calibration.  However,  for  actual  applications  it 
is  necessary  to  be  able  to  predict  the  water  column  manganese  concentrations  as  well. 

This  chapter  present  an  extension  of  the  manganese  and  iron  flux  model  which 
includes  the  overlying  water  as  an  integral  part  of  the  model.  This  formulation  is  the 
prototype  for  the  actual  coupled  water  colmnn-sediment  model  that  is  applied  to  two 
lakes  in  the  two  subsequent  chapters. 

The  addition  of  a  third  layer  to  the  model  is  conceptually  and  computationally 
straightforward.  Since  the  appHcation  is  initially  to  the  MERL  data  set,  the  water 
column  layer  is  modeled  as  one  completely  mixed  volume.  The  transport  terms  for 
this  layer,  in  addition  to  those  from  the  aerobic  layer,  is  simply  the  advective  inflow 
and  outflow. 

It  is  in  the  kinetic  formulation  that  some  choices  need  to  be  made.  The  question  is: 
what  controls  the  speciation  between  dissolved  and  particulate  manganese  in  the  water 
column.  For  the  sediment  aerobic  layer,  two  reactions  are  considered:  partitioning 
of  Mn(II)  to  sediment  solids,  and  the  oxidation  reaction  which  converts  Mn(II) 
to  Mn02  (5) .  Clearly  the  oxidation  reaction  is  important  since  the  overlying  water 
is  usually  oxic  and  that  reaction  occurs  fairly  rapidly.  The  question  is:  does  the 
partitioning  reaction  need  to  be  considered?  The  question  is  examined  next. 

7.1.1  Partitioning  of  Manganese  (II) 

The  principle  species  of  manganese  (II)  in  water  are  the  divalent  cation  Mn^+  and 
its  hydrolysis  product  Mn(PH)+  (see  Chapter  2  and  [Lindsay,  1979]).  One  expects, 
based  on  the  partitioning  behavior  of  other  divalent  metal  cations,  that  manganese 
(II)  should  partition  to  particles  in  general  and  to  hydrous  iron  oxides  in  particular. 
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In  principle  there  is  no  difficulty  in  allowing  both  the  partition  reaction 


Mud  •<->  Mup 


(7.1) 


and  the  oxidation  reaction 

Mn  {II)  ->  Mn02  (s)  (7.2) 


to  occur  simultaneously  in  the  water  column.  After  all,  this  is  exactly  the  formulation 
that  is  used  in  the  sediment  aerobic  layer.  The  issue  is:  can  the  partition  coefficient  be 
determined  in  some  independent  way,  or  can  it  be  shown  that  the  partitioning  reac¬ 
tion  is  insignificant  relative  to  the  oxidation  reaction  in  the  production  of  particulate 
manganese. 


An  estimate  of  the  expected  partition  coefficients  can  be  made  using  a  linear  free 
energy  relationship  [Brezonik,  1994].  The  idea  is  that  the  variation  in  the  equilibrium 
constant  for  a  series  of  chemical  species  reacting  with  one  ligand  can  be  used  to  predict 
the  variation  in  the  equilibrium  constant  for  another  ligand.  The  usual  choice  for 
sorption  reactions  is  the  equilibrium  constant  for  metal  hydrolysis  [Baes  and  Mesmer, 
1976] 


Kmoh  = 


[m2+]  [OH-] 
[MOH+] 


(7.3) 


which  can  be  compared  to  the  intrinsic  sorption  constants  to  hydrous  iron  oxide  [Dzom- 
bak  and  Morel,  1990] ,  and  the  sorption  constants  to  sediment  organic  carbon  [Mahony 
et  al.,  1996].  The  relationships  are  presented  in  Fig.  7.lA.Manganese  (II)  has  the 
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Fig.  7.1  (A)  Linear  free  energy  relationship  between  the  first  hydrolysis  constants  of  metal 

ions,  the  intrinsic  constants  of  sorption  to  hydrous  iron  oxide,  and  the  partition  coefficient  to 
organic  carbon  at  pH  =  7.  (B)  Partition  coefficients  versus  overlying  water  pH  for  the  MERL 
data  set. 


smallest  hydrolysis  constant  of  the  cations  listed  and,  therefore,  it  is  expected  to  have 
a  comparatively  small  partition  coefficient  to  hydrous  iron  oxide  and  particulate  or- 
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Table  7.1  Manganese  and  Cadmium  Partition  Coefficients 


pH 

Observed  ttmu 

Estimated  TMn,oe 

Measured  Ttod  nr. 

(L/kg) 

(L/kgOC) 

(L/kgOC) 

7 

lO'^ 

10®-10® 

105.40 

8 

10® 

10®-10^ 

105.96 

ganic  carbon.  Iron  (II)  is  in  the  middle  range  so  it  would  be  expected  to  have  sorption 
characteristics  that  are  similar  to  and 

Fig.  7.  IB  presents  the  observed  partition  coefficient  for  manganese  7rji^„  from  the 
MERL  experiment,  computed  from  the  overlying  water  particulate  [Mnp]  and  dis¬ 
solved  [Mrid]  manganese  concentrations  and  the  concentration  of  suspended  solids 
m 


'^Mn  — 


[Mrip]  /m 

[Mud] 


(7.4) 


This  computation  assumes  that  all  the  particulate  manganese  is  the  result  of  sorption 
onto  suspended  solids  and  not  a  result  of  the  oxidation  of  Mn  (II).  It  can  be  thought  of 
as  the  partition  coefficient  that  is  necessary  in  order  to  explain  the  observed  magnitude 
of  particulate  manganese.  The  results  for  pH  =  7  and  8  are  show  in  Table  7.1.  The 
magnitudes  of  -kmu  are  quite  large.  At  pH  =  7,  ttmu  =  10^  L/kg  (Fig.  7.1B).  If  the 
particles  were  between  1  and  10%  organic  carbon,  the  weight  fraction  of  organic  carbon 
would  be  in  the  range  /oc  =  0.01—0.1.  The  range  in  the  organic  carbon  normalized 
partition  coefficient 


T^Mn,OC 


Mn 

foe 


(7.5) 


would  be  10^-10®  L/kg  OC.  For  pH  =  8  the  range  would  be  t^muOC  =  10®-10^  L/kg 
OC. 

For  cadmium  the  measured  organic  carbon  normalized  partition  coefficients  are 
T^Cd,oc  —  10^'^  and  lO®'^®  L/kgOC  at  pH  =  7  and  8  respectively  [Mahony  et  al., 
1996].  Thus  'KMn,oc  would  need  to  be  at  least  as  large  as  Trcd,oc,  for  foe  =  0.1, 
and  a  factor  of  ten  larger  for  foe  =  0.01,  in  order  to  account  for  the  partition  co¬ 
efficients  shown  in  Fig.  7.1B.  But  manganese  has  a  hydrolysis  equilibrium  constant, 
10®-^'^,  which  is  approximately  one-half  order  of  magnitude  smaller  than  for  cadmium, 
103.92  jg  ujjikely  that  partitioning  alone  is  the  mechanism  that  generates 

particxilate  manganese.  Of  coiurse,  the  partition  coefficient  estimated  from  the  Hnear 
free  energy  relationship  could  be  used  in  the  model,  but  its  influence  would  be  small 
as  demonstrated  by  the  results  in  Table  7.1. 
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7.2  MODEL  FORMULATION 

The  model  formulation  is  shown  in  Fig.  7.2.  The  water  column  segment  is  linked  to 
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Fig.  7.2  Schematic  of  the  manganese  sediment  flux  model  including  the  overlying  water. 

the  aerobic  layer  by  the  transport  terms:  the  diflFusive  exchange  and  particle  settling. 
The  oxidation  of  Mn  {II)  is  included.  The  water  column  transport  is  represented  by 
an  advective  flow,  which  is  the  conventional  representation  for  a  completely  mixed 
lake  model,  for  example  [Thomann  and  Mueller,  1987;  Schnoor,  1996;  Chapra,  1997]. 
The  model  equations  that  result  for  the  water  column  concentrations  are: 

Layer  0  Mn{II): 


=  s  (/diMn(l)  -  Mn  (0))  -  KMn,oMn(0)  +  qo  {Mn  (0)^  -  Mn  (0))  (7.6) 
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Layer  0  MnOa  (s); 

dMn02{Q)  ,  . 

Ho - ^ - =  «Arn,oMn(0)  -  wiMn02  (0)  +  qo  (Mn02  (0%  -  Mn02  (0))  (7.7) 

where 

R'Mnfi  —  k^nfiHo  (7.8) 

the  usual  expression  for  a  reaction  velocity. 

The  terms  in  equations  (7.6-7.7)  represent  the  flux  from  the  sediment  to  the  over- 
lying  water 

J  [Mn]  =  s  (/diMn(l)  -  Mn  (0))  (7.9) 

the  oxidation  of  Mn  (II)  to  Mn02  (s) ,  the  inflowing  source  of  dissolved  manganese 
qoMn  (0)^  and  the  loss  by  outflow  qoMn  (0) .  The  flow  rate  qo  is  included  as  an  overflow 
rate  (m/d) 


where  Q  is  the  volumetric  flow  rate  (m^/d)  and  A  is  the  surface  area  of  the  water 
column  (m^),  so  that  the  units  of  qo  (m/ d)  are  consistent  with  the  other  mass  transport 
coefficients  in  the  equations. 

The  layer  1  Mn(l)  equation  (5.17a)  is  coupled  to  the  layer  0  equation  through 
Mn  (0)  in  eq.  (5.17a)  which  is  now  computed  using  eq.(5.17a).  In  addition  the  source 
term  in  the  layer  1  Mn02  (1)  eq.(5.17c)  becomes 

«4fn02  =  wiMn02  (0)  (7.11) 

which  couples  the  sediment  to  the  overlying  water  particulate  concentration  computed 
by  eq.(7.7).  Thus  the  sediment  manganese  mass  balance  equations  become 
Layer  1  Mn(II): 

dMn(V\ 

Hi — —  =  -s(fdiMn(l)-Mn(Q))  +  KLi2(fd2Mn(2)- fdiMn(l))  (7.12) 


^2 

-\-Wi2(fp2Mn{2)  -  /piMn(l))  -  W2Mn  (1)  -  -^^^f^^Mn(l) 
+Mn  (2)  ^1+  -  Mn  (1)  [Hx  +  Ff ) 

Layer  2  Mn(II): 

dMn{2)  ,  / , 

H2 — —  =  -KLi2(fd2Mn(2)  -  fdiMn(l))  (7.13) 

-wi2(fp2Mn(2)  -  fpxMn(l))  -  W2  (Mn{2)  -  Mn(l)) 
+«Mn,2Mn02(2)  -  Mn(2)  (iJs  +  Ht)  +  Mn(l)H^ 
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Layer  1  MnOi  (s): 

„rfMn02(l)  ,,  ,,,  ^ 

Hi - ^ -  =  — - — fdiMn{l) -W2Mn02{l) 

+wi2(Mn02{2)  —  Mn02(l))  +  wiMn02  (0) 

+Mn02(2)H+  -  Mn02{l)  (iJi  +  Ff ) 

Layer  2  Mn02  (s): 

„  dMn02(2) 

Hi - — -  =  -KMn,2Mn02(2) 


(7.14) 


(7.15) 


-wi2{Mn02{2)  -  Mn02{l))  -  W2  {Mn02{2)  -  Mn02(l)) 
-Mn02{2)  {h2  +  +  Mn02{l)Hi 

The  result  is  six  coupled  equations  in  the  six  unknown  concentrations  of  Mn{II) 
and  Mn02  (s)  in  the  three  layers. 


7.2.1  Analytical  Solutions 

The  steady  state  equations:  eqs. (7.6-7. 15)  with  dMn/dt  =  0  and  dMn02/dt  =  0,  can 
be  solved  analytically  for  the  general  case,  but  the  solutions  are  too  unwieldy  to  be 
useful.  For  the  special  case  where  the  loss  by  burial  in  layer  2  is  small  relative  to  the 
loss  by  outflow  in  layer  0  and  «;2  =  0  is  a  reasonable  approximation,  the  equations  can 
be  solved  and  reduced  to  simple  expressions. 

The  steady  state  equations  are 
Layer  0  Mn(II): 

0  =  s  {fdiMn{l)  -  Mn  (0))  -  /CMn,oMn(0)  +  qo  (Mn  (0)^  -  Mn  (0))  (7.16a) 

Layer  0  Mn02  (s): 

0  =  KMn,oMn(0)  -  wiMn02  (0)  +  qo  {Mn02  (0)^  -  Mn02  (0))  (7.16b) 

Layer  1  Mn(II): 

0  =  -s{fdiMn{l)-Mn{0))  +  KLi2(fd2Mn{2)-fdiMn{l))  (7.16c) 

+wi2(fp2Mn(2)  -  fpiMnil))  -  -^/^iMn(l) 

s 

Layer  1  Mn02  {s): 

0  =  fdiMn(l)  +  wi2{Mn02i2)  -  Mn02{l))  +  wiMn02  (0)  (7.16d) 
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Layer  2  Mn{II): 

0  =  -KLi2{fd2Mn{2)  -  fdiMn(l))  -  wi2ifp2Mn{2)  -  fpiMn{l))  (7.16e) 
Layer  2  Mn02  (s): 

0  =  -KMn,2Mn02{2)  -  wi2{Mn02{2)  -  Mn02{l))  (7.16f) 

These  six  simultaneous  equations  (7.16)  can  be  solved  analytically  to  yield  a  re¬ 
markably  simple  and  instructive  set  of  solutions.  The  layer  0  concentrations  are 


Mn(0)  =  Mn{0)i 


wi+qo 


+  90  +  «Mn,0 

MnO2(0)  =  Mn(0)  - 

Note  that 


-I-  Mn02  (0)j 


Wl 


Wi+qo  +  KMn,0 


(7.17) 


-+MnO2{0)i  (7.18) 

Wl  +  90  +  F,Mn,0  Wi+qo  +  KMnfl 


Mn  (0)  -f  MnOa  (0)  =  Mn  (0)^  -I-  Mn02  (0)^  (7.19) 

which  follows  from  a  mass  balance  aroimd  the  entire  water  column-sediment  system 

90  (Mn  (0)i  +  Mn02  (0)^)  =  90  (Mn  (0)  +  Mn02  (0))  (7.20) 


The  mass  flux  of  inflowing  manganese  must  be  equal  to  the  only  loss  terms  for  man¬ 
ganese,  namely  the  layer  0  outflows.  This  is  why  the  assumption  W2  =  0  simplifles  the 
situation  dramatically. 


7.2.1. 1  Limiting  Forms  The  concentrations  of  Mn{II)  and  Mn02  (s)  that  result  in 
the  water  column  (eqs.7.17-7.18)  are  controlled  by  ratios  involving  KMnfl,  Qo,  S'lid  wi 
which  are  all  parameters  for  water  column  processes:  oxidation,  outflow  and  settling 
to  the  sediment.  No  sediment  process  parameters  are  involved.  This  is  a  surprising 
result  which  requires  an  explanation. 

The  water  colunm  concentrations  are  linear  functions  of  the  inflowing  concentrations 
weighted  by  ratios  with  the  sum  of  KMnfl,  9o,  and  wi  in  the  denominator.  Consider 
flrst  what  happens  if  the  flow  rate  through  the  water  column  qo  is  large.  Then  the 
solutions  become 


Mn(0)  = 


Mn(0), 


Wl  -1-90 


Wl  +  90  4-  KMnfl 


-I-  Mn02  (0)j 


Wl 


Wl  +  9o  +  KMnfl 


Mn  (0)i  ^  -b  Mn02  (0)^  ^ 

as  90-^00  go  90 


(7.21) 


Mn(0), 
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Mn02  (0)  = 


as  go-^oo 


Mn  (0), - - +  Mn02  (0),  ■  + 

+  90  +  X’Mnfi  Wi+qo  +  KMnfi 

Mn  (0) .  +  Mn02  (0)^  ^ 

9o  go 


Mn02  {0\ 


(7.22) 

i.e.  there  is  no  change  in  the  water  column  concentrations;  they  remain  at  the  inflow 
values.  This  is  entirely  reasonable.  The  large  inflow  provides  an  infinite  supply  of 
Mn  {II)  and  Mn02  (s)  so  that  whatever  happens  in  the  sediment  has  no  effect. 
Consider  what  happens  if  the  settling  velocity  wi  is  large 


Mn(0) 


Mn(0), 


m+Qo 


Wi  +  go  +  KMnfi 


+  Mn02  (0)j 


toi 


^1  +  go  +  I^Mnfi 


^  Mn  (0).  ^  +  Mn02  (0).  ^ 

as  wi^oo  Wi  Wi 


(7.23) 


Mn  (0)j  +  Mn02  (0)^ 


MnO2(0)  = 


Mn(0)i 


K'Mnfi 


lOl  +  go  +  KMnfi 


+  Mn02  (0)j 


I^Mnfi  +  go 
^1  +  go  +  KMnfi 


as  wi—¥oo 


Mn  (0)i  +  Mn02  (0). 

wi  '  Wi 


0 

(7.24) 

This  is  a  more  interesting  result;  all  the  inflowing  manganese  appears  as  dissolved 
Mn{II). 

The  explanation  can  be  arrived  at  by  considering  the  constraint  imposed  by  mass 
balance.  Since  the  settling  velocity  is  becoming  large,  the  particulate  concentration 
in  the  water  column  is  approaching  zero.  Therefore  all  the  inflowing  manganese  must 
escape  via  the  dissolved  concentration  in  the  water  column  since — and  this  is  why  the 
simplification  102  =  0  is  so  powerful — there  is  no  other  loss  pathway.  So  the  dissolved 
concentration  increases  to  what  is  necessary  to  satisfy  mass  balance,  namely  the  total 
inflowing  concentration  of  manganese,  both  particulate  and  dissolved. 
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The  converse  situation  occurs  if  the  oxidation  rate  KMn,o  is  large 

Mn{0)  =  Mn{0),-  - +A/nO2(0), - ^ _ 

*^1  +  90  "I”  ^Mn,0  Wi  +  9o  +  ^Mn.,0 

^  Mn  (0)i  +  MnO^  (0),- 

as  «Mn,o-oo  KMnfi  ^  KMn.O 


(7.25) 


Mn02  (0)  = 


Mn{0\ 


Wi+qo  +  KMn,0 


+  Mn02  (0)i  + 

+  90  +  KMnfi 


Mn  (0)i  +  Mn02  (01 

->■  Mn  (0)^  +  Mn02  (0)^ 

(7.26) 

All  the  dissolved  manganese  in  the  water  column  oxidizes  to  particulate  manganese 
which  must  increase  in  concentration  to  the  total  inflowing  manganese  concentration 
so  that  the  outflowing  flux  equals  the  inflowing  flux. 

7 .2.1.2  Other  Layers  The  sediment  concentrations  of  Mn  (II)  and  Mn02  (s)  in  lay¬ 
ers  1  and  2  are  also  reasonably  concise  expressions 


Mn  (1)  =  KMnfiWi)  +  Mn02  (0)^  W\  (a  +  go  -f  KMnfi) 

sfdl  (wi  +  go  +  KMnfi) 


(7.27a) 


Mn02  (1)  =  Mn02  (2) 


ti'lZ  +  KMnfi 


(7.27b) 


Mn  (2)  =  A/n  (1)  +  MnO,  (2)  -  (7.27c) 

fd2KL\2  +  /p2^12  fd2KL12  +  fp2'^l2  ^ 

(Mn(0),  +  MnO, (0),)  ^ 


Mn02  (2) 


+9os  ^ — j^Mn  (0)^  +  wiMn02  (0), 
KMnfiS  (Wi  -I-  go  +  KMnfi) 


(7.27d) 
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The  importance  of  the  various  parameter  groupings  will  be  examined  using  the  nu¬ 
merical  results  presented  below  as  a  guide. 

7.2.1.3  Manganese  Flux  The  most  surprising  result  of  all  is  the  solution  for  the 
sediment  flux  of  Mn  (II) .  It  is  obtained  by  evaluating  the  deflning  equation 

J  [Mn]  =  s  (fdiMn(l)  -  Mn  (0))  (7.28) 

using  the  solutions  for  Mn  (1)  and  Mn  (0)  given  above  (eqs.7.17  and  7.27a)  and  using 
eq.(7.18)  to  substitute  for  the  resulting  expression.  The  result  is 

J  [Mn]  =  wiMnOi  (0)  (7.29) 

It  is  equal  to  the  flux  of  particulate  MnO^  (s)  that  settles  from  the  water  column  to 
the  sediment.  There  is  a  complete  recycling  of  the  particulate  flux  to  the  sediment  as 
a  dissolved  fliix  from  the  sediment. 

This  remarkable  prediction  appears  to  be  substantiated  by  the  MERL  experimental 
observations.  Fig.7.3  compares  the  manganese  flux  J[Mn]  to  the  particulate  man¬ 
ganese  concentration  in  the  overlying  water  Mn(0)p.  The  point  by  point  comparison 
(Fig.7.3A)  shows  no  apparent  relationship  although  the  linear  regression  line  included 
in  the  flgure  does  indicate  that  a  relationship  exists.  If  the  data  are  grouped  (Fig.7.3B) 
a  clear  relationship  emerges.  The  line  corresponds  to  a  settling  velocity  of  wi  —  1.0 
m/d. 


Fig.  7.3  (A)  The  observed  manganese  flux  J  [Mn]  versus  the  particulate  manganese  concen¬ 
tration  m  the  overlying  water  Mn(0)p.  Only  data  for  which  Mn(0)p  >  1/ig/L  are  included  in 
order  to  reduce  the  influence  of  analytical  errors  at  low  concentrations.  The  linear  regression 
line  IS  shown.  (B)  The  data  are  grouped  in  0.2  logio  units.  The  mean  ±  the  standard  deviation 
IS  shown.  The  line  corresponds  to  a  settling  velocity  of  wi  =  1  m/d. 

We  can  imderstand  the  flux  solution,  eq.(7.29),  by  examining  what  happens  if  no 
Mn02  (s)  forms  in  the  water  column.  This  happens  if  no  particulate  manganese  flows 
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into  the  water  column,  Mn02  (O)j  =  0,  and  no  manganese  dioxide  forms  by  oxidation 
of  Mn  (//) ,  so  that  KMn,o  =  0.  For  this  case  the  concentrations  in  the  overlying  water 
(eqs.7.17-7.18)  become 


MnO2(0)  =  0  (7.30) 

Mn  (0)  =  Mn  (0)^  (7.31) 


The  pore  water  concentration  in  the  aerobic  layer  equilibrates  with  the  overlying  water 
concentration  (eq.7.27c) 


(l)d  =  fdiMn  (1)  =  Mn  (0) 


(7.32) 


and  all  the  other  sediment  concentrations  adjust  accordingly  to  this  concentration. 
Since  there  are  no  mass  balance  consequences — ^these  concentrations  can  be  any  value 
whatever  without  a  need  to  remove  manganese  to  maintain  mass  balance — ^no  flux 
from  the  sediment  is  generated. 

Note  that  no  other  solution  is  possible.  If  the  pore  water  concentration  in  layer  1  was 
different  from  the  layer  0  concentration,  that  would  imply  that  a  flux  of  manganese  is 
present  either  from  or  to  the  overlying  water.  But  if  that  occurred,  then  the  sediment 
concentrations  would  either  be  continually  depleted  or  build  up  indefinitely,  neither 
of  which  is  a  steady  state  solution.  Therefore  the  zero  flux  case  is  the  only  allowable 
solution. 

If  Mn02  (0)  7^  0  so  that  there  is  a  flux  wiMn02  (0)  to  the  sediment,  then  a  flux 
from  the  sediment  is  necessary  for  mass  balance  to  be  maintained.  Since  the  only 
possible  flux  is  a  dissolved  flux  J  [Mn\  back  to  the  overlying  water — a  burial  flux  has 
been  eliminated  as  a  possibihty  by  assuming  W2  =  0— the  concentrations  adjust  to 
exactly  return  this  incoming  flux  (eq.7.29). 

This  remarkable  result:  that  the  sediment  flux  is  controlled  by  water  column  pro¬ 
cesses  and,  in  particular,  the  particulate  concentration  and  settling  velocity,  was  even¬ 
tually  noticed  after  many  numerical  simulations  were  performed.  It  provided  the  moti¬ 
vation  for  seeking  analytical  solutions  to  establish  the  precise  form  of  the  dependence. 
The  solutions  for  the  general  case,  u;2  ^  0,  are  too  complex  to  be  useful,  but  removing 
the  possibility  of  a  burial  flux  simplifies  the  solutions  so  that  the  relationship  between 
the  particulate  flux  to  the  sediment  and  the  recycled  dissolved  flux  emerges.  The 
utility  of  these  analytical  solutions  to  help  in  the  understanding  of  the  time  variable 
numerical  results  is  examined  next. 
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7.3  TIME  VARIABLE  MODEL 

The  complete  model  includes  the  overlying  water  and  the  sediment  mass  balance 
equations.  The  equations  are  discussed  first,  followed  by  the  a  sequence  of  numerical 
calculations  of  increasing  complexity  and  realism. 


7.3.1  Vector-Matrix  Equation 

The  equations  used  for  the  numerical  calculations  are  expressed  in  vector-matrix  form. 
They  comprise  the  two  layer  model  (eqs.  5.18-5.23)  with  the  addition  of  the  water 
column  as  the  third  layer.  The  manganese  concentration  state  variable  vector  is 

■Mn(l) 

Mn(2) 

Mn02(l) 

^  ’  Mn02  (2) 

Mn(0) 

Mn02  (0)  _ 

where  two  additional  rows  have  been  added  for  the  water  column  concentrations.  The 
matrix  A{t)  =  [aij]  has  elements 

an  =  (7.34a) 

-(i7i+^r) 


=  hi  +  H+ 


(7.34b) 


fll5 

=  s 

(7.34c) 

021 

=  fl2  +  102  + 

(7.34d) 

022 

=  -f21-W2-(^H2  +  H+"j 

(7.34e) 

024 

(7.34f) 

031 

=  0,(0)  A, 

(7.34g) 

O33 

=  -W12  -W2-  (^Hi  +  Hi'j 

(7.34h) 

O34 

=  W12  + 

(7.34i) 

036 

=r  Wl 

(7.34j) 
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<*43  =  W12  +W2+ Hi 

(7.34k) 

044  =  -KMn,2  “  W12  -  W2  -  (H2  +  H^^  + 

(7.341) 

®51  =  Sfdi 

(7.34m) 

055  =  -S  -  H’o/to  -  KMn,0  j  O2  (0) 

(7.34n) 

“65  =  KA/n.O  O2  (0) 

(7.34o) 

“66  =  —Ho/to  —  Wi 

(7.34d1 

where  /12  and  /2i  are  the  effective  mass  transfer  coefficients  from  layer  1  to  layer  2 
and  from  layer  2  to  layer  1  respectively  as  before  (eq.5.21a-5.21b)  and  aU  other  a^  ’s 
are  zero.  The  forcing  function  is  ^ 

0 

0 

=  0 

qoMn  (0)  ■ 

_  qoMn02  (0)- 

with  the  only  sources  being  the  Mn{II)  and  Mn02(s)  flowing  into  layer  1.  The 

diagonal  matrix  of  layer  depths,  H  (t) ,  now  includes  the  overlying  water  segment 
depth 

H(t)  =  diag[iyi(t),  H2(t),  H2(t),  Ho,  Ho]  (7.36) 

An  implicit  finite  diflference  equation  (eq.5.26)  is  used  to  solve  these  equations. 

7.3.2  Comparison  to  Steady  State  Results 

How  useful  are  the  steady  state  solutions  in  understanding  the  behavior  of  the  time 
variable  model?  The  first  step  is  to  check  the  results  of  the  time  variable  model  when 
It  IS  configured  to  be  a  steady  state  computation.  This  is  accomplished  by  disabling 
all  the  time  varying  elements.  Thus  all  the  parameters  that  cause  the  reaction  rates 
to  vary  with  temperature  (the  0’s)  are  set  equal  to  one.  This  includes  the  nutrient 
and  oxygen  model  parameters— for  example  the  diagenesis  reactions.  The  inflowing 
concentrations  Mn  (0)-  and  Mn02  (0)^  are  set  to  constant  average  values.  The  only 
time  variation  that  is  allowed  to  remain  is  that  induced  by  the  var5dng  overlying  water 
dissolved  oxygen  variation  O2  (0).  The  reason  is  to  ascertain  the  effect  of  just  this 
source  of  variation.  The  burial  velocity  is  set  to  zero  as  well  to  correspond  to  the 
steady  state  solutions. 
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The  time  variable  model  is  initialized  by  cycling  the  inputs  for  the  first  year  un¬ 
til  a  periodic  steady  state  is  achieved.  Then  average  concentrations  are  computed. 
In  addition,  annual  average  parameters  and  forcing  functions  are  computed  and  the 
steady  state  equations  are  evaluated  using  these  values,  which  are  presented  in  Table 
7.2.  This  procedure  is  followed  for  all  the  time  variable-steady  state  comparisons. 


Comparison  of  steady  state  and  tune  variable  solutions.  Annual  averages  for  the 
periodic  steady  state  initialization,  and  the  three  years  of  simulation.  (A)  Mass  balances: 
flux  in  (Jin)  and  flux  out  (J<,„t).  (B)  Overlying  water  concentrations.  (C)  Aerobic  layer.  (D) 
Anaerobic  layer.  The  subscript  ”ss”  indicates  the  steady  state  analytical  solution. 


The  results  are  compared  to  the  annual  average  fluxes  and  concentrations  in  Fig.7.4 
as  year  1.  The  subsequent  years  2-4  represent  the  results  of  the  simulation  of  the  three 
year  MERL  experiment.  The  mass  balances  are  checked  in  Fig.7.4A.  Since  the  year  one 
results  are  for  an  equilibrated  model,  the  inflowing  flux  balances  the  out-flowing  flux 
perfectly.  However,  the  year  2  and  3  time  variable  model  results  show  an  imbalance 
with  the  outflow  being  larger  than  the  inflow.  This  is  occurring  because  the  sediment 
is  releasing  more  manganese  than  is  settling  to  the  sediment  during  this  period.  The 
cause  is  the  change  in  overlying  water  conditions  at  the  onset  of  the  experiment.  Of 
course  the  steady  state  solutions  are  perfectly  in  balance  since  they  are  the  result  that 
would  be  expected  after  the  system  has  adjusted  to  the  changed  conditions. 

The  steady  state  and  annual  average  concentrations  for  the  overlying  water  (B), 
aerobic  layer  (C),  and  anaerobic  layer  (D)  are  somewhat  different,  reflecting  the  lack 
of  steady  state  conditions  in  the  time  variable  results.  However  the  magnitudes  are 
correct,  which  suggests  that  the  steady  state  solutions  can  be  used  to  guide  calibration 
and  to  explain  the  time  variable  results. 
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Table  7.2  Steady  State  Sediment  Model  Parameters 


Treatment 

pH 

T 

Ho 

to 

s 

^Mn,0 

90 

Wi 

Mn(0). 
Mn02  (0)i 

Mr^ss 

Mn  (0) 

MnO2(0),, 

Mn02  (0) 

fdl 

fpi 

Wi2 

Kl12 

J[Mn]ss 

J[Mn\ 

Mn(l) 
^”^2  (1).,., 
Mn02  (1) 
H2 

W2 

fd2 

fp2 

^Mn,2 

Mn(2),, 

Mn(2) 

Mn02  (2) 


0.000 

8.093 

10.644 

5.000“ 

27.400* 

0.089“ 

0.963“ 

0.182“ 
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9.332^^ 
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11.232'^ 

8.863'* 

5.581'* 

7.603'* 

0.000951 

0.999 

0.048“ 

0.500“ 

1.42x10"'^“ 
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0.“ 
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98.342** 

11.161** 
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2.000 
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27.400 
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0.182 

2.000 

9.332 
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10.927 

9.878 

5.886 

8.359 

0.000951 
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0.069 

0.500 

1.58x10-^ 
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8.625 

115.267 

91.045 

11.942 
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0. 

0.00285 
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114.090 
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0.171 
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11.371 
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0. 
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10.883 
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32.000 
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0.223 
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F/g-.  7.5  Comparison  of  annual  average  observations  and  time  variable  model  results.  (A,D,G) 
1981,  (B,E,H)  1982,  (C,F,I)  1983. 
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The  annual  average  results  from  the  time  variable  model  are  compared  to  the  ob¬ 
servations  in  Fig.7.5. 


Manganese  Flux 


Fig.  7.6  Manganese  fluxes  for  the  control,  2X,  8X,  and  32X  treatment. 


7.3.2.1  Effect  of  Dissolved  Oxygen  Variation  It  is  instructive  to  examine  the  time 
variable  results  in  more  detail  since  there  is  only  one  source  of  time  variation  in 
the  computation:  the  varjdng  overlying  water  dissolved  oxygen  concentration.  The 
variation  in  manganese  flux  from  the  sediment  J  [Mn\  is  shown  in  Fig.  7.6  for  the 
control  and  three  treatments.  The  variation,  which  is  not  insignificant,  is  due  in  part 
to  the  variation  in  s  =  SOD  fO^  (0)  and  the  variation  in  the  depth  of  the  aerobic  layer 
H\  that  results  from  the  variation  in  O2  (0) . 

Overlying  water  dissolved  oxygen  concentration  also  affects  the  rate  at  which  dis¬ 
solved  manganese  (II)  oxidizes  and  forms  particulate  manganese  dioxide  MnOa  (0). 
The  result  is  a  variation  in  the  particulate  manganese  concentration  shown  in  Fig.7.7. 
Since  the  flux  to  the  sediment  is  wiMn02  (0),  its  variation  also  induces  an  additional 
variation  in  J[Mn]  (Fig.7.6).  The  combination  of  the  varying  oxidation  rate  and  the 
varying  flux  from  the  sediment  produces  a  variation  in  the  overljdng  water  dissolved 
manganese  concentration  Mn  (0)  as  shown  in  Fig.7.7. 
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Fig.  7.7  Variation  in  overlying  water  dissolved  and  particulate  manganese.  The  definitions  of 
the  solid  and  dashed  lines  are  the  same  as  in  Fig.965.0008  and  Fig.965.0009  to  be  discussed 
subsequently. 


7.4  CALIBRATION 

The  calibration  of  the  three  layer  model  to  the  MERL  data  set  is  presented  in  this 
section.  The  particulate  and  dissolved  concentrations  MnO^  (0)^  and  Mn  (0)^  in  the 
inflow  to  the  water  column  for  the  three  years  are  obtained  from  the  average  particulate 
and  dissolved  concentration  in  the  control  tank’s  water  columns  (Fig.7.8).  The  idea  is 
that  the  control  tanks  are  replicating  the  bay  and  therefore  the  concentrations  in  the 
control  tanks  are  similar  to  what  is  in  the  inflow  to  the  other  tanks.  For  the  latter 
part  of  the  third  year  the  actual  inflowing  concentrations  were  determined  and  these 
are  included  in  Fig.7.8.  They  are  similar  to  the  control  tank  concentrations. 

The  pff  in  the  three  layers  of  the  model  is  shown  in  Fig.7.9.  The  method  used  to 
compute  the  aerobic  layer  pifi  is  presented  in  Chapter  6.  The  rest  of  the  inputs  have 
been  described  in  Chapter  5. 

7.4.1  Model  Results 

The  model  is  equihbrated  to  a  periodic  steady  state  by  using  the  average  water  column 
concentrations  of  the  control  tanks  for  the  inflowing  concentrations  and  cycling  the 
model.  The  same  input  is  used  for  each  cycle  with  the  ending  condition  for  the 
previous  year  is  used  as  the  initial  conditions.  The  idea  is  to  simulate  the  state  of  the 
sediments  before  the  nutrient  enrichment  experiment  commenced.  The  presumption 
is  that  the  behavior  of  the  first  year  of  the  controls  is  characteristic  of  the  average 
situation  in  Narragansett  Bay  in  general.  Once  the  computation  is  initialized,  the 
treatment  loadings  are  initiated  and  the  model  is  run  for  the  three  years. 

The  annual  average  data  emd  model  results  for  manganese  fluxes,  water  column 
particulate  and  dissolved  concentrations,  are  shown  in  Fig.7.10.  The  model  parameters 
are  Hsted  in  Table  5.1. 
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Particulate  Mn 


Dissolved  Mn 


Fig.  7.8  Concentration  of  particulate  and  dissolved  manganese  in  two  control  mesocosms. 
The  lines  are  fitted  to  both  data  sets  and  they  are  used  as  the  inflowing  concentrations. 
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Fig.  710  Comparison  of  average  annual  data  and  model  results:  1981  (A,D,G);  1982  (B  E  H)- 

1983  (C,F,I).  Manganese  flux  (A,B,C).  Water  coliunn  particulate  (D,E,F)  and  dissolved  (G  H  li 
manganese.  v  j  »  y 
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The  manganese  fluxes  (A-C)  increase  with  increased  nutrient  loading  from  the  con¬ 
trol  (OX)  to  the  most  heavily  loaded  (32X).  The  magnitude  and  pattern  are  reproduced 
for  the  first  and  third  year.  The  large  flux  increases  in  the  IX  to  4X  treatments  for 
the  second  year  (B)  are  not  reproduced.  These  correspond  to  the  large  increase  in 
particulate  concentration  (E). 

The  time  variable  computed  and  observed  concentrations  in  the  overlying  water  are 
shown  in  Figs.7.11  and  7.12. 

Two  sets  of  model  computations  are  presented.  The  dashed  lines  are  the  model 
computed  dissolved  Mn  (0)  and  particulate  Mn02  (0)  concentrations.  These  results 
correspond  to  the  assumption  that  all  the  particulate  manganese  is  filterable  and  would 
be  quantified  as  particulate  manganese.  The  solid  lines  are  the  concentrations  that 
are  computed  if  a  fraction  of  the  particulate  manganese,  /colloidal,  is  not  filterable.  The 
reason  for  this  modification  is  discussed  below  in  the  next  section. 

The  model  correctly  computes  the  seasonal  variation  of  dissolved  manganese  and 
the  pattern  of  increasing  concentrations  to  the  8X  treatment  and  decreasing  thereafter 
(Figs.7.11).  The  differences  between  the  dashed  and  solid  lines  are  not  large  and  the 
data  do  not  clearly  favor  one  over  the  other.  There  is  a  systematic  under-prediction 
of  the  dissolved  manganese  during  the  winter  and  spring  of  each  year.  We  have  no 
explanation  for  this  discrepancy.  The  increase  in  pHo  that  occurs  (Fig.7.9)  causes 
a  rapid  increase  in  the  rate  at  which  the  dissolved  Mn{0)  oxidizes  to  particulate 
Mn02  (0) .  The  particulate  manganese  settles  from  the  water  column  to  the  sediment 
where  it  remains  for  a  period  of  time  and  then  is  recycled  to  the  overlying  water  as  a 
flux. 

The  manganese  fluxes  which  result  are  shown  in  Fig.7.13.  The  model  reproduces 
the  observed  seasonal  variation  as  well  as  the  increase  in  flux  as  the  tanks  are  enriched 
from  IX  to  32X.  The  flux  to  the  water  column  causes  the  dissolved  concentration  to 
increase  which  produces  the  seasonal  pattern. 

As  a  result  of  the  formation  of  particulate  manganese  in  the  water  column,  its  set¬ 
tling  to  the  sediment,  and  its  release  to  the  water  column,  there  is  a  seasonal  variation 
in  the  total  manganese  concentration  in  the  water  column  as  shown  in  Fig.7.14.  It 
might  be  supposed  that  the  variation  in  total  manganese  concentration  is  also  influ¬ 
enced  by  the  variation  in  input  concentrations,  Fig.7.8.  However,  this  is  not  the  case. 
It  is  possible  to  suppress  the  seasonal  variation  in  the  inflowing  water  by  using  con¬ 
stant  input  concentrations,  which  are  the  yearly  average  of  the  concentrations  used 
previously  (Fig.7.8),  and  are  shown  in  Fig.7.15.The  resulting  total  manganese  concen¬ 
trations  are  shown  in  Fig.7.16.  The  concentrations  are  almost  indistinguishable  from 
the  results  using  the  time  variable  inputs  (Fig.7.14).  This  is  also  true  for  all  the  other 
model  outputs.  Thus  the  temporal  variation  in  total  manganese  in  the  water  column 
is  generated  by  the  migration  of  manganese  from  the  sediment  to  the  water  column 
and  back  again  during  the  annual  cycle. 


Fig.  7.11  Computed  dissolved  Mn(0)  (dashed  line)  and  non-filterable  Mn(0) 
fcoiioidaiMnOa  (0)  (solid  line)  manganese  in  the  water  column  compared  to  observations. 
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Overlying  Water  Concentrations  -  Particulate  Mn 


I  1 98 1 - 1 - 1 982 - 1 - 1 983 - !  I - 1 98 1 - 1 - 1 982 - 1 - 1 983 - 1 
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1  1 98 1  I  1 982 - 1 - 1 983 - 1  I - 1 98 1 - 1 - 1982 - ] - 1 983 - 1 


Fig.  7.12  Computed  particulate  MnO^iO)  (clashed  line)  and  filterable  MnOiiQ) 
fcoiioidaiMn02  (0)  (solid  line)  manganese  in  the  water  column  compared  to  observations. 
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S  concentration  in  the  water  column.  Comparison  of  observations 
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Inflowing  Concentrations 
Particulate  Mn 


Dissolved  Mn 


Fig.  7.15  Inflowing  concentrations  for  the  constant  concentration  case. 
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Overlying  Water  Concentrations  -  Total  Mn 


I - 1 981 - 1 - 1 982 - 1 - 1 983 - 1 


Tot 3.1  ni3iig3iiGSG  concontrs-tion  in  th.0  W3,t6r  column.  Compsxison  of  observations 
and  model  output.  Constant  inflow  concentrations. 
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7.4.2  Colloidal  Fraction 

The  suggestion  that  a  fraction  of  the  particulate  manganese  is  non-filterable  and  is 
quantified  as  dissolved  manganese  is  not  based  on  the  model’s  ability  to  reproduce 
the  observed  particulate  and  dissolved  manganese  concentrations.  In  fact,  the  dashed 
lines  (no  colloidal  manganese)  capture  the  magnitude  of  the  peak  concentrations  better 
than  the  solid  fines  (Figs.7.12)  and  the  use  of  a  colloidal  fraction  does  not  appear  to 
significantly  change  the  fit  to  the  dissolved  data  (Fig.7.11).  However  if  the  data  are 
analyzed  in  another  way  then  a  problem  appears. 

The  fraction  of  manganese  in  the  water  column  that  is  particulate  is  given  by 

,  _  MnO2(0) 

~  Mn02  (0)  +  Mn  (0) 

The  variation  in  ^  as  a  function  of  the  pH  is  shown  in  Fig.  7. 17.  The  upper  solid 


Particulate  Fraction 


Fig.  7.17  P^iculate  fr^tion  versus  pH  in  the  water  column.  Model  results  assnminpr 
f colloidal  =  0  (upper  lines)  and  fcoiioidai  =  1/2  (lower  lines).  The  dashed  lines  are  the  ranee 
m  model  results  for  that  pH  interval.  ^ 

and  dashed  fines  are  the  model  average  and  range  over  each  pH  interval,  assuming 
that  all  the  particulate  manganese  Mn02  (0)  is  filterable  and  would  be  measured  as 
particulate  manganese.  The  model  kinetics  predict  that  all  the  manganese  should 
be  particulate  /p  ~  1  above  pH  =  8.5.  However  the  observations  rarely  exceed  a 
particulate  fraction  of  fp  =  0.5.  Since  the  pH  dependency  of  the  oxidation  kinetics  is 
experimentally  observed,  another  explanation  for  this  discrepancy  is  required. 

The  hypothesis  is  that  the  particulate  manganese  is  composed  of  a  colloidal  frac¬ 
tion,  which  cannot  be  separated  by  filtration.  It  is  known  that  freshly  precipitated 
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manganese  dioxide  is  largely  in  the  colloidal  particle  size  range  so  that  this  assump¬ 
tion  is  not  without  support  [Pankow,  1979].  To  be  consistent,  the  colloidal  fraction  is 
assumed  not  to  settle  to  the  sediment  so  that  wi  is  replaced  by  (1  —  /colloidal)  in 
eqs.(7.6-7.15). 

The  lower  solid  and  dashed  lines  in  Fig.7.17  are  the  model  results  assuming  that 
50%  of  the  particulate  manganese  is  colloidal  (/colloidal  =  0.5)  and  therefore  would  not 
be  filterable.  The  measured  dissolved  Mn^  (0)  and  particulate  Mup  (0)  manganese 
concentrations— the  dashed  lines  in  Figs.7.11-7.12— are  computed  from  the  truly  par- 
ticulate  and  dissolved  concentration  using 

MudiO)  =  Mn(0)-f-/coiioidaiMnO2(0)  (7.38) 

Mnp{0)  =  Mn02  (0)  - /coiioidaiMn02  (0)  (7.39) 

The  supposition  that  a  colloidal  fraction  is  erroneously  measured  as  a  dissolved 
specie  is  commonly  used  as  an  explanation  for  the  observed  decrease  in  partition 
coefficient  as  the  concentration  of  particles  increase  [Gschwend  and  Wu,  1985]  although 
an  alternate  explanation  has  been  offered  [Di  Toro,  1985].  That  effect  is  not  much 
different  from  the  observations  in  Fig.7.17  where  the  particulate  fraction  does  not 
exceed  fp  =  0.5,  corresponding  to  a  lower  partition  coefficient,  whereas  the  expectation 
is  that  /p  — >•  1  which  corresponds  to  a  high  partition  coefficient. 

7.4.3  Effect  of  Anaerobic  Layer  Depth 

The  depth  of  the  anaerobic  layer  H2  is  an  important  parameter  in  the  sediment  model 
since  it  controls  the  time  constant — ^the  characteristic  time  for  the  model  to  equili¬ 
brate  to  new  conditions — for  conservative  species  such  as  phosphorus  [Di  Toro  and 
Fitzpatrick,  1993].  It  is  true  for  the  manganese  model  as  well.  The  reason  is  that  the 
anaerobic  layer  provides  a  large  storage  volume  for  manganese  and  achieving  steady 
state  requires  that  the  concentration  Mn  (2)  in  H2  be  equilibrated.  This  can  hap¬ 
pen  in  two  ways.  Either  the  flux  to  or  from  the  sediment  J  [Mn]  ,or  the  burial  flux 
W2  (Mn  (2)  -t-  Mn02  (2)) ,  provides  the  route  for  which  Mn  (2)  adjusts.  The  time  vari¬ 
able  model  results  (Fig.7.4)  presented  in  section  7.3.2  had  H2  =  1  cm.  As  a  result, 
steady  state  was  achieved  relatively  rapidly,  even  without  a  loss  by  burial  since  the 
comparisons  to  the  analytical  model  steady  state  solutions  are  made  with  W2  =  0. 
The  question  remains:  what  is  the  appropriate  H2  for  the  MERL  nutrient  exchange 
experiment. 

The  manganese  concentration  in  the  anaerobic  layer  near  the  end  of  year  three  of 
the  simulation  is  compared  to  observations  in  Fig.7.18.  These  simulation  are  per¬ 
formed  with  H2  =  10  cm.  There  is  a  notable  decline  in  observed  sediment  manganese 
concentration  at  the  higher  enrichment  levels  which  the  model  does  not  capture.  The 
same  results  for  H2  =  2  cm.  are  shown  in  Fig.7.19.The  decline  in  Mn  (2)  with  in- 
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f'S-  ^18  Sediment  and  pore  water  manganese  concentrations  at  the  end  of  year  3  Solid  line 
IS  the  model  computation.  H2  =  10  cm.  ' 


Fig  7.19  Sediment  and  pore  water  manganese  concentrations  at  the  end  of  year  3.  Sohd  line 
IS  the  model  computation.  H2  =  2  cm. 


creasing  nutrient  inputs  is  the  result  of  increasing  fluxes  of  Mn  to  the  overlying  water 
in  the  first  year  of  the  simulation  (Fig.7.5).  Since  the  anaerobic  layer  is  thinner  than 
in  Fig.7.18,  the  sediment  concentrations  respond  more  rapidly.  From  these  results  it 
appears  that  the  thiimer  layer  H2  =  2  cm.  is  more  representative  of  the  short  term 
response  observed  in  these  experiments. 

It  should  be  pointed  out  that  the  depth  of  the  anaerobic  layer  is  an  artificial  con¬ 
struct  that  is  the  consequence  of  the  assumed  two  layer  structure  of  the  model.  If  the 
model  had  been  constructed  using  multiple  layers  in  the  vertical,  then  the  parame¬ 
ters  that  control  the  long  term  response  are  the  mixing  coefficient  for  particles  and 
the  depth  of  particle  mixing.  Although  this  formulation  is  preferable,  it  can  add  an 
unacceptably  large  computational  burden  if  the  model  is  to  be  applied  as  a  part  of 
a  large  water  quality  model  with  many  individual  sediment  segments.  The  depth  of 
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the  active  layer  H2  is  a  surrogate  for  these  parameters  and  plays  a  similar  role  in  the 
model. 


8 

Application  of  Iron  Flux  Model  to  Onondaga 

Lake 

8.1  INTRODUCTION 

This  chapter  presents  the  application  of  a  iron  sediment  flux  model  to  Onondaga 
Lake.  The  model  is  constructed  in  exactly  the  same  way  as  the  manganese  sediment 
and  water  column  flux  model  described  in  the  previous  chapters.  Onondaga  Lake 
was  chosen  because  a  complete  eutrophication  water  column  data  set  is  available  and 
has  been  modeled,  and  because  iron  data  are  also  available.  The  hjrpolimnion  of  the 
lake  becomes  anaerobic  each  year,  releasing  significant  quantities  of  iron  to  the  water 
column.  Thus  Onondaga  Lake  provides  a  very  good  setting  to  test  the  validity  of  the 
iron  sediment  flux  model. 

8.1.1  Physical  Characteristics 

Onondaga  Lake  is  located  immediately  north  of  the  city  of  Syracuse,  in  Onondaga 
County,  in  Central  New  York  State.  The  lake  has  a  length  of  7.6  km.  and  a  width  of 
1.8  km.  The  outflow  from  the  lake  exits  through  a  single  outlet  at  the  northern  end 
and  enters  the  Seneca  River.  A  bathymetric  map  is  shown  in  Fig.  8.1  and  relevant 
physical  data  is  given  in  Table  8.1.  The  is  lake  comprised  of  two  basins,  commonly 
referred  to  as  the  south  and  north  basins,  which  are  separated  by  a  slight  saddle  region 
that  is  located  approximately  3.6  km  from  the  outlet.  The  north  basin  has  a  maximum 
depth  of  18.9  m  and  the  south  basin  has  a  maximum  depth  of  19.9  m. 

Table  8.1  Physical  characteristics  of  Onondaga  Lake 


Characteristic 

Value 

Surface  Area  (km^) 

12.0 

Volmne  (m^) 

131.  X  10® 

Mean  Depth  (m) 

10.9 

Maximum  Depth  (m) 

19.9 

Watershed  Drained  (km^) 

689. 

Detention  Time  (days) 

90.0 
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Bathymetric  Mop  of 


Fig.  8.1  Bath}Tiietric  map  of  Onondaga  Lake 


8.1.2  Tributaries 

Onondaga  Lake  receives  surface  runoff  from  a  drainage  basin  of  689  km^  which  is 
almost  entirely  located  in  Onondaga  County.  The  major  freshwater  inputs  are  Nine- 
mile  Creek,  Onondaga  Creek,  the  Metropolitan  Syracuse  Sewage  Treatment  Plauit 
(METRO),  Ley  Creek,  Bloody  Brook  and  Harbor  Brook.  Minor  inputs  include  Sawmill 
Creek,  Tributary  5A  and  the  East  Flume.  Flow  data  for  the  tributaries  for  the  period 
1971-1989  is  given  in  Table  8.2. 

8.1.3  Loadings 

As  a  result  of  urban  and  industrial  development,  Onondaga  Lake  has  been  used  for 
both  domestic  and  industrial  waste  disposal  for  over  100  years.  Major  industrial  load¬ 
ings  include  calcium,  alkalinity,  inorganic  carbon  and  chlorides  that  were  associated 
with  the  production  of  sodium  carbonate  (Na2C03),  commonly  known  as  soda  ash. 
The  production  of  soda  ash  began  on  the  shores  of  Onondaga  Lake  in  1884  and  con¬ 
tinued  until  1986. 
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Table  8.2  Ttibutary  flows  into  Onondaga  Lake  1971-1989 


TVibutary 

Drainage  Area 
km^ 

Average  Flow 
(m^/s) 

%  Contribution 

Onondaga  Creek 

285 

5.25 

31.2 

Ninemile  Creek 

298 

5.29 

31.5 

METRO 

2.99 

17.8 

Ley  Creek 

77 

1.32 

7.78 

Harbor  Brook 

29 

0.38 

2.30 

Other  inputs 

1.57 

9.40 

Total 

689 

16.80 

100.0 

The  Metro  wastewater  treatment  plant  was  upgraded  in  1979  to  secondary  treat¬ 
ment.  Tertiary  treatment  for  phosphorus  removal  was  added  in  1981.  By  the  early 
80’s  the  plant  capacity  was  80  MGD.  Loading  data  are  hsted  in  Table  8.3 


Table  8.3  Summary  of  loadings  to  Onondaga  Lake  1971-1989 


Constituent 

1970’s 

early  1980’s 

late  1980’s 

Organic  Carbon 

45,000 

18,000 

7,000 

Total  nitrogen 

8,000 

6,000 

total  phosphorus 

800 

400 

Silica 

400 

250 

Calcium 

1  X  10® 

0.6  X  10® 

0.2  X  10® 

Alkalinity 

0.25  X  10® 

0.25  X  10® 

0.25  X  10® 

Total  Inorganic  Carbon 

4.0  X  10^ 

Chloride 

1.8  X  10® 

1.8  X  10® 

0.4  X  10® 

8.1.4  Water  Quality 

The  water  quality  in  Onondaga  Lake  during  the  period  simulated  below  is  characteris¬ 
tic  of  a  highly  eutrophic  lake.  Due  to  the  large  nutrient  loadings,  summer  chlorophylla 
concentrations  can  exceed  100  /xg/L.  The  lake  is  stratified  in  the  summer  months 
from  June  to  September  with  temperature  differences  exceeding  18  °C.  Hypolimnitic 
hypoxia  and  anoxia  exist  for  9  to  10  months  per  year.  Epilimnion  dissolved  oxygen 
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concentrations  are  generally  above  6  mg/L  throughout  the  year.  Table  8.4  provides 
typical  siunmer  and  winter  concentrations. 

Table  8.4  Onondaga  Lake  water  quality  overview. 


Variable 

Summer 

Winter 

Epilimnion  Temp  (°C) 

22 

3 

Hjqjolinmion  Temp  (°C) 

7 

3 

Epilimnion  DO  (mg/L) 

13 

8 

Hypolimnion  DO  (mg/L) 

0 

1 

Chlorophyll_a  ( fig/L) 

40 

8 

Epilimnion  Salinity  (g/L) 

1.5 

1.5 

Hypolimnion  Salinity  (g/L) 

2.5 

1.7 

Epilimnion  pH 

7.8 

7.8 

Hypolimnion  pH 

7.3 

7.3 

8.1.5  Data  Sets 

A  substantial  data  set  is  available  for  Onondaga  Lake  dating  back  to  the  late  1960’s. 
A  compilation  and  review  of  all  data  collected  between  1968-1989  is  available  [Walker, 
1991].  Data  sources  include  Onondaga  County  Department  of  Drainage  and  Sanitation 
(D&S),  Onondaga  County  Department  of  Health  (DOH),  New  York  State  Department 
of  Environmental  Conservation  (DEC),  Upstate  Freshwater  Institute  (UFI)  and  the 
United  States  Geological  Survey  (USGS).  The  data  have  been  assembled  onto  a  single 
database  containing  186,000  water  quality  observations. 

8.1.5.1  Iron  and  Manganese  Dissolved  iron  data  were  collected  by  the  Department 
of  Drainage  and  Sanitation  in  the  years  1968-1975.  Dissolved  iron  data  were  also 
collected  by  the  Upstate  Freshwater  Institute  in  1980-1981, 1985-1986  and  1988-1989. 
Data  were  collected  at  station  41,  in  the  southern  basin  of  the  lake.  It  is  the  data 
from  this  station  that  will  be  used  subsequently.  Data  were  collected  at  1  m  depth 
intervals  on  a  weekly  basis.  In  1985  iron  data  was  recorded  from  May  to  November. 
In  1986  iron  data  for  March  through  November  are  available.  In  1987  no  iron  data 
were  recorded.  In  1988  data  were  recorded  from  May  to  October.  In  1989  iron  data 
are  available  from  June  to  August. 

Dissolved  manganese  data  were  collected  by  the  Department  of  Drainage  and  Sani¬ 
tation  in  the  years  1968-1970  which  is  the  only  manganese  data  available.  Because  of 
the  age  of  these  data  and  questions  over  their  validity,  manganese  was  not  considered. 
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8.2  MODEL  CONFIGURATION 

The  eutrophication  model  has  been  constructed  for  Onondaga  Lake  as  part  of  another 
project.  A  detailed  description  of  RCA,  the  water  quality  model  employed,  is  available 
[HydroQual,  1992].  The  model  is  direct  descendent  of  the  WASP  model  [Di  Toro 
et  al.,  1981].  For  Onondaga  Lake  the  version  of  RCA  used  had  32  state  variables, 
including  the  two  iron  state  variables  Fe(II)  and  FeOOH{s).  The  model  is  configured 
to  simulate  the  annual  cycle  of  phytoplankton  production,  nutrient  cycling,  and  the 
resulting  dissolved  oxygen  distribution.  A  sediment  flux  model  which  provides  the 
fluxes  of  nutrients  and  oxygen,  and  is  identical  to  that  used  for  the  MERL  simulations, 
is  coupled  to  the  water  column  [Di  Toro  and  Fitzpatrick,  1993].  A  hst  of  the  systems 
used  in  the  Onondaga  Lake  model  are  given  in  Table  8.5. 

8.2.1  Iron  Flux  Model  Equations 
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Fig.  8.2  Schematic  diagram  of  the  iron  flux  model 
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Table  8.5  Water  quality  model  systems 


Number 


System 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


Salinity 

Phytoplankton  -  diatoms 
Phytoplankton  -  summer  groups 
Particulate  Organic  Phosphorus  -  refractory 
Particulate  Organic  Phosphorus  -  labile 
Dissolved  Organic  Phosphorus  -  refractory 
Dissolved  Organic  Phosphorus  -  labile 
Total  Dissolved  Organic  Phosphorus 
Particulate  Organic  Nitrogen  -  refractory 
Particulate  Organic  Nitrogen  -  labile 
Dissolved  Organic  Nitrogen  -  refractory 
Dissolved  Organic  Nitrogen  -  labile 
Total  Ammonia 
Nitrate  +  Nitrate 
Biogenic  Silica 
Total  Silica 

Particulate  Organic  Carbon  -  refractory 
Particulate  Organic  Carbon  -  labile 
Dissolved  Organic  Carbon  -  refractory 
Dissolved  Organic  Carbon  -  labile 
Dissolved  Organic  Carbon  -  reactive 
Dissolved  Organic  Carbon  -  algal  exudate 
Oxygen  Equivalents 
Dissolved  Oxygen 
Total  Active  Metal 
Total  Inorganic  Carbon 
Alkalinity 
Calcitep 
Calcium 
Temperature 
Dissolved  Iron 
Particulate  Iron 
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The  formulation  of  the  sediment  flux  model  for  iron  exactly  parallels  that  developed 
for  manganese.  Iron  is  tracked  both  as  dissolved  and  particulate  FeOOH  (s) 
iron.  A  schematic  is  presented  in  Fig.8.2.  The  equations  for  the  aerobic  (1)  and 
anaerobic  (2)  sediment  layers  are: 

Layer  1  Fe{II): 


dFe{l) 


Layer  2  Fe{II): 

rr  dFe{2) 


-s  (/rfiFe(l)  -  Fe  (0))  +  KLn{fd2Fe{2)  -  /diFe(l)) 
+wi2(/p2Fe(2)  -  /piFe(l))  -  W2Fe  (1) 

-kpe,iHi  O2  (0)  fdiFeil) 

+Fe  (2)  -  Fe  (1)  (Hi  +  ) 


=  -Kn2{fd2Fe{2)-UiFe{l)) 

-wi2{fp2Fe(2)  -  /piFe(l))  -  W2  (Fe(2)  -  Fe(l)) 


(8.1a) 


Layer  1  FeOOH  (s): 
„  dFeOOH{l) 


+kFe,2H2FeOOH{2)  -  Fe(2)  (.^2  +  ^1+)  +  Fe(l)iff  (8.1b) 

=  kpciHi  O2  (0)  /diFe(l)  -  W2FeOOH{l) 

+wi2{FeOOH{2)  -  FeOOH{l))  +  wiFeOOH  (0) 
+FeOOH{2)H+  -  FeOOH{l)  (Hi  +  (8.1c) 


Layer  2  FeOOH  (s): 
„  dFeOOH{2) 


=  -kFe,2H2FeOOH{2) 

-wi2{FeOOH{2)  -  FeOOH{l))  -  W2  {FeOOH{2)  -  FeOOH{l)) 


-FeOOH{2)  {h2  +  Ht)  +  FeOOH(l)jyf  (8.1d) 

where  FeOOH  (0)  and  Fe  (0)  refer  to  the  particulate  and  dissolved  iron  concentrations 
in  the  bottom  most  water  column  segment,  designated  as  segment  0,  which  overlays 
the  sediment. 

The  water  colmnn  equations  for  iron  for  this  segment  are 
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Layer  0  Fe{II): 
dFe(O) 


^  s(/diFe(l)-Fe(0))-fcjre,o/^o(lo2(P"‘>-''))02(0)Fe(0) 


+Fo,_i  (Fe(-l)-Fe(O)) 


(8.2a) 


Layer  0  FeOOH  (s): 

„  dFeOOH{0) 
■no - T. - 


kpefiHo  ^io2(p^o-7)  j  (0)  Fe(0) 
+wi  {FeOOH  (-1)  -  FeOOH  (0)) 

KM,Fefi 


-kpe,lHo- 


FeOOH  (0) 


KM,Fefi  +  [O2  (0)] 

+E0-1  {FeOOH  (-1)  -  FeOOH  (0)) 


(8.2b) 


where  Fo -1  represents  the  vertical  mixing  coefficient  between  the  bottom  layer  and 
the  layer  immediately  above.  The  concentrations  in  the  water  column  are  numbered 
in  reverse  order  to  preserve  the  notation  used  previously.  Thus  layer  —1  refers  to  the 
layer  just  above  layer  0  in  the  water  column. 

An  additional  term  has  been  added  to  the  water  column  equations  (8.2)  to  include 
the  reduction  of  FeOOH  if  the  dissolved  oxygen  concentration  approaches  zero.  The 
rate  of  the  reduction  of  FeOOH  to  Fe{II)  is  given 


(8.3) 


where  kpe,\  is  the  rate  at  which  FeOOH{s)  is  reduced  to  Fe{II)\  KM,Fe,o  is  the  half 
saturation  constant  for  oxygen  for  this  reaction  which  is  set  to  a  small  number  so  that 
reduction  occurs  only  occurs  if  the  dissolved  oxygen  concentration  is  almost  zero.  The 
nitrogen,  oxygen,  phosphorus,  and  silica  portions  of  the  sediment  model  are  identical 
to  those  used  above. 


8.2.2  Segmentation  for  Onondaga  Lake 

Onondaga  Lake  is  represented  by  10  vertical  water  column  segments.  Each  segment 
has  a  depth  of  2  meters  and  a  surface  area  chosen  to  match  the  hypsometry  of  the 
lake.  Table  8.6  gives  the  depth,  surface  area  and  volume  of  each  vertical  segment.  The 
lake  is  assumed  to  be  laterally  well  mixed,  consistent  with  observations  [Effler,  1996]. 
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Table  8.6  Vertical  segment  depths,  areas  and  volumes 


Segment 

No. 

Average 

Depth 

m 

Siurface  Area 

m^  xlO® 

Volume 

m^  xlO® 

1 

1 

11.9 

23.8 

2 

3 

11.9 

22.3 

3 

5 

10.4 

19.9 

4 

7 

9.49 

18.3 

5 

9 

8.81 

16.9 

6 

11 

8.12 

15.0 

7 

13 

6.91 

13.0 

8 

15 

6.11 

10.9 

9 

17 

4.76 

8.61 

10 

19 

3.85 

3.85 

8.2.3  Iron  Loading 

Data  for  the  composition  of  particles  entering  Onondaga  Lake  in  1981  have  been 
reported  [Yin  and  Johnson,  1984].  The  loading  rate  in  1981  for  a  “Fe-Mn-rich”  class  of 
solids  is  reported  that  contains  Fe{OH)z  (s),  Fe^Oz  (s),  FeOOH  (s),  MnOz  (s),  and 
other  Fe  and  Mn  compounds.  These  are  shown  in  Table  8.7.  Since  the  concentration  of 
manganese  is  much  less  than  iron  in  typical  soil  particles  [Lindsay,  1979]  it  is  assumed 
that  these  loading  rates  represent  only  iron  inputs.  Using  the  average  tributary  flows 
an  estimate  of  the  influent  concentration  can  be  obtained.  This  concentration  was 
then  used  to  estimate  daily  loads  for  the  calibration  years  1985-1989. 

Table  8. 7  Iron  loading  to  Onondaga  Lake  for  1981 


Tributary 

Flow  1981 
m^/s 

Fe  Load 
tons/yr 

Onondaga  Creek 

5.1 

315 

Ninemile  Creek 

4.9 

63 

METRO 

3.1 

2 

Ley  Creek 

1.2 

22 

Harbor  Brook 

0.4 

2 

Total 

14.7 

405 
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This  loading  represents  all  particulate  iron  compounds.  In  order  to  obtain  the 
reactive  iron  load  it  is  necessary  to  reduce  the  total  load  to  accoimt  for  the  unreactive 
iron  compounds.  As  there  is  no  simple  way  to  do  this  without  sequential  extraction 
information  [Canfield,  1989],  the  load  was  reduced  as  part  of  the  calibration  procedure. 
A  reduction  of  40%  was  required  in  order  to  match  the  observations.  This  reduction 
was  applied  to  all  model  year  runs.  All  infiuent  iron  was  assumed  to  be  particulate 
since  inflows  to  the  lake  were  aerobic  and  the  oxidation  rate  of  ferrous  iron  is  very 
rapid  (Fig.2.4). 

8.2.4  Dissolved  Oxygen 

The  original  development  of  the  Onondaga  Lake  model  focused  on  calcitep  precipita¬ 
tion  and  deposition,  with  secondary  emphasis  on  eutrophication.  Although  the  model 
produced  a  reasonable  fit  of  the  dissolved  oxygen  profiles,  it  did  not  always  reproduce 
the  exact  timing  of  the  anoxia  in  Onondaga  Lake.  An  example  of  model-data  compar¬ 
ison  for  1986  is  shown  in  Figure  8.3.  The  model  is  unable  to  reproduce  the  observed 
hypolimnetic  anoxia  from  February  through  May.Because  of  the  strong  relationship 
between  iron  flux  from  the  sediment  and  the  overlying  water  dissolved  oxygen,  it  was 
decided  to  use  the  measured  dissolved  oxygen  data  rather  than  the  model  results  to 
specify  when  anoxia  was  present.  This  was  done  for  three  reasons. 

1.  The  dissolved  oxygen  data  were  very  complete  and  thorough,  with  weekly  dis¬ 
solved  oxygen  profiles  measured  at  1  m  depth  intervals. 

2.  Because  dissolved  oxygen  is  coupled  to  so  many  other  systems,  achieving  a  better 
dissolved  oxygen  fit  to  the  data  would  have  essentially  involved  recalibrating  the 
existing  water  quality  model,  which  would  have  been  time  consuming. 

3.  Because  the  iron  only  responds  to  the  dissolved  oxygen  but  does  not  affect  it, 
at  least  at  these  iron  concentrations  for  which  the  oxygen  consumed  by  iron 
oxidation  is  small,  there  is  no  loss  of  a  feedback  mechanism  by  specifying  the 
dissolved  oxygen  directly. 

8.2.5  Calibration  Results 

The  model  was  run  for  5  years,  beginning  January  1,  1985  to  December  31,  1989.  Iron 
loading  was  supplied  via  the  tributaries  based  on  the  calculation  procedure  described 
above.  The  values  of  the  constants  used  for  the  nutrients  and  oxygen  sediment  model, 
and  iron  model  in  the  water  column  and  sediment  are  given  below  in  Tables  8.9  and 
8.8.  Note  that  the  same  oxidation  and  reduction  rates  are  used  for  both  the  sediment 
and  water  column.  This  parallels  the  procedure  followed  in  the  manganese  model. 
The  model  was  cycled  through  20  years  of  simulations — ^the  5  year  sequence  repeated 
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LEGEND  Data:  T  T  Monthly  Avg  G  Mange  Model:  -  Surface 

I  X  Surf  G  Bot  Respectively  Bottom 

Fig.  8.3  Calibration  results  for  RCA  water  quality  model  for  1986.  Monthly  averages  for 
surface  layer:  model  (solid  line),  data  (a)  ±  range;  and  bottom  layer  model  (dotted  line),  data 
(t)  ±  range. 
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four  times — in  order  to  generate  equilibrium  initial  conditions.  These  were  then  used 
as  the  initial  conditions  for  the  final  five  year  calibration  run. 

Table  8.8  Constants  for  the  iron  flux  model  -  Onondaga  Lake 


Water  Column 

Name 

Value 

Unit 

Description 

10.0 

day“^ 

Oxidation  rate 

Gpe,! 

1.08 

- 

Temperature  coefficient 

fefe.2 

4.0 

day“^ 

Reduction  rate 

KM,Fefi 

0.2 

mg  O2/L  Oxidation  half  saturation  constant 

Sediment 


ATTiTe,! 

15.0 

L/kg 

Aerobic  layer  partition  coefficient 

T^Fe,2 

3000.0 

L/kg 

Anaerobic  layer  partition  coefficient 

O^^Crit.Fe 

2.0 

mg/L 

Partition  coefficient  half  saturation  constant 

^Fe.l 

10.0 

day“^ 

Aerobic  layer  oxidation  rate 

^Fe,l 

1.08 

- 

Temperature  coefficient 

kFe,2 

4.0 

day“^ 

Anaerobic  layer  reduction  rate 

^Fe,2 

1.08 

- 

Temperature  coefficient 

Figure  8.4  shows  the  iron  calibration  results  for  1986  using  the  fitted  dissolved  oxy¬ 
gen  for  the  months  of  March  to  November  when  iron  data  were  recorded.  Also  shown 
are  temperature  and  dissolved  oxygen.  The  measured  data  average  (open  symbols) 
and  range  (bars)  and  the  model  computed  monthly  averages  (solid  line)  and  range 
(shading)  are  compared.  It  is  interesting  to  note  how  early  in  the  year  anoxia  started 
and  how  long  it  persisted.  In  fact,  in  1986  the  bottom  of  the  lake  was  essentially 
anoxic  or  hypoxic  (DO  <  2  mg/L)  for  the  whole  year. 

The  dissolved  iron  concentrations,  shown  in  the  third  row  of  panels  in  Fig.  8.4, 
increase  as  a  consequence  of  the  iron  flux  being  generated  by  the  sediment.  The  data 
for  March  were  actually  collected  on  March  30  and  so  the  right  edge  of  the  shading 
should  be  compared  to  the  data,  to  which  it  compares  reasonably  well.  There  is  a 
continual  buildup  of  dissolved  iron  due  to  an  iron  flux  from  the  sediment. 

The  particulate  iron  concentrations  are  shown  in  the  last  row  of  Fig.8.4.  When 
anoxia  is  not  present  the  particulate  iron  will  exhibit  an  increasing  concentration  with 
depth,  due  to  the  settling  of  particulate  iron  into  the  deeper  parts  of  the  lake  where 
the  volume  is  decreasing.  As  the  particulate  iron  settles  into  the  anoxic  zone  it  is 
reduced  to  form  .  Thus  the  particulate  concentration  will  initially  increase  and 
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Table  8.9  Sediment  Model  Parameters  for  Onondaga  Lake 


mi 

(kg/L) 

0.5 

m2 

(kg/L) 

0.5 

W2 

(m/d) 

6.85  X  10“® 

H2 

(m) 

0.10 

Dd 

Dp 

Ddo 

^Ddo 

(mVd) 

- 

(mVd) 

- 

(mVd) 

- 

5.0E-03 

1.08 

1.2  X  lO-'^ 

1.15 

1.0  X  10"^ 

1.08 

KNHi,! 

TnH4 

^NHi 

Km,nh4 

Km,nh4,02 

(m/d) 

(L/kg) 

- 

(mg  N/L) 

- 

(mg  O2/L) 

0.2 

1.00 

1.123 

0.728 

1.125 

0.74 

I^NOsA 

K-N03,2 

^NOs 

(m/d) 

(m/d) 

- 

0.15 

0.25 

1.080 

K'H2S,dl 

«K25,pl 

T^H2S,2 

Km,H2S,02 

(m/d) 

(m/d) 

(L/kg) 

(L/kg) 

- 

(mg  O2/L) 

0.2 

0.4 

100.0 

100.0 

1.15 

4.0 

ksi 

[Si].at 

AtTSi,! 

7rsi,2 

^Si 

Km,p  Si 

i/d) 

(mg  Si/L) 

- 

(L/kg) 

- 

(mg  Si/L) 

0.5 

30. 

6.0 

175.0 

1.100 

5.0  X  lO'* 

A7rpo4,i 

TP04,2 

^2  (0),r,PO4 

^ M,Dp,02 

O2  (0),,.si 

- 

(L/kg) 

(mg  O2/L) 

(mg  O2/L) 

(mg  O2/L) 

65.0 

225.0 

2.0 

4.0 

1.0 

then  decrease  in  the  anoxic  zone.  The  model  computations  are  shown  in  Figure  8.4. 
Unfortunately  there  are  no  data  with  which  to  compare  the  model  results. 

From  March  to  May  the  Fe^+  concentration  at  the  hypolimnion  continues  to  in¬ 
crease.  However,  once  the  overlying  water  column  concentration  reaches  that  of  the 
sediment  pore  water,  there  is  no  longer  a  gradient  and  the  sediment  flux  of  iron  from 
the  sediment  to  the  overlying  water  ceases.  At  this  point  the  hypolimnion  concentra¬ 
tion  remains  approximately  constant.  It  can  be  seen  from  Fig. 8. 4-8.5  that  this  occurs 
near  the  end  of  May  and  beginning  of  June  when  the  bottom  Fe^'^  concentration 
reaches  about  ~  0.8  mg/L.  At  this  point  two  things  happen.  First,  the  depth  of  the 
thermocline  increases,  beginning  from  a  depth  of  13  m  from  the  water  surface  and 
climbing  to  10  m  and  then  to  8  m  from  the  water  surface  over  a  30  day  period.  Sec¬ 
ondly,  although  the  concentration  at  the  bottom  will  remain  constant,  the  continued 
settling  and  subsequent  reduction  of  particulate  iron  will  continue.  This  source  from 
the  overl5dng  oxic  water  eventually  will  cause  a  peak  to  occur  in  dissolved  iron,  not 
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Fig.  8.5  Calibration  results  for  Onondaga  Lake.  June— August,  1986;  concentrations  versus 
depth.  Columns  are  each  month.  Rows  are  temperature,  dissolved  oxygen,  dissolved  iron, 
particulate  iron. 
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at  the  bottom,  but  at  a  location  just  below  the  thermocline.  This  behavior  is  clearly 
shown  in  the  months  of  June  through  September. 

The  model  overestimates  the  amount  of  slightly  in  August,  but  more  notice¬ 
ably  in  September  and  October  (Fig.8.6).  One  reason  for  this  may  be  the  formation 
and  subsequent  precipitation  of  iron  sulfide,  FeS.  This  mechanism  is  discussed  in 
more  detail  in  Section  8.2.6.  Note  that  the  Fe^-^  data  for  November  1986  are  only  for 
one  day,  November  3. 

A  temporal  summary  of  the  calibration  for  1986  is  shown  in  Figure  8.7.  The  top  plot 
is  a  temporal  plot  of  epilimnion  and  hypolimnion  temperatures.  The  depth  of  the  ther¬ 
mocline  was  estimated  at  8  m.  The  symbols  represent  the  average  observed  epilimnion 
and  hypolimnion  concentrations.  The  solid  line  represent  the  model  monthly  averages 
with  the  shading  representing  the  range.  As  can  be  seen,  the  model  is  reasonably  well 
calibrated  to  the  observed  temperature. 

Figme  8.7  also  contains  a  similar  plot  for  dissolved  oxygen.  Note  the  persistence  of 
anoxia  in  the  lake  for  1986.  This  was  not  originally  reproduced  by  the  model  and  so 
had  to  be  specified. 

The  bottom  plot  in  Figure  8.7  shows  Fe^+  concentrations.  Since  Fe^+  is  rapidly  ox¬ 
idized  in  aerobic  waters,  it  is  only  present  during  anoxic  conditions  in  the  h3rpolimnion. 
The  model  reproduces  both  the  average  and  the  range  of  the  data.  The  exception  is 
in  September  when  both  the  model  range  and  average  overpredict  the  data,  and  in 
October  when  the  model  average  is  close  but  the  range  is  too  large.  As  mentioned 
previously,  this  could  be  the  effect  of  precipitation  of  FeS,  and  will  be  discussed  in 
Section  8.2.6. 

In  1987  no  iron  data  were  collected.  Figs.  8.8-8.9  show  the  calibration  results  for 
1988  for  the  months  of  May  through  October.  It  is  interesting  to  note  that  the  data 
for  both  1988  and  1989  exhibit  much  more  short  term  variation  than  in  1986.  It  is 
not  clear  why  this  happens,  and  although  the  model  fits  the  monthly  concentration 
averages,  it  cannot  reproduce  the  large  variations.  The  observed  dissolved  oxygen  data 
does  not  suggest  any  reaeration  events  occurring,  which  would  cause  large  variations 
in  Fe2+. 

Figure  8.10  is  a  temporal  plot  of  the  1988  calibration  for  the  epilimnion  and  hy¬ 
polimnion.  Again  the  large  ranges  in  are  noticeable. 

Iron  data  were  collected  for  only  three  months:  Jime  through  August  in  1989.  The 
calibration  is  shown  in  Figure  8.11,  and  the  temporal  summary  for  the  epilimnion  and 
hypolimnion  are  presented  in  Figme  8.12.  The  model  demonstrates  good  agreement 
with  the  data. 

It  is  of  interest  to  determine  the  extent  to  which  the  model  predictions  are  dependent 
on  an  accmrate  representation  of  the  dissolved  oxygen  concentrations.  This  is  examined 
in  the  next  three  figures  8.13—8.15  which  show  temporal  plots  of  the  iron  calibration 
using  the  dissolved  oxygen  concentrations  calculated  by  the  model.  Note  the  marked 
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Fig.  8.6  Calibration  results  for  Onondaga  Lake.  September-Noverber  1986;  concentrations 
versus  depth.  Columns  are  each  month.  Rows  are  temperature,  dissolved  oxygen,  dissolved 
iron,  particulate  iron. 


Fig.  8.11  Calibration  results  for  Onondaga  Lake.  June-August,  1989 


142  APPLICATION  OF  IRON  FLUX  MODEL  TO  ONONDAGA  LAKE 


effect  on  the  1986  calibration.  The  onset  of  iron  generation  from  the  sediment  is 
delayed  until  the  dissolved  oxygen  is  zero  at  the  sediment-water  interface.  The  model 
predicted  the  period  of  anoxia  in  1988  and  1989  and  so  there  is  little  difference  between 
the  results  where  the  period  of  anoxia  is  specified  (Figs.  8.7,  8.10  and  8.12). 

8.2.6  Iron  Sulfide  Precipitation 

The  suggestion  was  made  above  that  the  decline  in  dissolved  iron  in  the  latter  portions 
of  1986  was  due  to  the  precipitation  of  iron  sulfide.  The  sediment  model  computes 
sulfide  fluxes  and  the  water  column  model  tracks  the  oxygen  equivalents  of  dissolved 
sulfide.  It  is  possible,  therefore,  to  see  if  the  sulfide  concentrations  are  increasing  to  a 
point  where  precipitation  is  possible. 

Figures  8.16,  8.17,  and  8.18  are  plots  of  sulfide  and  dissolved  iron  for  the  months 
of  May  through  October  for  1986,  1988  and  1989  respectively.  For  1986  (Fig.  8.16) 
the  buildup  is  reproduced  until  October  when  premature  destratification  in  the  model 
reduces  the  sulfide  concentration.  The  computed  sulfide  concentrations  for  1988  (Fig. 
8.17)  increase  more  rapidly  than  is  observed,  and  the  data  for  the  two  months  in  1989 
suggest  the  same  overprediction.  The  magnitude  of  the  sediment  flux  of  sulfide  during 
anoxia  is  influenced  by  the  magnitude  of  the  sulfide  partition  coefficients  tth.si  and 
^H2S,2-  These  are  empirical  parameters  that  were  calibrated  using  the  Chesapeake 
Bay  data  set  [Di  Toro  and  Fitzpatrick,  1993].  Since  these  parameters  have  not  been 
calibrated  for  Onondaga  Lake,  the  results  are  encouraging. 

The  sulfide  and  iron  data  can  be  used  to  make  a  qualitative  analysis  of  the  likelihood 
that  sulfide  precipitation  is  the  cause  of  the  model  overprediction  of  dissolved  iron  and 
sulfide.  The  excessive  iron  concentrations  are  occurring  primarily  in  September  and 
October  of  1986  (Fig.  8.16).  There  is  no  indication  that  the  sulfide  concentrations 
are  increasing  rapidly  enough  to  initiate  precipitation  of  iron  sulfide.  It  is  possible,  of 
course,  that  FeS  (s)  is  being  formed  continuously  and  that  its  addition  as  a  mechanism 
would  improve  the  model’s  predictive  ability.  However,  a  qualitative  examination  of 
the  model-data  discrepancies  does  not  suggest  that  this  is  a  first  order  effect. 


8.2.7  Sediment  Concentrations 

Some  limited  data  are  available  for  both  pore  water  dissolved  iron  and  sediment  iron 
concentrations  in  Onondaga  Lake.  The  data  are  summarized  in  Table  8.10,  which  also 
shows  the  model’s  calculated  annual  average  concentrations.  Sediment  iron  data  are 
available  for  1986  [Yin  and  Johnson,  1984].  Pore  water  dissolved  iron  data  are  available 
for  August  and  November  1991  [Cornwall,  1993].  The  computed  and  observed  pore 
water  concentrations  are  in  the  same  order  of  magnitude  and  the  range  is  similar. 
The  solid  phase  comparison  is  also  reasonable  when  the  observed  concentrations  are 
corrected  for  the  reactive  fraction,  as  shown. 
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Fig.  8.17  Iron  and  sulfide  calibration  for  Onondaga  Lake.  August-October,  1988 


Fig.  8.18  Iron  and  sulfide  calibration  for  Onondaga  Lake.  August-October,  1989 
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Table  8.10  Compaxison  of  observed  and  computed  sediment  and  pore  water  iron  concentra¬ 
tions:  annual  average  (range) 


Phase 

Model 

Data 

Pore  water  (mg  Fe/L) 

0.88 

0.58 

(0.02-2.78) 

(0.13-3.34) 

Solid  phase  (mg  Fe/g) 

4.0 

5.6“ 

(3.0-5.5) 

(3.6-7.8)“ 

“Reactive  iron  is  estimated  as  40%  of  total  iron  consistent  with  the  reactive  fraction  assumption  used 
to  estimate  the  loading 


8.2.8  Summary 

The  iron  sediment  and  water  column  model  formulated  in  this  chapter  can  reproduce 
the  available  observations  of  water  column  and  sediment  concentrations  in  Onondaga 
Lake  for  three  years  with  reasonable  fidelity.  The  model  is  initialized  to  a  periodic 
steady  state  and  then  run  for  5  years,  beginning  in  1985  and  ending  in  1989.  The 
same  oxidation  and  reduction  kinetic  rate  constants  are  used  in  both  the  water  col¬ 
umn  and  sediment,  consistent  with  the  procedure  followed  for  the  manganese  model. 
The  seasonal  development  of  the  vertical  distribution  of  iron  in  the  water  column  is 
correctly  reproduced.  An  initial  increase  occurs  at  the  sediment-water  interface  at  the 
beginning  of  anoxia  (March-May  1986),  followed  by  an  almost  uniform  concentration 
in  the  hypolimnion  (Jtme),  and  then  the  development  of  a  maximum  concentration  at 
the  thermocline  (July-Sept.),  and  the  subsequent  decline  at  overturn. 

The  predicted  sulfide  concentrations  are  also  in  reasonable  agreement  with  observa¬ 
tions  without  any  further  calibration  of  the  partitioning  parameters  for  sulfide.  There 
does  not  appear  to  be  any  strong  evidence  for  the  precipitation  of  iron  sulfide  if  the 
criteria  used  is  the  discrepancy  between  model  and  data.  This  does  not  rule  out  the 
possibihty  and,  therefore,  the  addition  of  this  mechanism  would  be  a  logical  next  step. 
It  would  also  remove  the  need  to  specify  the  partition  coefiicients  for  sulfide  and  iron, 
thereby  removing  a  degree  of  freedom  from  the  model. 
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Application  to  Croton  Reservoir 


9.1  INTRODUCTION 

The  eutrophication  model  for  the  Croton  Reservoir  was  constructed  as  a  part  of  a 
study  of  the  water  supply  system  for  the  City  of  New  York  [Metcalf  et  al.,  1995].  The 
iron  and  manganese  models  developed  in  this  study  were  incorporated  in  that  model 
and  the  results  are  presented  below. 

9.1.1  Physical  Characteristics 

The  Croton  Reservoir  is  located  approximately  40  miles  north-east  of  New  York  City, 
in  Westchester  and  Putnam  Counties,  in  lower  New  York  State.  The  reservoir  rep¬ 
resents  part  of  the  drinking  water  supply  system  for  New  York  City  and  supplies  on 
average  about  140  million  gallons  per  day  of  drinking  water  to  the  city.  The  Croton 
Reservoir  was  created  by  the  impoimdment  of  the  Croton  River  by  the  New  Croton 
Dam  At  the  northwest  end  of  the  reservoir  is  the  Muscoot  Dam  which  delineates 
the  Muscoot  Reservoir.  Approximately  4  km  upstream  from  the  Croton  Dam,  in  the 
Croton  Reservoir  (sometimes  called  the  New  Croton  Reservoir),  is  a  submerged  dam 
that  formed  the  original  Croton  Reservoir  at  this  site. 

The  reservoir  is  oriented  along  a  southeast-northwest  axis,  with  the  New  Croton 
Dam  located  at  the  southeast  end.  A  bathymetric  map  of  the  Croton  Reservoir, 
which  also  outlines  the  model  surface  segmentation,  and  sampling  sites,  is  shown  in 
Fig.9.1.  The  length  of  the  Croton  Reservoir  is  approximately  14.6  km.  The  maximum 
depth  when  full  is  36  m.  The  surface  area  when  full  is  92  x  10®  ft^  and  the  volume 
when  full  is  4.14  x  10®  ft®.  The  average  residence  time  of  the  Croton  Reservoir  is  91 
days.  The  direct  drainage  area  is  56  square  miles. 

At  the  Muscoot  Dam  the  depth  of  the  Croton  Reservoir  is  about  10  m  depending 
on  the  water  level  in  the  reservoir.  The  depth  increases  to  16  m  at  sampling  site  6. 
Upstream  of  the  submerged  dam,  at  sampling  site  5,  the  depth  is  20  m  .  At  site  4,  just 
after  the  submerged  dam,  the  depth  is  24  m.  At  the  New  Croton  Dam  the  maximum 
depth  of  the  reservoir  is  40  m. 

9. 1.1.1  Tributaries  Most  of  the  inflow  to  the  Croton  Reservoir,  85-90%,  comes  from 
the  adjoining  Muscoot  Reservoir.  The  rest  is  inflow  from  tributaries,  groimdwater 
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Fig.  9.1  New  Croton  Resevoir  segmentation. 


recharge,  and  surface  runoff.  The  main  tributaries  are  the  Kisco  River  and  Himter 
Brook.  The  largest  inflows  occur  in  the  months  January  through  April.  The  smallest 
inflows  occur  in  the  months  June  through  September.  In  order  to  maintain  adequate 
stream  conditions  in  the  Croton  River  downstream  of  the  reservoir,  a  minimum  outflow 
of  5.5  MGD  is  required  to  be  maintained.  A  summary  of  the  inflow  for  1993  is  shown 
in  Table  9.1. 

9.1.2  Loadings 

The  nutrient  loadings  to  the  Croton  reservoir  are  from  the  following  sources:  runoff, 
wastewater  treatment  plants,  septic  systems,  atmospheric  deposition,  and  the  Muscoot 
Reservoir  inflow.  Phosphorus  loading  is  of  most  interest  in  the  Croton  and  Muscoot 
Reservoirs  because  it  is  the  limiting  nutrient.  Phosphorus  loading  from  the  Muscoot 
Reservoir  for  1993  was  estimated  to  be  42  kg/day.  The  two  tributaries,  the  Kisco 
River  and  Hunter  Brook,  together  contribute  an  estimated  additional  6  kg/day  to  the 
Croton  Reservoir.  Atmospheric  deposition  is  estimated  to  contribute  0.6  kg/day  to 
the  Croton  Reservoir. 

9.1.3  Water  Quality 

The  Croton  reservoir  is  mesotrophic  with  intermediate  levels  of  chlorophyll.  It  is  highly 
stratifled  in  the  summer  months  from  June  to  September  with  temperature  differences 
exceeding  18°C  .  In  1993  hypolimnion  hypoxia  and  anoxia  at  Stations  5  and  6  lasted 
from  June  through  November.  Epilimnion  dissolved  oxygen  concentrations  at  these 
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Table  9.1  Inflow  To  Croton  Reservoir  in  1993  (MGD)“ 


Month 

Muscoot  inflow 

Rimoff 

Other  tributaries 
and  groimdwater 

January 

472 

3 

52 

February 

243 

8 

53 

March 

665 

22 

50 

April 

796 

8 

60 

May 

249 

3 

31 

Jme 

153 

7 

26 

July 

126 

0 

24 

August 

149 

23 

27 

September 

144 

11 

36 

October 

198 

10 

39 

November 

184 

14 

34 

December 

283 

26 

53 

Average 

305 

11 

39 

“[Metcalf  et  al.,  1995] 


stations  were  above  8  mg/L  all  year.  Table  9.2  shows  some  relevant  water  quality  data 
on  a  seasonal  basis  for  Stations  5  and  6. 


Table  9.2  Croton  Reservoir  water  quality  overview  for  Stations  5  and  6  for  1993 


Variable 

Summer 

Winter 

Epihmnion  Temp  (°C) 

24 

6 

Hypolimnion  Temp  (°C) 

7 

6 

Epilimnion  DO  (mg/L) 

9.0 

11 

Hypohmnion  DO  (mg/L) 

0.0 

11 

Chlorophyll  a  (^tg/L) 

6.0 

3.0 

Epilimnion  Chlorides  (mg/L) 

34 

34 

Hypolimnion  Chlorides  (mg/L) 

34 

34 

Epihmnion  pH 

8.5 

7.2 

Hypolimnion  pH 

6.7 

7.2 

9.1.4  Data  Sets 

Extensive  water  quality  monitoring  of  the  Croton  and  Muscoot  Reservoirs  has  been 
carried  out  by  the  New  York  City  Department  of  Environmental  Conservation  since 
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1976.  The  parameters  that  have  been  measured  include  both  iron  and  manganese. 
The  location  of  the  sampling  sites  is  shown  on  Figure  9.1. 

9. 1.4.1  Iron  and  Manganese  Two  sites  in  particular  have  useful  iron  and  manganese 
data,  Stations  5  and  6.  Station  5  is  located  just  upstream  of  the  submerged  original 
Croton  Dam.  Station  6  is  located  approximately  4.5  km  downstream  of  the  Muscoot 
Reservoir. 

For  1993,  the  data  at  Stations  5  and  6  was  collected  monthly  from  April  to  Decem¬ 
ber.  While  many  parameters  such  as  dissolved  oxygen  and  temperature  were  measured 
at  1  m  depths,  the  sampling  protocol  for  iron  and  manganese  was  different.  Two  grab 
samples  were  taken  at  a  site:  one  sample  above  the  thermocline  and  one  sample  below, 
near  the  reservoir  bottom.  Only  dissolved  iron  and  mangjmese  were  analyzed. 

Monthly  measurements  of  total  iron  and  total  manganese  were  made  at  the  tribu¬ 
tary  inflows  to  the  Muscoot  and  Croton  Reservoirs.  These  data  are  used  to  generate 
the  daily  iron  and  manganese  loads  to  the  system.  As  all  tributaries  are  aerobic,  the 
influent  iron  and  manganese  loads  are  assumed  to  be  particulate  concentrations.  At 
present  there  are  no  sediment  iron  or  manganese  data  available  at  the  Croton  Reser¬ 
voir. 

9.2  MODEL  SETUP 

For  the  model  of  the  Croton  and  Muscoot  Reservoirs  the  water  quality  model,  RCA, 
has  34  state  variables.  The  first  32  variables  are  the  same  as  used  in  the  Onondaga 
Lake  model  and  are  listed  in  Table  8.5.  An  additional  2  variables  were  added  for 
manganese.  Like  iron,  manganese  is  modeled  as  dissolved  Mn  {II)  and  particulate 
MnO^  (s)  manganese  as  in  Chapter  7.  Both  the  manganese  and  iron  sediment  model, 
discussed  in  the  previous  Chapter  8  are  included.  This  application  provides  a  test  of 
both  models. 

9.2.1  Manganese  Flux  Model  Equations 

The  manganese  flux  model  is  virtually  identical  to  that  applied  to  the  MERL  mesocosm 
data  set  (Chapter  7).  The  only  change  is  the  addition  of  the  reduction  reaction  of 
Mn02  (5)  to  Mn  {II)  if  the  hypolimnion  becomes  anoxic.  A  schematic  is  presented  in 
Fig.7.2.  The  equations  for  the  aerobic  (1)  and  anaerobic  (2)  sediment  layers  are: 

Layer  1  Mn{II): 


dMn{l) 


-s  (/diMn(l)  -  Mn  (0))  -b  KLi2{fd2Mn{2)  -  /diMn(l)) 
+wi2{fp2Mn{2)  -  fpiMn{l))  -W2Mn{l) 
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F/^.  P.2  Schematic  diagram  of  the  manganese  flux  model 


(0, /«Mn(l) 

+Mn{2)H+  (9.1a) 

Layer  2  Mn{II): 

dMn(2) 

=  -iirLi2(/d2Mn(2)-/diMn(l)) 

—‘Wi2ifp2Mn{2)  -  fpiMn{l))  —  ty2  (Mn(2)  —  Mn(l)) 
+K'Mn,2Mn02{2)  —  Mn{2)  ^  +  Mn(l)H^  (9.1b) 


Layer  1  Mn02  (s): 

dMnC>2(l) 

■^1 - - 


k2 


O2  (0)  /<iiMn(l)  -  W2Mn02{l) 


s 
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+wi2{Mn02{2)  -  MnOaCl))  +wiMn02  (0) 

+Mn02{2)H^  -  Mn02(l)  (Hi  +  (9.1c) 


Layer  2  Mn02  (s): 

=  -KMnfiMn02{2) 

-wi2(Mn02{2)  -  Mn02(l))  -  W2  {Mn02{2)  -  Mn02(l)) 
-Mn02{2)  (h2  +  i?i+)  +  Mn02il)Hr  (9.1d) 

where  Mn02  (0)  and  Mn  (0)  refer  to  the  particulate  and  dissolved  manganese  concen¬ 
trations  in  the  bottom  most  water  column  segment,  designated  as  segment  0,  which 
overla}^  the  sediment. 

The  water  column  equations  for  manganese  for  this  segment  are 
Layer  0  Mn{II): 


=  s{fdiMn{l)-Mn{0))-KMn,o[lO‘^^”°-'^^)02(0)Mn{0) 


+Eo-i  (Mn(-l)-Mn(O)) 


(9.1e) 


Layer  0  Mn02  (s): 

^jq2(pHo-7)^  Mn(0)  -I-  wi  {Mn02  (-1)  -  MnC>2  (0)) 


-kMnfiHo 

KM,Mn,0  +  [O2  (U}J 
+E0-1  (Mn02  (-1)  -  Mn02  (0)) 


(9.1f) 


where  kMnfi  is  the  rate  at  which  Mn02{s)  is  reduced  to  Mn{II);  KM,Mn,o  is  the 
half  saturation  constant  for  oxygen  for  this  reaction;  and  jBo,-i  represents  the  vertical 
miving  coefficient  between  the  bottom  layer  and  the  layer  immediately  above.  The 
concentrations  in  the  water  column  are  numbered  in  reverse  order  to  preserve  the 
notation  used  previously.  Thus  layer  -1  refers  to  the  layer  just  above  layer  0  in 
the  water  colunm.  The  remaining  portions  of  the  sediment  model — the  nutrient  and 
oxygen  components — are  identical  to  those  used  above. 
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9.2.2  Model  Segmentation  for  Croton  and  Muscoot  Reservoirs 

The  surface  segmentation  for  the  Croton  and  Muscoot  model  is  shown  in  Figure  9.1. 
The  model  has  55  surface  segments,  32  in  the  Croton  Reservoir  and  23  in  the  Muscoot. 
Vertically  the  model  has  10  layers  and  uses  a  stretching  coordinate  system  to  adjust 
the  layer  depth  to  fit  the  varying  bathymetry. 

9.2.3  Model  Calculations  and  Results 

9.2.3.1  Dissolved  Oxygen  It  is  critical  that  the  timing  of  the  onset  of  anoxia  be 
calculated  correctly  since  this  event  triggers  the  release  of  iron  and  manganese  from 
the  sediment.  As  with  the  Onondaga  Lake  application,  it  was  decided  that  actual 
dissolved  oxygen  data  would  be  used  to  specify  when  anoxia  was  present.  This  allows 
us  to  examine  how  the  iron  and  manganese  model  would  perform  if  the  eutrophication- 
dissolved  oxygen  model  was  perfectly  calibrated  to  the  dissolved  oxygen  data. 

9.2.3.2  Calibration  Scenario  The  model  was  run  for  1993,  beginning  January  1  to 
December  31,  a  total  of  365  days.  In  order  to  ensure  that  the  sediment  had  reached 
a  steady  state  the  model  was  cycled  for  several  years  before  using  the  results  as  the 
calibration  output. 

The  values  of  the  constants  used  for  the  nutrients  and  oxygen  sediment  model  are 
given  below  in  Table  9.5  The  values  of  the  constants  used  for  the  iron  and  manganese 
models  in  the  water  column  and  sediment  are  given  below  in  Tables  9.3  and  9.4.  The 
nutrients,  oxygen,  and  iron  model  values  are  identical  to  those  used  for  Onondaga 
Lake  with  one  exception  noted  below. 

9.2.4  Calibration  Results 

As  mentioned  previously,  the  iron  and  manganese  data  for  Croton  Reservoir  consisted 
of  samples  at  two  depths,  which  were  collected  approximately  once  a  month.  One 
sample  was  taken  in  the  epihmnion,  which  was  always  oxic  and  hence  yielded  negligible 
dissolved  iron  and  manganese.  The  hypolimnion  sample  was  taken  near  the  bottom. 
Therefore  no  comparison  can  be  made  to  the  computed  vertical  profile. 

Figure  9.3  shows  the  calibration  results  for  site  5.  The  surface  and  bottom  model 
(lines)  and  data(symbols)  are  shown  for  temperature  and  dissolved  oxygen.  The  bot¬ 
tom  dissolved  iron  and  dissolved  manganese  concentrations  are  compared  to  the  model 
computations  in  the  bottom-most  layer  of  the  model.  The  extended  period  of  anoxia 
is  clearly  evident,  from  early  June  until  early  November.  For  these  simulations  the 
observed  time  for  the  both  onset  and  end  of  anoxic  was  used  to  correct  the  bottom 
layer  dissolved  oxygen  computed  by  the  model. 
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Tsbie  9.3  Constants  for  the  iron  flux  model  -  Croton  Reservoir 


Water  Column 

Name 

Value 

Unit 

Description 

kpe,! 

10.0 

day"^ 

Oxidation  rate 

^Fe,l 

1.08 

- 

Temperature  coefficient 

kFe,2 

4.0 

day“^ 

Reduction  rate 

KM,Fe,2 

0.2 

mg  O2/L 

Oxidation  half  saturation  constant 

Sediment 


ATTfe,! 

15.0 

L/kg 

Aerobic  layer  partition  coefficient 

T^Fe,2 

100.0 

L/kg 

Anaerobic  layer  partition  coefficient 

02,Crit,Fe 

2.0 

mg/L 

Partition  coefficient  half  saturation  constant 

kpe,! 

10.0 

day~^ 

Aerobic  layer  oxidation  rate 

Ope,! 

1.08 

Temperature  coefficient 

kFe,2 

4.0 

day“^ 

Anaerobic  layer  reduction  rate 

&Fe,2 

1.08 

Temperature  coefficient 

The  iron  data  shows  very  large  concentrations  during  anoxia,  between  4  and  6 
mg/L.  The  results  for  site  6,  which  are  shown  in  Figure  9.4,  show  markedly  lower  iron 
concentrations,  ~  1  mg/L  for  July  and  August,  but  ~  4  mg/L  for  early  September. 
By  comparison,  Onondaga  Lake  bottom  iron  concentrations  were  <  1  mg/L. 

In  order  to  reproduce  this  diflference,  it  is  necessary  to  imderstand  what  controls  the 
maximum  hypolimnion  iron  concentration.  The  following  observation  is  useful.  Within 
~  14  days  from  the  onset  of  anoxia,  the  bottom  concentration  will  be  approximately 
the  same  as  the  pore  water  iron  concentration.  For  Onondaga  Lake  this  was  ~  1 
mg/L.  Hence,  in  order  for  the  model  to  calculate  4  mg/L  in  the  bottom  layer  of  the 
Croton  reservoir,  the  sediment  partition  coefficient  for  iron  must  be  reduced  in  order 
to  increase  the  pore  water  iron  concentration.  The  results,  which  are  shown  in  Figures 
9.3  and  9.4,  was  achieved  by  changing  the  anaerobic  layer  partition  coefficient  from 
7fFe,2  =  3000  L/kg,  which  was  used  for  Onondaga  Lake,  to  7rjre,2  =  100  L/kg. 

It  is  not  clear  why  the  pore  water  iron  concentration  in  the  Croton  reservoir  is  larger, 
and  the  partition  coefficient  is  lower,  than  in  Onondaga  Lake.  This  is  the  principle 
drawback  associated  with  using  an  empirical  partition  coefficient  in  the  sediment  flux 
model  to  determine  the  pore  water-solid  phase  speciation  .  There  is  no  guarantee  that 
it,  or  any  other  of  the  partition  coefficients  used  in  the  model,  are  transferable  from 
one  application  to  another. 
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Table  9.4  Constants  for  the  manganese  flux  model  -  Croton  Reservoir 

Water  Column 


Name 

Value 

Unit 

Description 

^Mn,l 

10.0 

day“^ 

Oxidation  rate 

1.08 

- 

Temperature  coefficient 

kMn,2 

4.0 

day“^ 

Reduction  rate 

KM,Mn,2 

0.2 

mg  O2/L 

Oxidation  half  saturation  constant 

Sediment 

AtTmti,! 

200.0 

L/kg 

Aerobic  layer  partition  coefficient 

T^Mn,2 

8000. 

L/kg 

Anaerobic  layer  partition  coefficient 

02jCritjMn 

2.0 

mg/L 

Partition  coefficient  half  saturation  constant 

^Mn,! 

0.01 

day“^ 

Aerobic  layer  oxidation  rate 

^Mn,l 

1.08 

- 

Temperature  coefficient 

^Mn,2 

4.0 

day“^ 

Anaerobic  layer  reduction  rate 

&Mn,2 

1.08 

Temperature  coefficient 

The  manganese  concentrations,  by  contrast,  are  about  the  same  at  site  5  and  site 
6.  As  can  be  seen,  the  model  produces  a  good  fit  to  the  data  with  no  need  for  further 
adjustments. 

9.3  SUMMARY 

The  Croton  reservoir  was  used  as  a  field  test  of  the  iron  and  manganese  model  because 
data  are  available  and  the  eutrophication  model  had  been  developed  previously.  The 
model  was  run  for  1993,  after  equilibrating  the  sediment.  The  calibration  for  iron 
involved  changing  the  anaerobic  layer  sediment  partition  coefficient  from  the  value 
used  for  Onondaga  Lake  in  order  to  compute  the  higher  iron  concentrations  observed 
in  the  Croton  reservoir.  This  difference  in  dissolved  iron  between  the  two  lakes  is 
imexplained  at  present.  The  calibrated  manganese  model  was  able  to  reproduce  the 
observed  data  without  any  further  adjustments. 
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Table  9.5  Sediment  Model  Parameters  for  the  Croton  Reservoir 


mi 

m2 

W2 

H2 

(kg/L) 

(kg/L) 

(m/d) 

(m) 

0.5 

6.85  X  10-® 

0.10 

Dd 

Dp 

Gdj. 

Ddo 

^£><10 

(mVd) 

- 

(mVd) 

- 

(mVd) 

- 

5.0E-03 

1.08 

1.2  X  lO"'* 

1.15 

1.0  X  lO-'* 

1.08 

Km,nh4 

Km,nh4,02 

(m/d) 

(L/kg) 

(mg  N/L) 

- 

(mg  O2/L) 

0.2 

1.00 

1.123 

0.728 

1.125 

0.74 

I<'N03,2 

OnOs 

(m/d) 

(m/d) 

- 

0.15 

0.25 

1.080 

f^H2S,dl 

I^H2S,pl 

GH2S 

Km,H2S,02 

(m/d) 

(m/d) 

(L/kg) 

(L/kg) 

- 

(mg  O2/L) 

0.2 

0.4 

100.0 

100.0 

1.15 

4.0 

ksi 

[Si].at 

Arrsi,! 

7tSi,2 

Osi 

Km,ps\ 

(/d) 

(mg  Si/L) 

- 

(L/kg) 

- 

(mg  Si/L) 

0.5 

30. 

6.0 

175.0 

1.100 

5.0  X  10^ 

A7rpo4,i 

7rpo4,2 

^2  (0)^_po4 

^ M,Dp,02 

O2  (0).,si 

— 

(L/kg) 

(mg  O2/L) 

(mg  O2/L) 

(mg  O2/L) 

65.0 

225.0 

2.0 

4.0 

1.0 

65.0 
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